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This review summarizes the program in the physics of the top quark being pursued at 
Fermilab's Tevatron proton-antiproton collider at a center of mass energy of 1.96 TeV. 
More than a decade after the discovery of the top quark at the two collider detectors 
CDF and DO, the Tevatron has been the only accelerator to produce top quarks and to 
study them directly. 

The Tevatron's increased luminosity and center of mass energy offer the possibility 
to scrutinize the properties of this heaviest fundamental particle through new measure- 
ments that were not feasible before, such as the first evidence for electroweak production 
of top quarks and the resulting direct constraints on the involved couplings. Better mea- 
surements of top quark properties provide more stringent tests of predictions from the 
standard model of elementary particle physics. In particular, the improvement in mea- 
surements of the mass of the top quark, with the latest uncertainty of 0.7% marking 
the most precisely measured quark mass to date, further constrains the prediction of the 
mass of the still to be discovered Higgs boson. 
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1. Introduction 




The existence of a third and most massive generation of fundamental fermions was 
unveiled in 1975 with the discovery of the t lepton at SLAC-LBL [1]. In 1977, 
the discovery of the bottom quark [2] at Fermilab extended the knowledge of a 
third generation into the quark sector and immediately raised the question of the 
existence of the top quark as the weak isospin partner of the bottom quark. 

To remain self consistent, the standard model (SM) of elementary particle 
physics required the existence of the top quark, and electroweak precision mea- 
surements offered increasingly precise predictions of properties such as its mass. 
The top quark's large mass prevented its discovery for almost two decades, but by 
1994 it was indirectly constrained to be 178 ± 11 GeV/c^ [3]. After mounting 
experimental evidence [4-9], the top quark (t) was finally discovered in 1995 at 
Fermilab by the CDF and DO collaborations [10,11] in the mass range predicted by 
the standard model. The completion of the quark sector once again demonstrated 
the enormous predictive power of the SM. 

By now, the mass of the top quark is measured to be 172.4 ± 1.2 GeV/c^ [12], 
marking the most precisely measured quark mass and the most massive fundamen- 
tal particle known to date. The consequent lifetime of the top quark in the SM of 
w 5 • 10^^^ s is extremely short, suggesting that it decays before hadronizing. This 
makes it the only quark that does not form bound states, allowing the study of 
an essentially bare quark with properties such as spin undisturbed by hadroniza- 
tion [13]. 

The measurement of top quark pair {tt ) production probes our understanding 
of the strong interaction and predictions from perturbative QCD, while the decay 
of top quarks and the production of single top quarks reflect the electroweak inter- 
action. Measuring other properties of the top quark, such as its electric charge, the 
helicity of the W boson in t ^ Wb decay, the branching fraction B{t ^ Wb), etc., 
and comparing these with predictions of the SM is a very powerful tool in searching 
for new physics beyond the standard model. 

The top quark can also be used to constrain the mass range of the last yet to 
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be observed particle of the standard model, the Higgs boson, because their masses 
and the mass of the W boson are linked through radiative corrections [14]. The 
Higgs boson is a manifestation of the Higgs mechanism [15-20], implemented in 
the standard model to provide the needed breaking of the electroweak symmetry to 
which the top quark may be intimately connected because of its large mass. 

Because of its fairly recent discovery, the top quark's properties have not yet 
been explored with the same scrutiny as those of the lighter quarks. However, in the 
ongoing data taking at Fermilab's Tevatron proton-antiproton collider, an integrated 
luminosity of more than 4 fb~^ has already been recorded by each of the collider 
experiments CDF and DO, corresponding to an increase of about a factor seventy 
relative to the data that was available for the discovery of the top quark. The new 
data can be used to refine previous measurements to higher precision that starts 
to become limited by systematic rather than statistical uncertainties. In addition, 
measurements that have never been performed become feasible, such as the first 
evidence for electroweak production of single top quarks and the consequent first 
direct measurement of the CKM matrix element \Vtb\, recently published by DO 
[21,22] and CDF [23]. 

This article is intended to provide an overview of the current status of the top 
quark physics program pursued at the Tevatron. Results available until the LHC 
startup in September 2008 have been included, utilizing samples of data of up to 
2.8 fb~^ in integrated luminosity. Previous reviews of the top quark are available in 
Refs. [24-28]. The outline of this article is as follows: The second chapter provides a 
brief introduction to the standard model, with emphasis on the special role played 
by the top quark. Chapter 3 describes production and decay modes for top quarks 
in the framework of the standard model. Chapter 4 outlines the experimental setup 
used for the measurements described in the following sections. Chapter 5 presents 
studies of the production of top quarks, including measurements of cross section 
that form the basis for other measurements of top quark characteristics. Chapter 
6 elaborates on the different results for top quark decay properties, followed in 
Chapter 7 by a discussion of measurements of fundamental attributes of the top 
quark, such as its charge and mass. The final chapter (8) contains a brief summary 
of the achievements to date. 

2. The Standard Model and the Top Quark 
2.1. A brief overview of the standard model 

The standard model of elementary particle physics describes very successfully the 
interactions of the known fundamental spin J = \ fermion constituents of matter 
through the exchange of spin J = 1 gauge bosons. 

As shown in Table [l] both quarks and leptons occur in pairs, differing by one 
unit of electric charge e, and are replicated in three generations that have a strong 
hierarchy in mass. The fermion masses span at least 11 orders of magnitude, with 
the top quark being by far the heaviest fundamental particle, which may therefore 
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Table 1: The known fundamental fermions and their masses [12,29]. 



fermion 


electric 


generation 


type 


charge [e] 


1. 


2. 


3. 


quarks 


^3 


up (u) 
1.5 - 3.3 MeV/c2 


charm (c) 
1.27 t^il GeV/c2 


top (t) 
172.4 ± 1.2 GeV/c2 


1 
3 


down (d) 
3.5 - 6.0 MeV/c2 


strange (s) 
104 til MeV/c2 


bottom (6) 
4.20 torn GeV/c2 







< 2 eV/c2 


< 0.19 MeV/c2 


< 18.2 MeV/c2 


leptons 




(95% C.L.) 


(90% C.L.) 


(95% C.L.) 




-1 


e 

0.511 MeV/c^ 


M 

105.658 MeV/c^ 


T 

1777 MeV/c2 



provide further insights into the process of mass generation. The origin of this 
breaking of the flavor symmetry and the consequent mass hierarchy is still not 
understood but can be accommodated in the standard model as shown below. 

The forces among the fundamental fermions are mediated by the exchange of the 
gauge bosons of the corresponding quantized gauge fields, as listed in Table [2l The 
gravitational force is not included in the framework of the standard model, and will 
not be considered, as its strength is small compared to that of the other interactions 
among the fundamental fermions at energy scales considered in this article. 

The standard model is a quantum field theory based on the local gauge symme- 
tries SU{3)qcd X SU{2)l X U{1)y- The theory of the strong interaction, coupling 
three different color charges ("red", "green" and "blue") carried by the quarks 
and the eight massless gauge bosons (gluons), is called Quantum Chromodynam- 
ics (QCD), and is based on the gauge group SU{3)qcd [30-35]. This symmetry 



Table 2: The known fundamental interactions and their properties [29]. Gravitation 
is shown separately as it is not included in the SM of elementary particles. 



interaction 
type 


couples 
to 


affected 
particles 


exchange 
boson 


mass 
[GeV/c2] 


charge 
[e] 


spin 


strong 


color 
charge 


quarks, 
gluons 


gluon (g) 








1 


weak 


weak 
charge 


quarks, W^, 
leptons, Z° 




80.4 
91.2 


+ 1,-1 



1 
1 


electro- 
magnetic 


electric 
charge 


electrically 
charged 


photon (7) 








1 


gravitation 


mass, 
energy 


all 


graviton 
(unobserved) 
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is exact, and the gluons carry both a color and an anticolor charge. At increas- 
ingly short distances (or large relative momenta), the interaction gets arbitrarily 
weak (asymptotically free), thereby making a perturbative treatment viable. Via 
the strong interaction, quarks can form bound color-singlct states called hadrons, 
consisting of either a quark and an antiquark (mesons) or three quarks respectively 
antiquarks (baryons). The fact that only color- neutral states and no free quarks 
are observed is referred to as the confinement of quarks in hadrons. Due to its 
large mass, the top quark decays faster than the typical hadronization time of QCD 
{J^top ^ ^QCd), and it is therefore the only quark that does not form bound states. 
Its decay hence offers the unique possibility to study the properties of essentially a 
bare quark. 

The theory of clcctroweak interactions developed by Glashow [36], Salam [37] 
and Weinberg [38] is based on the SU{2)l x U{1)y gauge group of the weak left- 
handed isospin T and hypercharge Y. Since the weak {V — A) interaction only 
couples to left-handed particles, the fcrmion fields 4* arc split up into left-handed 
and right-handed fields '^l,r= ^(ITTs)* that are arranged in weak isospin T = ^ 
doublets and T = singlets: 

/c\ / A Ur Cr tR 

\dj ^ \sj^ \bj ^ dR Sr bR 

fl^t^\ / i^r \ l^eR I'lXR VrR 

V M / i V / i eR Ur TR 

In the doublets, neutrinos and the up-type quarks (m, c, t) have the weak isospin 
T3 = -|-|, while the charged leptons and down- type quarks {d,s,b) carry the weak 
isospin T3 = — i. The weak hypercharge Y is then defined via electric charge and 
weak isospin to be F = 2Q — 2T3. Consequently, members within a doublet carry 
the same hypercharge: y = — 1 for leptons and F = | for quarks, as implied by the 
product of the two symmetry groups. 

The SU{2)l x U{1)y gauge group does not accommodate mass terms for the 
gauge bosons or fermions without violating gauge invariance. A minimal way to 
incorporate these observed masses is to implement spontaneous electroweak sym- 
metry breaking (EWSB) at energies around the mass scale of the W and Z boson, 
often referred to as the "Higgs mechanism" [15-20], by introducing an SU(2) dou- 
blet of complex scalar fields $ = ($+, $°)'^. When the neutral component obtains 
a non-zero vacuum expectation value v/^/2 ^ 0, the SU{2)l x U{1)y symmetry 
is broken to U{1)QED^ giving mass to the three clcctroweak gauge bosons VK^, Z° 
while keeping the photon massless, and thereby leaving the electromagnetic symme- 
try U{1)qed unbroken. Prom the remaining degree of freedom of the scalar doublet, 
we obtain an additional scalar particle, the Higgs boson. 

The Higgs mechanism also provides fermion masses through fermion Yukawa 
couplings to the scalar field, with masses given by m/ = Xfv/V2, for a Yukawa cou- 
pling constant A/ for each massive fermion in the standard model. With a Yukawa 
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coupling close to unity, the top quark may play a special role in the process of mass 
generation. 

The mixing of flavor eigenstates in weak charged-current interactions of quarks 

is described by the Cabibbo-Kobayashi-Maskawa (CKM) matrix [39,40]. By conven- 
tion, this is a 3 X 3 unitary matrix Vckm that operates on the negatively-charged 
flavor states d, s and b: 

d'\ /VudVusVub\ /d\ /d\ 
s' = Vcd Vcs Vcb] \s \ = Vckm \s \ . (1) 
b'J ^ \Vtd Vs VtJ \bj ^ \bj ^ 

This complex matrix can have 18 independent parameters. However, to conserve 
the probability, this matrix has to be unitary, which means that there are only nine 
free parameters. An additional five out of the nine can be absorbed as phases in the 
quark wave functions. This results in four independent parameters in total - three 
real Euler angles and one complex phase, the latter implementing CP violation in 
the standard model. Since the CKM matrix is not diagonal, charged current weak 
interactions can have transitions between quark generations ( "generation mixing" ) 
with coupling strengths of the boson to the physical up and down type quarks 
given by the above matrix elements. 

From experimental evidence [41-47], neutrinos also have mass, which has led, 
among other things, to the introduction of an analogue leptonic mixing matrix, 
the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix [48,49]. It contains four 
independent parameters as well if one assumes that neutrinos are not Majorana 
particles. 

In summary, the standard model of elementary particle physics is a unitary, 

renormalizable theory [50,51], that can be used to perturbatively calculate processes 
at high energies. It incorporates 25 parameters that have to be provided through 
measurement: 

• 12 Yukawa couplings for the fermion masses 

• 8 parameters for the CKM and PMNS mixing matrices 

• 3 coupHng constants as,g,g' of SU{3)qcd, SU{2)l and U{1)y, respectively 

• 2 parameters from EWSB: v,mH- 

At currently accessible energy scales, the standard model describes successfully 
the interactions of fundamental fermions and gauge bosons, with only the Higgs 
boson remaining to be observed. For a more detailed introduction to the standard 

model, the reader is referred to corresponding textbooks, Refs. [52-54] on elementary 
particle physics and topical reviews such as Ref. [55]. 

2.2. The need for the top quark in the standard model 

The existence of the top quark was postulated well before its discovery mainly for 
three reasons. The first argument reflects the desire to have the standard model 
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Fig. 1: A "problematic" fermion triangle diagram that could introduce an anomaly. 

correspond to a renormalizablc theory. When expressed via a perturbation series - 
usually depicted in Feynman diagrams with first order "tree" diagrams and higher 
order "loop" terms - certain loop diagrams cause divergences that have to cancel 
exactly to ensure that the theory is renormalizablc. One example is the fermion 
triangle diagram, as shown in Fig.[TJ The contribution for each such diagram is pro- 
portional to c;^(5j, with c;^ being the weak neutral current axial coupling strength, 
and Qf the electric charge for the respective fermion in the loop. Since = T3 and 
neutrinos do not contribute, for the total strength of the anomaly to be cancelled, 
an equal number of lepton flavors and quark-doublets N families , and quarks in three 
colors {Nc = 3) are required [52]: 



Consequently, the discovery of the r lepton already called for an additional quark 
doublet to be present to keep the standard model renormalizablc. 

The second argument results from the fact that transitions that change the 
flavor but not the charge of a fermion (m <-!■ c f or 0? <-!■ s <-!■ 6) are observed to be 
strongly suppressed. The absence of such flavor-changing neutral currents (FCNC) 
for two quark generations could be accommodated through the GIM mechanism [56] 
by postulating the existence of the charm quark - and thereby completion of the 
second quark doublet - years before its discovery. This mechanism can be applied 
in a similar way for three quark generations, requiring a sixth quark as a partner 
of the b quark to complete the doublet. 

The third argument comes from the experimental confirmation that the b quark 
is not a weak isospin singlet but is part of an isospin doublet carrying the weak 
isospin = — i and electric charge Qb = —^e. The electric charge of the b quark 
was measured first at the electron-positron storage ring DORIS at DESY operating 
at the T and T' resonances through a measurement of the cross section for resonant 
hadron production ah [57-59]. The integral over at is related to the electronic 
partial width Fge, the hadronic partial width Th, the total width Ttot and the 
resonance mass Mr via / an dM = ^tt'^T h / {M^T tot) ■ Assuming that the total 
width is dominated by the hadronic partial width (F/j w Ttot), a measurement of 




) 



0. 



(2) 
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the integrated cross section and the resonance mass provides the electronic partial 
width Fee of the T and of the T'. In the framework of non-relativistic quarkonium 
potential models [60,61], this partial width can then be related with the bound 
quark's charge. 

The weak isospin of the b quark was measured via the forward-backward asym- 
metry ApB in the process e+e" bb ^ At^+ hadrons with the JADE detector at 
PETRA [62]. The asymmetry originates from electroweak interference effects and 
is defined as the difference between the number of fermions produced in the for- 
ward direction (with polar angle 9 < 90°) and the number of fermions produced 
backward {6 > 90°), divided by their sum. Aps is proportional to the ratio of 
the weak axial to the electric charge and vanishes for a weak isospin singlet. For a 
T3 = — i, Q = — ie 6-quark, the predicted asymmetry is —25.2%, in good agreement 
with the measurement of — 22.8± 6.0 (stat.) ± 2.5 (syst.)%. 

As a result of these measurements, the top quark's weak isospin and electric 
charge within the standard model were assigned to be T3 ~ +i, Qi = +|e, well 
before its discovery. The mass of the top quark, being a free parameter in the 
standard model, could not be predicted. Nevertheless, the mass of the top quark 
can indirectly be constrained by precision electroweak measurements. 

2.3. Top quark mass from precision electroweak measurements 

As discussed above, the standard model comprises a set of free parameters that arc 
a priori unknown. However, once these are measured, all physical observables can 
be expressed in terms of those parameters. To make optimal use of the predictive 
power of the theory, it is therefore crucial to measure its input parameters with 
highest possible precision, and thereby probe the self-consistency of the SM and 
any contributions beyond its scope. Being a renormalizable theory, predictions for 
any observable can be calculated to any order and checked with experiment. 

Electroweak processes depend mainly on three parameters: the coupling con- 
stants g and g' of SU{2)l and U{1)y, respectively, and the Higgs vacuum expecta- 
tion value V. Since these input parameters have to be obtained from experiment, it 
is best to substitute them with the most precisely measured quantities of the elec- 
tromagnetic fine structure constant a (using electron-positron annihilations into 
hadrons at low center of mass energies to measure hadronic vacuum polarization 
corrections [63-65]), the Fermi constant Gp (from the muon lifetime [66,67]) and 
the mass of the Z boson mz (from electron-positron annihilations around the Z 
pole [68]). 

With these input values, the theoretical framework can be used to predict other 
quantities such as the W boson mass. Given precision measurements, the W boson 
mass is sensitive to the mass of the top quark and the mass of the Higgs boson 
through higher order radiative quantum corrections [69-71]. 

The most precise electroweak measurements to date have been performed at the 
Large Electron-Positron (LEP) Collider [72, 73] at CERN by the four experiments 
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H 
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(b) 

Fig. 2: (a) Radiative contributions to the W and Z propagators from the top quark, 
(b) Radiative corrections to the W and Z propagators from the Higgs boson. 

ALEPH [74,75], DELPHI [76,77], L3 [78,79] and OPAL [80,81], and at the Stanford 
Linear CoUider (SLC) [82,83] by the SLD experiment [84,85]. The LEP experiments 
have analyzed ^17 miUion Z decays, and a sample of ~600 thousand Z bosons 
produced with longitudinally polarized electron beams was analyzed by SLD. 
Defining the electroweak mixing angle 0w via the vector boson masses: 

!!^ = l-sm'0w, (3) 

the W boson mass can be expressed as [70]: 

n 7ra 1 



where the radiative correction Ar is a directly observable quantum effect of elec- 
troweak theory that depends on a,mw,mz,mH and mj. The contributions from 
single-loop insertions containing the top quark and the Higgs boson, as depicted in 
Fig. El are [86]: 

3V2Gf cot^ Ow 



Ar*°P = ■ (for > nib) (5) 



167r2 V mfy 6 



(for ruH > mw)- (6) 



Thus, a precise measurement of W and Z boson masses provides access to the mass 
of the top quark and the Higgs boson. The top quark contribution to radiative 
corrections is large, primarily because of the large mass difference relative to its 
weak isospin partner, the b quark. While the leading top quark contribution to Ar 
is quadratic, it is only logarithmic in mass for the Higgs boson. Consequently, the 
constraints that can be derived on the mass of the top quark are much stronger 
than for the Higgs boson mass. 

In 1994, the most stringent constraints on mt were based on preliminary LEP and 
SLD data, combined with measurements of niw in proton- antiproton experiments. 
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LEP + SLD + Colliders + vq 




Fig. 3: (a) distributions of the standard model fit to precision electroweak data 
versus top quark mass for various Higgs boson masses [3]. (b) Comparison of the 
indirect top quark mass measurements via radiative corrections (shaded area) with 
the direct measurements from the Tevatron (points) versus time [68]. 



and neutral to charged-current ratios obtained from neutrino experiments, yielding 
178 ± 11 til GeV/c^ [3]. As illustrated in Fig. |3(a)[ the central value and the first 
set of uncertainties are from a fit of the SM to precision electroweak data, 
assuming uih = 300 GeV/c^. The second set of uncertainties stems from the impact 
of varying tuh between 60 and 1000 GeV/c^. 

The good agreement between predicted and observed values of rrit^ shown in 
Fig. |3(b)| as a function of time [68], is one of the great successes of the SM. The latest 
prediction from precision electroweak data yields rut = 179lg^ GeV/c^ without 
imposing constraints on mu [14], and is in excellent agreement with the current 
world average of mt = 172.4 ± 1.2 GeV/c^ [12]. 

This success of the SM also gives greater confidence in the predictions for mn- 
Since the precision of the prediction depends crucially on the accuracy of mw and 
rrit, it provides a strong motivation for improving the corresponding measurements. 
The current constraints on the mass of the Higgs boson will be discussed in Sec- 
tion [7331 

More details on precision electroweak measurements can be found in topical 
reviews such as given in Refs. [86,87]. 



3. Production and Decay of Top Quarks 

The production of top quarks is only possible at highest center of mass energies 
•y/i, set by the scale of mj. The energies needed for production of top quarks in 
the SM are currently (and will be at least for the next decade) only accessible at 
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hadron colliders. The Tevatron proton-antiproton collider started operation at y/s = 
1.8 TeV in 1987 for a first period of data taking ("Run 0") that lasted until 1989, 
with the CDF experiment recording about 4 pb^^ of integrated luminosity. The 
next data taking period from 1992 until 1996 at y/s = 1.8 TeV (the so-called Run I) 
was utilized by both the CDF and DO experiments and facilitated the discovery of 
the top quark. For the currently ongoing data taking that started in 2001 (Run II), 
the center of mass energy has increased to ^/s — 1.96 TeV. The Tevatron will lose its 
monopoly for top quark production only with the turning-on of the Large Hadron 
Collider (LHC) that will provide proton-proton collisions at -y/s = 14 TeV. 

In the framework of the standard model, top quarks can be produced in pairs {tt) 
predominantly via the strong interaction and singly via the electroweak interaction. 

3.1. Top quark pair production 

While hadron colliders provide the highest center of mass energies, the collision of 
hadrons complicates the theoretical description and prediction for processes such as 
tt production because of the composite nature of the colliding particles. These diffi- 
culties can be handled through the QCD factorization theorem [88,89] that provides 
a way to separate hadron collisions into universal long-distance (small momentum 
transfer) phenomena and perturbatively calculable short-distance phenomena. The 
latter processes involve therefore large square of momentum transfers Q^, and con- 
sequently the production of particles with large transverse momenta or large mass. 
The two components are set apart by introducing a factorization scale in the 
calculation. 

Using this approach, the proton can be described by a collection of partons 
(quarks, antiquarks, gluons) that interact at a low energy scale Aqcd < 1 GeV, 
whereas the elementary collisions between partons of the proton (or antipro- 
ton) occur on a "hard" energy scale characterized by large transverse momenta 
> 0(100 GeV). 

Consequently, the partons participating in any hard process (a, b) can be consid- 
ered quasi-free, and the partonic cross section of interest aa+b^x{s, q;s(/^^), /i^) can 
be calculated using perturbative QCD, independent of the type of hadrons contain- 
ing the partons. (The hatted variables denote parton quantities.) The regularization 
of divergences in higher order calculations (such as ultraviolet divergences from loop 
insertions, where the infinite range of four-momentum in the loop causes infinities 
in the integration from high momentum contributions) requires the introduction of 
a renormalization scale along with the corresponding running coupling constant 
as{^J.'jl). The leading order Feynman diagrams for tt production are shown in Fig.Hl 

The partons within the incoming proton (or antiproton) cannot be described by 
perturbative QCD, as the soft energy scale corresponding to small inherent momen- 
tum transfers implies large as(Q^) couplings. The distribution of the longitudinal 
momentum of the hadron among the partons is described by Parton Distribution 
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Fig. 4: Leading order Feynman diagrams for top quark pair production: (a) quark- 
antiquark annihilation and (b) gluon-gluon fusion. 

Functions (PDFs): /a//i(2;, Mf); corresponding to the probabiUty to find a given 
parton a inside hadron A with momentum fraction x when probed at an energy 
scale fj,^. Collinear and soft (infrared) singularities that arise in the perturbative 
calculation of the partonic cross section discussed above are absorbed in these PDFs. 

The factorization theorem is used to describe the tt production cross section via 
an integral over the corresponding hard scattering parton cross section, folded with 
the parton distribution functions of the incident hadrons as follows: 

aA+B^tt{Vs,mt) = ^ I aa+b^ti{s,as{HR),I^R,I^F,mt) 

a,b=g,q,q 

X fa/A{Xa,IJ'F)fb/B{Xb, IJ.%)dXadXb. (7) 

The hadrons A and B correspond to proton and antiproton in case of the Tevatron 
and to protons in case of the LHC. 

The physical cross section o■J^_^_g^ll{y^, mt) that would emerge from the evalu- 
ation of the full perturbation series does not depend on either of the two arbitrary 
scales for factorization and renormalization that had to be introduced for 

the calculation. However, the parton distribution functions and the partonic cross 
section do depend on these scales, and hence the result of any finite order calcu- 
lation will as well. This dependence gets weaker with the inclusion of higher order 
terms in the calculation. In practical application, both scales are usually set to the 
typical momentum scale of the hard scattering process, such as the transverse 
momenta of the produced particles or the mass of the produced particle, so that 
for ti production, typically ij,f = Hr = = mt- The scale dependence of the result 
is then usually tested by varying the central scale by a factor of two; the resulting 
variations are interpreted as systematic uncertainties that should not be mistaken 
as Gaussian in nature. 

The PDFs have to be determined experimentally, for example via deeply inelas- 
tic lepton scattering on nucleons, so that they can be extracted from the measured 
cross sections using perturbative calculations of the (hard) partonic cross sections. 
Once the parton densities fa/AixajQ^) have been measured as a function of mo- 
mentum fraction Xa at a scale Q'^, their value at a different scale can be predicted 
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Fig. 5: (a) CTEQ61 [93] parton distribution functions with their uncertainty bands, 
for the rrit mass scale (Q^ — (170 GeV)-^) for (anti-) up quarks, (anti-) down quarks 
and gluons in the proton, (b) Relative uncertainties on the PDFs shown in (a). 

perturbatively using the DGLAP evolution equation [90-92]. Since PDFs are uni- 
versal and do not depend on the process they were derived from, they can be used to 
predict cross sections in other hard scattering processes. For consistent application, 
it is important that the PDFs are derived to same perturbative order and with the 
same renormalization scheme as the calculation of any prediction. 

The PDFs are extracted from global fits to the available data, as is done, for ex- 
ample, by the CTEQ [94], MRST [95], GRV [96], Alekhin [97], HI [98] and ZEUS [99] 
groups. Different PDFs are based on different data, different orders of perturba- 
tion theory, renormalization schemes and fitting techniques - see, for example, the 
overview given in Ref. [100]. One commonly used set of PDFs derived at NLO, us- 
ing the MS renormalization scheme [101], is CTEQ61 [93], which incorporates the 
Tevatron Run I data on jet production, especially important for the gluon distribu- 
tion. CTEQ61 also includes an error analysis based on different sets of PDFs that 
describe the behavior of a global function for the fit around its minimum. The 
resulting error on the PDF (A/) can be obtained by summing over the variations 

along/against each PDF "eigenvector" for every free parameter in the global fit: 

A/ = ±i(i];ir(/+-/r)2)i 

Figure[5]shows the most important parton distributions within protons for tt pro- 
duction at the Tevatron or LHC, and their corresponding uncertainties. (For an- 
tiprotons, quarks and antiquarks have to be interchanged in Fig.[5l) All PDFs vanish 
at large momentum fractions x, and the gluon density starts to dominate over the 
valence-quark densities near x = 0.13. There is no flavor symmetry between the 
up and down quark distributions, neither on the valence nor the sea quark level 
(the latter is best seen at low Q^). At a;- values below 0.1, typical relative uncertain- 
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ties on the PDFs of valence quarks and gluons are ^5%. At larger a;-values, these 
uncertainties increase drastically, especially for gluons. 

To produce a top quark pair, the squared center of mass energy at the par- 
ton level s = XaXbS must at least equal (2mt)^. Assuming Xa ~ Xb = x yields as 
threshold for tt production: 



Since large momentum fractions are required for tt production at the Tevatron, the 
process is dominated by quark- antiquark annihilation (Fig. [4] (a)) of the valence 
quarks. For Run I energies, quark- antiquark annihilation contributes roughly 90% 
of the total ti production rate, and for Run II energies this fraction is w85% [25]. 

At the LHC, gluon-gluon fusion dominates (Fig. |4] (b)) with a contribution 
of about 90% [102], because a small momentum fraction suffices for tt production. 
This means that proton-proton collisions at the LHC have production cross sections 
comparable to pp rates, thereby obviating the need for the major technical challenge 
of producing an intense antiproton beam. 

The increase in the center of mass energy by «10% between Run I and Run II at 
the Tevatron and the correspondingly smaller minimum momentum fraction provide 
an increase in the ti production rate of 30%. At the LHC, the rate increases by 
roughly a factor of 100 compared to that of the Tevatron. 

The highest-order complete perturbative calculations for heavy quark pair pro- 
duction have been available at next-to-leading order (NLO) - to order - since 
the late 1980s from Nason et al. [103] and Beenakker et al. [104,105]. These cal- 
culations can be refined by the inclusion of large logarithmic corrections [106-108] 
from soft-gluon emission that are particularly important for the production of heavy 
quarks close to the kinematic threshold (s k, 4m^, a; — > 1). The contributions of 
these logarithms are positive at all orders when evaluated at the heavy quark mass 
scale and their inclusion therefore increases the production cross section above the 
NLO level. 

The impact of soft-gluon resummation on the tt production cross section has 
been studied by Berger and Contopanagos [109-111], Laenen, Smith and van Neer- 
ven [112,113] and Catani, Mangano, Nason and Trentadue [114,115] at the leading 
logarithmic (LL) level. Studies including even higher level corrections as carried 
out by Cacciari et al. [116], based on work by Bonciani et al. (BCMN) [117], and 
Kidonakis and Vogt [118,119] are summarized in Tabled 

In the case of ti production at the Tevatron, the inclusion of leading and next- 
to-leading logarithmic (NLL) soft-gluon resummation affects the cross sections only 
mildly by C(5%) (indicating production occurs not too close to the threshold), 
while significantly reducing the scale dependence of the predictions by roughly a 
factor of two to a level of «5% [117]. At the LHC, it production takes place even 
further away from the kinematic threshold, but since gluon fusion dominates there. 




2mt 



0.192 @ Tevatron Run I, Vs = 1.8 TeV 
0.176 @ Tevatron Run II, ^ = 1.96 TeV 
0.025 @ LHC, = 14 TeV 



(8) 
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Tabic 3: NLO cross section predictions including soft-gluon rcsummations beyond 
LL for tt production at the Tevatron and the LHC for a top quark mass of 175 
GeV/c^. For the different sources of the quoted uncertainties please refer to the 
text. 



Hadron Collider 


Processes 




^tt [pb] 
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90% qq 


tt 


c -|q+0-52 


Cacciari et al. [116] 


{pp. ^~ = 1.8 T(>\') 


10'/ gg — 


It 


.").24± ().:-!l 


Kidouakis (-/ al. [119] 


Tevatron Run II 


85% qq 


tt 




Cacciari cl al. [116] 


{pp. ^/T: = 1.90 TeV) 


to/( (jg 


IT 


().77± 0.42 


Kidouakis cj al. [119] 


LHC 


10% qq 


tt 


QQQ + 52 

8dd_39 


Bonciani et al. [117] 


{pp, ^=IA TeV) 


90% gg ^ 


tt 


873^^8 


Kidonakis et al. [118] 



the enhancement of the total production rate from soft-gluon resummation and the 
reduction of scale dependence stay at the same level as at the Tevatron. 

The results of Cacciari et al. [116] for the Tevatron use the NLO calculation 
with LL and NLL resummation at all orders of perturbation theory as carried out 
by Bonciani et al. (BCMN) [117], but are based on the more recent PDF sets with 
error analysis CTEQ6 [94] and MRST2GG1E [95] and also MRST2G01 [120], which 
includes varied ag values in the PDF fit. The updated PDFs cause an increase in the 
central values of about 3% relative to Ref. [117]. While the central values are very 
similar for the MRST2001E and CTEQ6 PDFs, the uncertainties for CTEQ6 are 
almost twice as large as for MRST2001E, unless the variations in ag in MRST2001 
are also included. For the determination of the uncertainty on the cross section, 
Cacciari et al. combine linearly the uncertainty due to scale variation by a factor 
of two with the PDF uncertainty evaluated at that scale. As central values, the 
CTEQ6M results are chosen, and the maximum lower (upper) uncertainties given 
stem from the CTEQ6 PDF variation (the variation in MRST2001). The PDF 
uncertainties and as variation contribute about 45% and 80% respectively to the 
total quoted uncertainty, including the scale variations, which emphasizes the im- 
portance of considering ag uncertainties in PDF fits. The PDF uncertainties are 
in turn dominated by the uncertainty of the gluon PDF at large x values, causing, 
for example, the gluon fusion contribution to the total production rate to fluctu- 
ate between 11% and 21% for — 1.96 TeV. Despite the large uncertainties on 
the tt production rate, the ratio of production cross sections for the two center 
of mass energies at the Tevatron is very stable and predicted with high precision: 
cr(1.96 TeV)/cr(1.8 TeV) = 1.295 ± 0.015 for top quark masses between 170 and 
180 GeV/c2 [116]. 

A prediction for the tt production rate at the LHC applying the same level 
of soft-gluon resummation is given by Bonciani et al. [117] using the MRS(R2) 
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PDF [121]. Since no PDF uncertainties were available for Ref. [117], the quoted 
uncertainty in Table [3] comes from changing the scale by factors of two alone. Since 
gluon fusion is the dominant contribution to the total rate, uncertainties on the 
gluon PDFs alone lead to an uncertainty of w 10% on the total production cross 
section [122]. 

The studies performed by Kidonakis and Vogt [118, 119] consider soft-gluon 
corrections up to next-to-next-to-next-to leading logarithmic (NNNLL) terms at 
NNLO in a truncated resummation, resulting in a reduced sensitivity of «3% to 
scale variations. For the Tevatron, the it production cross section is evaluated using 
MRST2002 NNLO [95] and CTEQ6M NLO [94] parton densities. Two different 
parton-level conditions are considered for the scattering process: (i) one-particle 
inclusive (IPI) and (ii) pair-invariant mass (PIM) kinematics [102]. While both sets 
of PDFs give very similar results, the variations from the difference in kinematics 
are significant. Consequently, the average of IPI and PIM kinematics for both PDFs 
is used as the central value in Table [3l while the separate averages over the PDFs 
for IPI and PIM are quoted as uncertainties. For the predicted LHC rate, which 
is dominated by gluon fusion, the IPI kinematics is considered more appropriate, 
and the value in Table [3] gives the corresponding result based on MRST2002 NNLO 
PDFs, using scale changes by factors of two for estimating the uncertainties. 

All results in Table [3] are evaluated for a top quark mass of 175 GeV/c^, and 
the Run II values serve as the main predictions for CDF and DO. To improve 
comparability of the uncertainties on the different predictions, the calculation by 
Kidonakis and Vogt has an additional uncertainty obtained from the maximum 
simultaneous changes in scale and PDFfE added in quadrature with the uncertainty 
due to the dependence on kinematics [123]. 

In spring 2008, Cacciari et al. [124] and Kidonakis et al. [125] updated their 
predictions using more recent PDFs such as CTEQ6.6M [126], which had only lit- 
tle impact on the results. In addition, Moch and Uwer [127] have now performed 
a complete NNLL soft-gluon resummation and provide an approximation of the 
NNLO cross section also based on CTEQ6.6M. To illustrate the dependence of the 
predictions on the top quark mass, Fig.[B]shows the central values and uncertainties 
from References [124,125,127] for Tevatron Run II versus mt- An exponential form, 
as suggested in Ref. [114], is applied in a fit to the central values and uncertain- 
ties for Kidonakis et al., while third-order polynomials, as provided by the authors, 
are used for the other references. The total uncertainties are obtained by linearly 
combining the provided uncertainties. 

For the current world-averaged top quark mass of 172.4 ± 1.2 GeV/c^ [12], the 
predicted tt production cross section is 7.2tQ g pb for Cacciari et al, l.'it.^^ pb 
for Kidonakis et al., and 7.5lQy pb for Moch et al.. An additional uncertainty of 
±0.3 pb arises from the uncertainty on the top quark mass for all three predictions. 
It should be noted that these predictions based on MRST 2006 NNLO PDFs [128] 

''The PDF uncertainties in this case stem from CTEQ6 sets "129" and "130" alone. 
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Fig. 6: Dependence of the tt production cross section on rrit in Run II of the Teva- 
tron. The current world-averaged top quark mass and the resulting expected tt pro- 
duction cross section are indicated by the vertical and horizontal lines, respectively. 
The predictions are based on the CTEQ6.6M PDFs [126]. 



yield about 6% higher central values and exhibit smaller uncertainties from PDFs. 

A precise measurement of the tt production cross section provides a test of the 
predictions for physics beyond the SM. Together with a precise mass measurement, 
the self-consistency of the predictions can also be examined. Because tt production 
is a major source of background for single top production (to be discussed in the 
next section) , standard model Higgs boson production and many other phenomena 
beyond the SM, its accurate understanding is crucial for such studies. 

Figure [7] illustrates the production rates of various processes versus center of 
mass energy for proton-antiproton collisions below ^/s = 4 TeV and for proton- 
proton collisions above ^/s — 4 TeV. As can be appreciated from the plot, tt pro- 
duction is suppressed by ten orders of magnitude relative to the total interaction 
rate at the Tevatron and eight orders of magnitude at the LHC. While the LHC is 
often referred to as a "top-factory" because of the increased production cross section 
by two orders of magnitude, extraction of the signal from the large background is 
a challenge at both hadron colliders, requiring efhcicnt triggers and selection meth- 
ods. The ti cross section measurements performed in Run II of the Tevatron will be 
described in Section [Ql 



3.2. Single top quark production 

In addition to the strong pair production discussed in the previous section, top 
quarks can also be produced singly via the electroweak interaction through a Wtb 
vertex (see Fig. [8]). Wts and Wtd vertices are strongly CKM suppressed (see Sec- 
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Fig. 7: Cross sections for various processes at hadron colliders as a function of center 
of mass energy [122]. (Jtoppairs denotes the ti production cross section. 



tion l3.3.T)) . There are three different production modes, classified via the virtuality 
(negative of the square of the four-momentum q) of the participating W boson 

• The Drell-Yan-like s- channel production proceeds via quark- antiquark annihila- 
tion into a time-like virtual W boson {q^ > [rat + nib)^ > 0), as illustrated in 
Fig.[8(a)| qq' tb [129,130]. 

• In the t-channel "flavor excitation" process, a space-like virtual W boson {q^ < 0) 
couples to a 6 quark from the nucleon's sea to produce a top quark, as shown in 

q't. A higher order contribution of 0{as) comes from gluon 
which is also referred to as W-gluon fusion for 



Fig. 8(b) for qb 



splitting, as depicted in Fig 
qq tq'b [131-134]. 

In associated production, an on-shell W boson (q^ 



) is produced together 



(b) 



{ 



(d) 
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Fig. 8: Representative Feynman diagrams for electroweak single top quark produc- 
tion: (a) s-channel, (b,c) t-channel and (d,e) associated production. 



with a top quark from a b quark and a giuon, as illustrated in Figs. [8(d)] and 8(e) 
tor gb-^tW [135-140]. 

In the above discussion, charge conjugate processes are implied for each pro- 
duction mode, and q represents a light-flavor quark. All three modes differ in both 
their initial and final states, and the processes are simply denoted as s-channel {tb), 
t-channel (tq, tqb) and associated (tW) production. The corresponding signatures 
can be used to discriminate between the production modes: The s-channel is char- 
acterized by an additional b quark accompanying the top quark, the t-channel by a 
forward light quark, and associated production by the decay products of the W bo- 
son in addition to those of the top quark. Due to the incoming b quark and gluon, 
the t-channel and tVF-channel rates are especially sensitive to the corresponding 
PDFs, which are known with less precision than the PDFs for the valence quarks of 
the proton. The measured cross sections will therefore provide further constraints 
on the b quark and gluon PDFs. 

The cross sections for all three modes have been evaluated at NLO, including 
radiative corrections of 0{as): s-channel [141-143], t-channel [143-145], and tW- 
channel [146,147]. (The most recent references provide differential distributions.) 
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Table 4: Cross sections for s-channel {tb) , t-channel {tq) and associated {tW) single- 
top production at NLO (Sullivan), and NLO, including soft-gluon resummations 
(Kidonakis), expected at the Tevatron and at the LHC, for a top quark mass of 175 
GeV/c^. (For sources of the quoted uncertainties see the text.) 
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Subsequent calculations also include top quark decay at NLO for s-channel [148- 
150], i-channel [148,150,151], and fVF-channel [152], and latest NLO calculations 
include higher-order soft-gluon corrections up to NNNLO at NLL accuracy [153- 
155]. 

Table U summarizes the expected single-top production cross sections at the 
Tevatron and the LHC for the NLO calculations by Sullivan [156] (based on the 
work of Harris et al. in [143]) and NLO results including soft-gluon resummations 
by Kidonakis [153] (based on his work in Refs. [154,155] and matching to the exact 
NLO results of Harris et al. [143] and Zhu [146,147]). Both results use current PDFs 
and include corresponding uncertainties. 

While top and antitop production are identical at the Tevatron for all production 
modes, at the LHC this is only the case for associated production. Consequently, 
the results given for the Tevatron include both top and antitop production but are 
given separately for the LHC. 

The NLO results of Sullivan are based on CTEQ5M1 PDFs [157] for their central 
values. The uncertainties for PDFs are derived from CTEQ6M [94], and added in 
quadrature with uncertainties from changes in scale by the usual factors of two, 
changes in top quark mass by 4.3 GeV/c^ (using an older world-averaged rrit = 
178 ± 4.3 GeV/c^ [158]), and uncertainties in b quark mass and as, the latter two 
being negligible. The rate dependence on the top quark mass is approximated as 
linear and is especially important for the s-channel, since a change from 175 GeV/c^ 
to the current world-averaged mt — 172.4 GeV/c^ raises the rates at the Tevatron 
by 7% for the s-channel and 5% for the t-channel. The observed uncertainties in 
scale are reduced relative to LO results, and amount to 4-6% at the Tevatron and 
2-3% at the LHC. 
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The NLO calculations of Kidonakis that include higher order soft-gluon cor- 
rections provide single-top production cross sections based on MRST2004 NNLO 
PDFs [159]. The quoted values are obtained by matching the NLO cross section to 
the results of Harris et al. [143] and Zhu [146,147], and including the additional 
soft-gluon corrections up to NNNLO. Exceptions are the tW rate at the Tevatron, 
where no corresponding NLO result is available, and the given value is therefore not 
matched, and the f-channel rate at the LHC, where no soft-gluon corrections are 
considered and an updated NLO result with the quoted PDFs is given instead. The 
uncertainties given are derived from varying the scale by a factor of two, and adding 
in quadrature PDF uncertainties derived using the MRST2001E NLO PDFs [95]. 
No uncertainty in mt is included. At the Tevatron, the t-channel uncertainty is dom- 
inated by the uncertainty in PDFs, and corrections from soft-gluon resummations 
relative to LO are small (~5%). In contrast, the soft-gluon corrections have a large 
effect (>60%) for the s-channel at the Tevatron and scale uncertainties dominate 
over those from PDFs. 

At the Tevatron, t-channel production dominates the total rate of single top 
quark production with a contribution of ~65%, followed by s-channel production 
at Ri30%. Associated {tW) production at the Tevatron contributes only ss5% to the 
total rate, and is usually neglected. At the LHC, t-channel production again dom- 
inates at ~74%, followed now by associated production at ^23%, while s-channel 
production contributes only Ri3% because of the missing contribution from valence 
antiquarks in the collisions, which will make it difficult to discriminate this chan- 
nel from background. Despite being an elcctroweak process, single top production 
has a cross section of the same order of magnitude as ti production (of 0(40%) of 
the it rate at both the Tevatron and the LHC). With only one heavy top quark 
to be produced, single top production is accessible at smaller and therefore better- 
populated momentum fractions of the partons. Furthermore, no color matching is 
required for the production. The fact that the observed yields of single-top and 
tt-pairs are consistent with theory is a major triumph of the SM. 

The measurement of single top production offers a check of the top quark's weak 
interaction, and direct access to the CKM matrix element \Vtb\, as the cross sections 
in all three production modes are proportional to |Vt(,p. The polarization of the top 
quark at production is preserved due to its short lifetime and provides a test of 
the V — A structure of the weak interaction via angular correlations among the 
decay products [160-163]. All three production modes provide different sensitivity 
to various aspects of physics beyond the standard model (BSM) [164], which makes 
their independent reconstruction a desirable goal. The s-channel is sensitive to the 
existence of new charged bosons (such as W or charged Higgs) that couple to the 
top-bottom weak-isospin doublet, an effect that could be detectable through an 
enhancement of the observed cross section. Such effects would not be observed in 
the tW mode, where the W boson is on-shell, or in the t-channel, where the virtual 
W boson is space-like and cannot go on-shell as in the s-channel. The f-channel 
production rate could be enhanced via FCNC processes involving new couplings 
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between the up-type quarks and a boson (Higgs, gluon, photon, Z). This would be 
hard to observe in the s-channel, since there is no b quark in the final state, which is 
essential for discrimination of the signal in that production mode. Finally, the tW 
channel is the only mode that provides a more direct test of the Wtb vertex since 
the W boson appears in the final state. 

A thorough understanding of single top quark production will also facilitate the 
study of processes exhibiting a similar signature such as SM M^-Higgs production 
or BSM signals to which single top production is a background process. Despite a 
production rate similar to that of it, the signature for single top quark production 
is much harder to separate from background, which has delayed first measurements 
until very recently. The current analyses at the Tevatron provide first evidence for 
production of single top quarks, and this will be described in Section [5.81 



3.3. Top quark decay 

3.3.1. Top quark C KM matrix elements 

Since the mass of the top quark is larger than that of the W boson, decays t — > Wq, 
with q being one of the down-type quarks d, s, b, are dominant. The contribution 
of each quark fiavor to the total decay width is proportional to the square of the 
respective CKM matrix element Vtq. Utilizing the unitarity of the CKM matrix 
and assuming three quark generations, the corresponding matrix elements can be 
constrained indirectly at 90% confidence level to [165]: 

\Vtd\ = 0.0048 - 0.014 (9) 
\Vts\ = 0.037 - 0.043 (10) 
\Vtb\ = 0.9990 - 0.9992. (11) 

Consequently, the decay t Wb is absolutely dominant and will be considered 
exclusively throughout this article, unless noted to the contrary. Potential deviations 
from the SM decay will be discussed in Section [6l 

It should be noted that the above constraints on the CKM matrix elements would 
change dramatically (especially Vtb) if there were more than three quark generations. 
Assuming the unitarity of the expanded matrix, the limits become [165]: 

\Vtd\ = - 0.08 (12) 
\Vts\ = - 0.11 (13) 
\Vtb\ = 0.07 - 0.9993. (14) 

It is therefore important to constrain these matrix elements through direct mea- 
surements, as outlined below. 

The Vtd and Vts matrix elements cannot be extracted from lowest-order (tree 
level) top decays in the framework of the standard model, but can be inferred from 
B-meson mixing, as shown in Fig. [9l While all up-type quarks can contribute in the 
depicted box diagrams, the contribution from the top quark is dominant [166]. The 
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Fig. 9: Feynman "box" diagrams for B'-j — and i?° — B'^ mixing. 

oscillation frequency given by the mass difference between heavy and light mass 
eigenstates, Amd for B^ — B^ and Atos for _B° — B^ oscillations, is proportional 
to the combination of CKM matrix elements iV^jVidp and iVjJVtsP, respectively. 
The mass difference for the B° - system is Am^ = 0.507 ± 0.004 ps'^ [167]. 
Using CKM unitarity and assuming three generations, yielding \ Vtb\ ~ 1, translates 
into Vtd — (7.4 ± 0.8) • 10~^ [168], where the uncertainty arises primarily from the 
theoretical uncertainty on the hadronic matrix element, which is obtained from 
lattice QCD calculations. To reduce these theoretical uncertainties, a measurement 
of the ratio, in which some uncertainties cancel, is more desirable (Am^i/Ams oc 
I VtdP/l Vtsp). With the recent first measurement of Am^ in i?°-oscillations by DO 
and CDF at the Tevatron [169,170], yielding 17 ps"^ < Am^ < 21 ps'^ at 90% C.L. 
and Atos = (17.3llQ '^g(stat.) ± 0.07(syst.)) ps^^, this ratio has now been measured 
for the first time as iVtd/VtsI = 0.208tQoo2(expt.)^QQQg(theor.). These results are 
in good agreement with SM expectations. 

The direct measurement of the Vth matrix element without assuming three quark 
generations and unitarity of the CKM matrix is only possible via single top quark 
production (described in Section [3. 2|) . because the production rate in each channel 
is proportional to |Vtbp. One way to assess the relative size of \Vth\ compared to 
\Vtd\ and \Vts\ is to measure the ratio R of the top quark branching fractions, which 
can be expressed via CKM matrix elements as 



B(^ ^ Wb) ^ \Vtb^ , . 

B{t ^ Wq) I Vtb ? + I Ft. P + I Vu V ^ ^ 

Assuming three generation unitarity, the denominator in the above expression equals 
one, and constraints on |Vtfc| can be inferred. The current status of these measure- 
ments is discussed in Section [6.21 

The most precise extraction of the top quark CKM matrix elements proceeds 
via global fits to all available measurements, imposing the SM constraints of three 
generation unitarity, as done by the CKMfitter [171] or UTfit [172] groups. The 
CKMfitter update for summer 2008 yields [171]: 

\V,d\ = 0.00853tO;S (16) 
\V,,\ = 0.999146tO:S^. (18) 
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3.3.2. Decay width of the top quark 

The decay width of the top quark in the SM, including first-order QCD corrections, 
can be expressed as follows [173,174]: 

where the above formula assumes ml/rn^ — > 0, ^ and ignores corrections 

2 

of O ( "° ™^ ) and 0{a^). While the above QCD corrections lower the width by 
«10%, first-order electroweak corrections increase the width by 1.7% [175,176]. 
However, the electroweak correction is almost cancelled when the finite width of the 
W boson is taken into account, thereby decreasing the width again by 1.5% [177, 
178]. Corrections to the top quark width of 0{al) have also been evaluated [179,180] 
and reduce the width by 2%. Including all these effects, the decay width is predicted 
to a precision of ^1%. The other SM decays, t Wd and t Ws, contribute 
negligibly to the total decay width Tt = S^Ftg because of proportionality to |Vtd|^ 
and \Vts\^. 

Equation [12] yields the top width to better than 2% accuracy, and the width 
increases with rrit. For as{mz) = 0.1176 and Gp = 1.16637- 10~5(GeV/c^)"^ [47], 
Ft is 1.02/1.26/1.54 GeV for top quark masses of 160/170/180 GeV/c^. 

The resulting lifetime of the top quark tj — hT^^ « ?i(1.3 GeV)^^ is approx- 
imately 5 • 10~^^ s, and significantly shorter than the hadronization time Thad ~ 
^^QCD ~ ^(0.2 GeV)^^ « 3 • 10"^'' s. As a consequence, the top quark decays 
before it can form hadrons, and in particular there can be no tt bound states (topo- 
nium), as was already pointed out in the 1980s [13, 181, 182]. Nevertheless, although 
the top quark can generally be considered as a free quark, residual non-perturbative 
effects associated with hadronization should still be present in top quark events, and 
the fragmentation and hadronization processes will be influenced by the color struc- 
ture of the hard interaction. 

In electron-positron annihilation, top quark pairs are produced in color singlet 
states, so that hadronization before decay depends mainly on the mass of the top 
quark and collision energy. In hadronic tt production, t and t are usually produced in 
color octet states and form color singlets with the proton and antiproton remnants. 
The energy in the color field (or in the string when using the picture of string 
fragmentation) is proportional to the distance between top quark and the remnant. 
If a characteristic length of about 1 fm is reached before the top quark decays, 
light hadrons can materialize out of the string's energy. The possibility for such 
string fragmentation depends strongly on the center of mass energy in the hadron 
collisions. For Tevatron energies, this can be neglected [183], while it may be more 
important at LHC energies, where top quarks are produced with sizeable Lorentz 
boosts. Since heavy quarks have hard fragmentation functions and the fractional 
energy loss of the top quarks is therefore expected to be small, it will be difhcult to 
experimentally establish these effects directly, even at the LHC. In case no string 



, o Gf tni ( Tow \ / TOT,, \ 
Vt^\Vthf ^^[\-^\ 1 + 2^ 
87rV2 V \ / 



Top Quark Measurements at the Tevatron 25 



? 



'W 



'W 



longitudinal (Wq) left-handed right-handed {W+) 

(forbidden if — 0) 

Fig. 10: Angular momentum conservation in top quark decay does not allow right- 
handed W bosons when h quarks are assumed massless. 



fragmentation takes place before the top quark decays, long-distance QCD effects 
will still connect the decay products of the top quark. 

With top quark mass measurements aiming at uncertainties of < 1 GeV/c^, it 
becomes more and more important to assess the impact of such non-perturbative 
effects on the measurements. One example that may play an important role in 
this context is the possibility of color reconnections before hadronization, and the 
corresponding modeling of the underlying event (beam-remnant interactions) [184, 
185]. 



3.3.3. Helicity of the W boson 

Top quark decay in the framework of the standard model proceeds via the left- 
handed charged current weak interaction, exhibiting a vector minus axial vector 
{V — A) structure. This is reflected in the observed helicity states of the W boson, 
which can be exploited to examine the couplings at the Wtb vertex [186-188]. 

The emitted b quark can be regarded as massless compared to the top quark, and 
hence expected to be predominantly of negative helicity (left-handed), meaning that 
its spin points opposite to its line of flight. The emitted W boson, being a massive 
spin-1 particle, can assume any of three helicities: one longitudinal (Wq) and two 
transverse states (W^, left-handed and W+, right-handed). To conserve angular 
momentum in the t Wb decay, the spin projection of the W boson onto its 
momentum must vanish if the b quark's spin points along the spin of the top quark, 
while a left-handed W boson is needed if the b quark's spin points opposite to the 
spin of the top quark. In the limit of a massless b quark, a right-handed W boson 
cannot contribute to the decay, as illustrated in Fig.[TOl For the decay of an antitop 
quark, a left-handed W boson is forbidden. 

At lowest "Born" -level, the expected fractions of decays with different W boson 
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helicities, taking the finite b quark mass into account, are given by [189]: 

* (1 -y2)2+ 2.2(1 _ 2a;2 + y2) ~ i_^ 2x2 ^ ' 



f =r /r - + ^ 2x2 

* (1 -?;2)2+ 2.2(1 _ 2x2 + y2) ~ J' (1_ 2-2)2(1 + 2x2) ^ ^ 

where the scaled masses x = mw/rrit, y — mh/mt and the "Kallen"-type function 
A = 1 + x'' + 2/'* — 2x2^/2 — 2x2 — 2y'^ were used. Inserting mt — 175 GeV/c2, 
mw — 80.419 GeV/c2, and a pole mass of rrib — 4.8 GeV/c2, the partial helicity 
rates are found to be [189]: 

/o = 0.703, /_ = 0.297, /+ = 0.00036. (23) 

TOb ^ results in a reduction of /o and /_ at the per mill level. Including one-loop 
QCD corrections [189], electroweak one- loop corrections and finite width corrections 
[190] leads basically to a cancellation of the last two corrections, as was already 
mentioned in Section [3.3.2[ and a change in partial helicity rates /o and /_ at the 
1-2% level. The right-handed helicity fraction /+ remains at the per mill level with 
these corrections included. Consequently, any observation of /+ at the percent level 
would signal the presence of physics beyond the standard model. 

Using a more general extension to the standard model Wtb interaction La- 
grangian, assuming that both the W boson and b quark are on-shell, leads to [186]: 



^/2 



-/i.e. (24) 



where Pr{l) are the right- and left-handed chiral projectors Pr{l) = ^(1 i 7^) 
and ia'^'^ = — ^[7^,7*^]. This model-independent extension has four form factors 

/i/2^' ^^'^ includes the standard model as a special case, with /f = 1 (left-handed 
vector coupling) and the other form factors (right-handed vector, and left- and right- 
handed tensor couplings) vanishing. These four general couplings for the Wtb vertex 
can be determined by measuring four observables sensitive to this interaction: the W 
helicity fractions /o and /+ in ti events, and the single-top production cross sections 
in the s- and f-channel. This model-independent determination of the parameters of 
the Wtb vertex can in turn be used to distinguish between different models proposed 
for EWSB [191]. 



4. Experimental Setup 

This section describes the experimental ingredients that are needed to study top 
quarks. Since this review focusses mainly on results obtained at Run II of the 
Tevatron, only the corresponding accelerator and detector setups are discussed. 
The experimental setup for Run I can be found, for example, in Ref. [192]. 



Top Quark Measurements at the Tevatron 27 




TEVATRON ' 




p (1 TeV) 
- 

p (1 TeV) 
D0 DETECTOR 



p ABORT 



Fig. 11: The Fcrmilab Run II accelerator complex consisting of the Tevatron pp 
collider and its pre-accelerators [193,196]. 



The Tevatron collider is discussed in the first part of this chapter, followed by a 
description of the two general-purpose detectors CDF and DO surrounding the two 
interaction points where protons and antiprotons are brought to collision. Subse- 
quently, the reconstruction and identification of particles produced in the collisions 
are briefly reviewed before the resulting experimental signatures of top quark events 
are described. Finally, the Monte Carlo (MC) simulation tools needed to model in- 
teractions in the detectors are considered. 



4.1. The Tevatron collider 

The Tevatron collider is part of the Fermi National Accelerator Laboratory (Fermi- 
lab, FNAL) in Batavia, Illinois, close to Chicago. Until the Large Hadron Collider 
at CERN starts operation, the Tevatron remains the particle accelerator with the 
highest center of mass energy worldwide. Here, 36 bunches of protons and antipro- 
tons with a spacing of 396 ns are brought to collision at -^s = 1.96 TeV at the two 
interaction points where the multi-purpose detectors CDF and DO reside. As illus- 
trated in Fig. [HI the Tevatron is the final stage in a chain of eight pre-accelerators 
and storage rings [193-196]. 

Beam protons are generated using a magnetron surface plasma source that pro- 
duces ions from hydrogen gas [197]. The ions are then accelerated to 750 keV 
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in a Cockcroft- Walton electrostatic accelerator, followed by a linear accelerator, 
bringing the ions to 400 MeV. Using a carbon stripping foil, both electrons are re- 
moved from the H~ ions, and the resulting protons are then accelerated to 8 GeV 
within 33 ms in the first of five synchrotrons, called the "Booster" , which has a ring 
circumference of 475 m. All these components are often referred to as the "Proton 
Source" . 

Acceleration continues in the oval Main Injector synchrotron that has a circum- 
ference of 3.3 km. Depending on their ultimate use. protons are either brought from 
8 GeV to 120 GeV within 2 s for fixed-target operation (including the production of 
antiprotons) or to 150 GeV within 3 s for injection into the Tevatron. With a ring 
radius of 1 km, the Tevatron is the final and largest synchrotron at Fermilab. It 
accelerates protons and antiprotons in a single beam pipe from 150 GeV to 980 GeV 
in about 85 s. 

For the production of antiprotons, 120 GeV protons from the Main Injector 
arc directed every two seconds at a Nickel target, producing one 8 GeV antiproton 
for every 0(50, 000) incident protons, in total 0(10^) per pulse. These antiprotons 
are focused into a beamline using a Lithium lens and separated from the other 
produced particles with a pulsed dipole magnet used as a charge-mass spectrom- 
eter. Once transferred to the "Dcbuncher" ring, the large momentum spread of 
the antiprotons is reduced using radio-frequency bunch rotation [198] and stochas- 
tic cooling [199], before the beam is passed to the "Accumulator" ring, where the 
antiprotons arc collected ("stacked") and cooled further. For collider operation, ap- 
proximately 30,000 such cycles are needed. The Debuncher and Accumulator are 
both 8 GeV rounded, triangle-shaped, concentric synchrotrons with circumferences 
of 505 m and 474 m, respectively, and together with the target station are referred 
to as the "Antiproton Source" . 

To operate at optimal stacking rates, antiprotons are transferred every few hours 
from the Accumulator to the "Recycler" (an 8 GeV storage ring housed in the Main 
Injector tunnel), providing both stochastic and electron cooling [200] and thereby 
improved beam quality. As the name implies, the Recycler was originally planned 
to reuse antiprotons from the Tevatron, but this was abandoned in favor of large 
stashes (6 • 10^^) of antiprotons of high beam quality [201]. The 8 GeV antiprotons 
from either Accumulator or Recycler arc accelerated in the Main Injector to 150 GcV 
for injection into the Tevatron, where they together with the protons are ramped 
up to 980 GeV for collisions. The bunch spacing of 396 ns corresponds to a collision 
rate of 2.5 MHz. With 36 out of 53 bunches being filled, the average rate is reduced 
to 1.7 MHz. 

A typical store at the Tevatron contains about 9 • 10^° antiprotons and 26 • 10^° 
protons per bunch. Characteristic r.m.s. bunch dimensions are 45 cm (50 cm) in 
the longitudinal and 16 fim (28 fim) in the transverse direction for antiprotons 
(protons) [29], respectively. A store lasts typically 16 to 24 hours for data tak- 
ing (governed by the remaining instantaneous luminosity versus the one achievable 
with a new store) before the beam is dumped and the Tevatron refilled (within two 
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Fig. 12: Left: Peak luminosities achieved at Run II Tevatron versus time [202]. Right: 
Integrated luminosity delivered by the Tevatron and recorded by the DO experiment 
in Run II versus time [203]. 



to three hours). The increased antiproton stacking rates achieved recently provide 
shorter overall turnaround times and store durations while raising initial luminosi- 
ties, thereby enabling the maximization of the delivered luminosity per time period. 

The Tevatron is performing well and keeps setting new world records on peak 
luminosity at a hadron collider. As of July 2008, the record is 3.2-10'^'^ cm^^s^^ [196]. 
For comparison, during Run I, from 1992 to 1996 at a center of mass energy of 
1.8 TeV, the record peak luminosity was 0.2 • 10^^ cm^^s"^, and both experiments 
recorded an integrated luminosity of wO.l fb^^, respectively. As illustrated in Fig. 
I12| since the beginning of Run II in 2001, each experiment has recorded more 
than 4 fb~^, and up to half of the total Run I luminosity is now collected by the 
experiments in one single week. The analyses discussed in this review utilize datasets 
up to an integrated luminosity of 2.8 fb^^. 

Until the currently scheduled end of Run II in October 2009, the Tevatron is ex- 
pected to deliver more than 6 fb~^ to each experiment, with possible improvements 
on that value crucially dependent on the achievable antiproton stacking rates [195]. 
An extension of Tevatron running into 2010 is currently being discussed and could 
increase the integrated luminosity by an additional 2 fb~^. 



4.2. The collider experiments 

Both general-purpose detectors CDF and DO follow the generic layout of a collider 
detector in having their subdetectors arranged symmetrically in layers around the 
interaction point and beam pipe (see Fig. llSp . The inner detectors are arranged 
in concentric cylindrical layers, with charged-particle tracking systems of low mass 
surrounded by solenoidal magnets defining the core. These are enclosed by electro- 
magnetic and hadronic calorimeters that provide energy measurements and identifi- 
cation of electrons, photons and hadrons. The outer systems are dedicated to muon 
detection, relying on the penetration capabilities of muons past all inner detectors. 
Both detectors use a right-handed coordinate system with the origin at the center 
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Fig. 13: Cross section views of the CDF detector (top, [204]) and the DO detector 
(bottom, [205]). 



of the detector and the z-axis pointing along the direction of the proton beam. 
The transverse plane is spanned by the y-axis, which points vertically upwards, 
and the a;-axis that points away from the center of the Tevatron. Positions in the 
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transverse plane are frequently described using the aziniutlial angle cf) with respect 
to the X-axis, </> = arctan^, and radius r = + y'^- Based on the polar angle 
6 relative to the ^;-axis, the pseudo-rapidity rj is defined as 77 = — ln(tan|). For 
niassless particles (or in the ultra-rclativistic case where masses can be neglected), 
the pseudo-rapidity is equivalent to the rapidity y = \ \d{{E +Pz)/{E — pz)], which 
is additive under parallel Lorentz transformations, resulting in Lorentz-invariant 
rapidity differences Ay. The distance of two objects in the r] — (j) plane is usually 
denoted as Ai? = \/ l^rf + IS.cj)^. To differentiate between variables calculated with 
respect to pp collision point and the center of the detector, the latter is often denoted 
with a subscript "det" to indicate origin of the detector coordinate system. In this 
review, unless indicated to the contrary, r] is used to refer to r?dct- 

Enclosing the luminous region which exhibits a Gaussian width of approximately 
25 cm, both CDF and DO have their innermost silicon microstrip trackers. These 
provide vertexing and tracking capabilities extending to pseudo-rapidities of |?7| < 2 
and 1 77! < 3 for CDF and DO, respectively. CDF complements its tracking system 
with a cylindrical open-cell drift chamber that provides 96 track measurements for 
|?7| < 1, while DO utilizes a scintillating- fiber tracker consisting of eight cylindrical 
layers with two overlapping 835 /xm diameter fiber doublets each, providing coverage 
for |?7| <3 1.7. Both tracking systems arc enclosed by superconducting solenoidal 
magnets that provide magnetic fields of 1.4 T (CDF) and 1.9 T (DO) along the 
beamline for measuring transverse momenta [pr) of charged particles. 

Supplemental particle identification systems are placed inside and also outside 
of the magnet for the CDF and DO detectors. Within the magnet, CDF employs a 
Time-of-Flight detector based on plastic scintillator panels covering \ri\ ^ 1, mainly 
to discriminate low-energetic {p < 1.6 GeV/c) charged pions from kaons (for tagging 
heavy- flavor) . Outside of the magnet, CDF uses scintillator tiles for early sampling 
of electromagnetic showers to improve electron and photon identification in the 
central detector. DO uses central {\r]\ < 1.3) and forward (1.5 < |r/| < 2.5) preshower 
detectors consisting of several layers of plastic scintillator strips to enhance electron 
and photon identification. 

Sampling calorimeters with an inner electromagnetic and an outer hadronic sec- 
tion enclose all inner subdetectors, providing energy measurement and identifica- 
tion capabilities for photons, charged leptons and hadrons. CDF uses Icad/iron- 
scintillator sampling devices covering pseudo-rapidities \r]\ ^ 3.6, while DO uses 
mainly depleted uranium (U^^*) as absorber material and liquid argon as active 
medium for nearly compensating calorimetry within I77I <, 4.2. Between the cen- 
tral and cndcap calorimcter-cryostats (1.1 < \t]\ < 1.4), layers of scintillating tiles 
provide additional sampling of showers for DO. 

The outermost subdetectors serve to identify muons. These devices are based on 
the fact that muons rarely interact or radiate, but rather traverse the calorimeter 
as minimum ionizing particles that rarely generate electromagnetic or hadronic 
showers. CDF and DO employ scintillators and drift tubes for muon detection within 
\r]\ < 1.5 and |r/| < 2, respectively. DO has in addition 1.8 T solid-iron toroidal 
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magnets between the detection layers to provide stand-alone measurements of muon 
momentum that are independent of the central tracking system. 

The luminosity for CDF and DO is measured, respectively, using Cherenkov and 
plastic scintillation counters covering 3.6 < |ry| < 4.6 and 2.7 < \ri\ < 4.4, respec- 
tively. To select events of interest from the effective 1.7 MHz bunch-crossing rate, 
both experiments employ three-level trigger systems of dedicated hardware at the 
initial levels and commercial processor farms at the later (higher) levels. Based on 
information from tracking, calorimetry and muon systems, events are recorded at a 
rate of approximately 100 Hz for storage and further processing. 

More detailed descriptions of the CDF and DO detectors can be found in Refs. 
[204,206,207] and [205], respectively. 

4.3. Object reconstruction 

To analyze top quark events and to study properties of the top quark, the funda- 
mental objects resulting from top quark decays must first be reconstructed. This 
section gives a brief overview of the objects to be considered, and how they can be 
reconstructed in the detectors. More information on such reconstruction, specific 
to the CDF and DO experiments, can be found, for example, in Refs. [22,208]. As 
will become clear, analyses of the top quark utilize all detector components and 
therefore need a thorough understanding of their performance and calibration. 

4.3.1. Primary vertices 

The point of the primary hard scatter is referred to as the primary vertex and is de- 
termined through a fit of well-measured emerging tracks and beamline constraints 
to a common origin. With increasing luminosity, the average number of interac- 
tions per bunch crossing increases as well, leading to the reconstruction of multiple 
primary- vertex candidates, only one of which will generally be compatible with the 
hard interaction of interest. The selection of the primary vertex can be based, for 
example, on the presence of an energetic lepton, the (maximal) scalar sum of pt 
of associated tracks, or the (lowest) compatibility with being a "minimum-bias" 
interaction, based on track templates [22,207]. 

The primary vertex is the origin of all objects produced in the interaction, 
including those from top quark and their subsequent W boson decays, both of 
which cannot be separated from the primary vertex within the detector resolution 
due to the extremely short lifetimes (see also Section 17. 2p . The primary vertex is 
also used as the origin of coordinates for evaluating kinematic variables of the pp 
collision. 

4.3.2. Charged leptons 

Leptonically decaying W bosons are a source of isolated energetic charged leptons 
that can be measured well with the tracking, calorimeter and muon systems de- 
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scribed in the previous section. Such leptons are part of the event signatures for 
several top quark decay modes (see Section and key in selecting these events at 
the trigger stage. However, r leptons play a special role due to their decay charac- 
teristics. They decay leptonically 35% of the time, yielding electrons or muons (and 
two neutrinos) that on average have far lower momenta than e or /i from direct 
W iiy decay, but are otherwise hard to distinguish due to the relatively short t 
lifetime. Consequently, such decays are usually included in the event selections for 
electrons and muons which are then referred to as "leptonic" final states of the W 
boson. Decays of r into hadrons (and a neutrino) are treated separately and are 
discussed further below. 

In the context of this review, the term "leptons" refers only to electrons and 
muons, unless specified to the contrary. Their reconstruction proceeds as follows: 

• Electrons leave a track in the inner tracking system and form showers, mainly in 
the electromagnetic part of the calorimeter. These are reconstructed as clusters 
of energy deposition in the electromagnetic calorimeter and matched to recon- 
structed tracks. Further requirements include selections based on the fraction of 
energy deposited in the electromagnetic calorimeter, isolation from other energy 
depositions in a cone defined by the electron candidate, the transverse and lon- 
gitudinal distribution of the shower, and the E/p ratio of cluster energy and the 
reconstructed track momentum. 

• Muons are identified via their characteristic penetration and minimum ionization 
along their path in the calorimeter. They are reconstructed by matching central 
tracks to track segments in the outer muon system and must be consistent with 
originating from a primary vertex. Cosmic-ray muons are suppressed via timing 
requirements. Isolation of muon trajectories can be imposed both on the charged- 
track candidate and within the calorimeter. 

The misidentification of isolated leptons has different origins for electrons and 
muons. Assuming contributions from hadrons "punching through" the calorimeter 
are negligible, false isolated muons arise mainly from muons emitted in semileptonic 
decays of heavy quarks where an associated jet (see below) is not reconstructed or 
the muon emerges outside of the cone of the jet. While such semileptonic heavy- 
flavor decays also contribute false isolated electrons, significant contributions arise 
here as well, for example, from jets with large electromagnetic fractions that can 
mimic electrons, photon conversions to e^e~, or photons that overlap a random 
track. Such instrumental background processes are usually estimated from data, as 
realistic simulation of these effects is quite difficult. 

The energy scale and resolution for leptons can be assessed by studying the 
reconstructed mass of the Z boson in Z ££ events. Z boson decays are also useful 
for studying lepton identification efficiencies with the "tag and probe" method in 
which one lepton is required to be well- identified (the "tag"), thereby providing a 
reasonably pure sample of Z bosons, while the second lepton serves as a probe for 
the efficiency being studied. 
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Fig. 14: Illustration of the evolution of a calorimeter jet from an initial parton [209]. 
The dashed lines represent the jet cone. 

4.3.3. Quark and gluon jets 

The hadronization of quarks and gluons leads to collimated showers of hadrons, 
referred to as "jets". The jet axis is highly correlated with the original parton's 
direction. While it is not possible to differentiate between quark and gluon jets on 
a per-object basis, they can be distinguished on a statistical basis because of small 
differences in shower shape (gluon jets tend to be wider and contain more "soft" 
particles). 

Jets are reconstructed from their energy depositions in the cells of the calorimeter 
using cone algorithms [210,211] that combine cell energies within a cone of fixed 
radius Ai? (see Fig. [T4t . The size of the cone radius is a compromise between 
collecting a high fraction of the original parton's energy and resolving the energy 
depositions of close-by partons, especially in busy tt events. DO uses a cone size of 
AR = 0.5, while CDF uses AR = 0.4. 

The measured jet energies are converted into particle-level energies through jet 
energy scale (JES) corrections [212,213] that take into account effects such as the 
presence of energy depositions not originating from the hard scattering process 
(noise in the calorimeter, multiple interactions, . . . ), particles within the jet cone 
that deposit energy in the calorimeter but outside of the cone or vice versa, due to 
their curved trajectories in magnetic fields and showering effects, and the calorimeter 
response accounting, for example, for nonlinearities or energy loss in uninstrumented 
regions of the detector. 

The electromagnetic calorimeter scale determined from resonances such as Z — s- 
ee, as described above, can be transferred to the full calorimeter by requiring a 
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balance in pT in photon + jet events. The intercalibration of the calorimeter is then 
complemented with dijet events. A precise JES calibration is a challenging task 
involving highly complex procedures to ensure understanding of all contributions 
and their systematic uncertainties. DO has achieved a JES precision at the 1-2% 
level over a wide kinematic range [213]. For this level of precision, strictly speaking, 
the JES is only applicable to photon + jet samples, and additional uncertainties 
need to be taken into account when transferring the JES, for example, to top quark 
samples [214]. The first direct measurement of the JES for h quarks at the Tevatron 
based on Z ^ 66 events has recently been performed by CDF, reaching a precision 
of better than 2% [215]. 

4.3.4. r jets 

T leptons decay into hadrons (and a neutrino) 65% of the time, with Ri77.5% of 
these decays yielding a single charged particle ("1-prong decays") and ~22.5% 
three charged particles ( "3-prong decays" ) [29] . These "hadronic" r decays are re- 
constructed as jets (that often also contain tt" mesons) and can be discriminated 
statistically from quark and gluon jets via their narrow shower shape and low track 
multiplicity within the jet cone [216,217]. 

4.3.5. b jets 

The identification ( "tagging" ) of 6 jets is a very powerful tool for separating the top 
quark signal from its background processes, which typically exhibit little heavy- 
flavor content. Also, the combinatorics for reconstructing top quark events from 
their final-state objects can be reduced using this additional information. There are 
in general two different approaches to identifying B hadrons formed from 6 quarks: 

• Lifetime Tagging: Due to their lifetime of about 1.5 ps and the boost from 
top quark decay, B hadrons can travel several millimeters before they decay. The 
resulting charged particle tracks therefore originate from (point to) a vertex dif- 
ferent than the primary one. This can be exploited by searching for secondary ver- 
tices significantly displaced relative to the primary event vertex (secondary-vertex 
tagging) or by requiring significant impact parameters relative to the primary 
vertex for tracks, without reconstructing a secondary vertex {impact-parameter 
tagging). A probability can also be calculated for a jet to come from the primary 
vertex based on the impact parameters of all its associated tracks (jet-probability 
tagging), or a combination of all information from the above tagging algorithms 
into a neural network response can be used (NN tagging). The two latter methods 
yield continuous output variables that can be used as input for further multi- 
variate analysis or for selecting analysis-specific values as compromises between 
6-tagging efficiency and the fraction of light-quark jets that are misidentified as 
6 jets. 
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• Soft-Lepton Tagging: This tagging is based on semileptonic decays of b and 
c hadrons with branching fractions of and 10%, respectively. With two b 

quarks and two W bosons per ti decay, and the fact that about one third of 
the W boson decays yield charm quarks (cs), the fraction of events containing a 
soft (low pt) lepton in a jet is about 40% per lepton flavor (e, fi). The isolation 
criteria used for leptonic W boson decays do not work for these leptons, and their 
reconstruction within jets is quite challenging, especially for electrons. 

Although the mistag rate for lifetime-tagging is usually very small for light-quark 
{u, d, s) and gluon jets, this is not the case for charm jets. For example, a typical 
operating point for DO's NN tagger yields a 6-tag efficiency of «50% and a mistag 
rate for light jets of ~0.5%, while it is ~10% for c jets [214]. More information 
on 6-tagging algorithms and their application in top quark analyses, including the 
performance for 6 jets, c jets and light-quark or gluon jets, can be found in Rcfs. [207, 
208,218] for CDF and [22,214,219] for DO. 

4.3.6. Neutrinos 

Neutrinos are not detected directly because of their negligible interaction cross 
section. Since the energy component along the beam axis at a hadron collider is 
unknown, only the transverse momentum carried away by neutrinos (or any other 
"invisible" particles) can be inferred from momentum conservation in the transverse 
plane. This "missing" transverse momentum (^t) is calculated from the vector sum 
of transverse energy depositions in the calorimeter, corrected for the energy scale 
of reconstructed electrons and jets and for the momenta of reconstructed muons 
(corrected for energy loss in the calorimeter). The resolution therefore depends 
strongly on the other objects present in the event. Taking this into account, for 
example, through selections on the significance of jix rather than its absolute value, 
improves performance. 

4.4. Top quark event signatures in the standard model 

Having discussed the reconstruction of the fundamental objects from the initially 
occurring particles, the experimental signatures of top quark events can be exam- 
ined. As noted in Section I3.3.1[ in the framework of the standard model the top 
quark decays dominantly into a W boson and a b quark. Consequently, the observed 
final states are defined by the decay modes of the W boson. 

W bosons decay into two fermions, either leptons (a charged lepton-neutrino 
pair iDe,^ = e,/i,r), with equal probability per lepton flavor at lowest order in 
perturbation theory, or into quark-antiquark pairs qq' with q = u,c and q' = c?, s, b. 
At Born level, the "hadronic" decay widths are enhanced over the leptonic modes 
by a color factor of three (taking the three possible quark colors into account), 
and scaled by the appropriate squared CKM matrix element Similar to 

top quark decay, the off-diagonal CKM matrix elements are greatly suppressed. 
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Fig. 15: Illustration of the SM tt decay modes and their branching fractions via the 
possible W boson decays. For each decay mode, a hh quark pair from ti decay is 
also present. 

and therefore only the ud and cs decay modes are considered in the following, 
contributing approximately 95% of the hadronic decay width [29] . 

In summary, W bosons decay leptonically with a branching fraction of wl/9 per 
lepton flavor and ^1/3 for each of the hadronic decays {ud and cs). The resulting 
possibilities for tt decays are illustrated in Fig. [151 where also the nomenclature for 
the different decay modes is introduced. These branching fractions do not yet take 
account of the leptonic decays of r leptons and their contributions to final states 
involving electrons and muons, as discussed in Section [4.3l This is considered in the 
following description of the four basic ti event classes: 

(i) Dilepton channels: Both W bosons decay leptonically {£vi^ £ = e, resulting 
in a final state comprised of two isolated high-p^ leptons, corresponding to 
the two neutrinos, and two b jets. Including leptonic r decays, this channel has 
a branching fraction of approximately 6.5%, shared wl:l:2 by the ee, /z^ and 
e/x final states. While these channels give samples of highest tt signal purity, 
they suffer from limited statistics due to the small branching fractions. 

(ii) Lepton + jets channels: One W boson decays leptonically and the other one 
hadronically, yielding a final state containing one isolated high-p^ lepton, jix, 
and four jets. Including leptonic r decays, these channels exhibit a branching 
fraction of approximately 34.3%, shared about equally by the e + jets and fi + 
jets final states. These channels represent the best compromise between purity 
of sample and available statistics. 

(iii) All-hadronic channel: Both W bosons decay to qq' pairs, resulting in a six-jet 
final state. With a branching fraction of «45.7%, this channel yields the high- 
est statistics of tt events but also suffers from large background from multijet 
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production. 

(iv) Hadronic r channels: Final states where at least one W boson yields a 
charged r lepton that in turn decays into hadrons (and a neutrino) are called 
hadronic r channels, and together comprise a branching fraction of wl3.5%. 
The nature of the decay of the second W boson is used to differentiate between 
the r + jets, r + lepton, and rr final states, which contribute with 9.5%, 3.6% 
and 0.5% branching fractions, respectively. The corresponding experimental 
signature has four/two/two jets, ^t, one/one/two hadronic r final states, and 
no/one/no isolated high-p^ lepton. The identification of r ^ hadrons makes 
these final states especially challenging to reconstruct. More inclusive sample 
selections, requiring, for example, leptons and isolated tracks or fix and (&- 
tagged) jets provide significant fractions of r — s- hadrons events, without their 
explicit reconstruction. 

In all of the above final states, two of the jets are b jets from tt decay. The rr 
final state of ti production remains the only channel that has not yet been explicitly 
analyzed. All others are discussed in Section [5T] Properties of top quarks have been 
extracted mainly from the first three of the above channels, and especially from the 
lepton + jets channel. 

Full kinematic reconstruction of tt events is possible in the all-hadronic final 
state since there are no high-pr neutrinos present. In the lepton + jets channel, a 
twofold ambiguity arises from the determination of the neutrino Pz when constrain- 
ing the invariant mass of the lepton and missing neutrino to mw, while the dilepton 
channel is kinematically underconstrained because two neutrinos contribute to ^t- 
The unknown assignment between partons and reconstructed objects in tt events 
leads to various possible combinations. The combinatorics can be reduced through 
identification of 6 jets. In particular, when both 5 jets are identified, only four combi- 
nations remain to be considered in the lepton -I- jets channel (including the neutrino 
Pz ambiguity), and six combinations in the all-hadronic channel. 

The experimental signature for single top quark production is based on the top 
quark decay mode and the production channel: in the s- (tb-) channel, the top quark 
is produced with an additional b jet, while in the t- (tqb-) channel a forward light- 
quark jet accompanies top quark production, sometimes along with another b jet 
from the gluon splitting into bb (see Fig. |8]). The W boson from top quark decay 
is usually required to decay leptonically {Ivijl = e, fj.) so as to suppress multijet 
background. Consequently, the final state signature of single top quark production 
contains an energetic isolated electron or muon, -^t and at least two jets, with at 
least one of them being a b jet. 

The large mass of the top quark makes it less likely that it is produced with 
large kinetic energy at the Tevatron. Its decay products therefore tend to be emitted 
at central rapidities, non-planar with good angular separation, and are character- 
ized by a large sum of transverse energies Ht- Event selections usually require 
the channel-characteristic objects (leptons, ^t, and (6 tagged) jets) to be present 
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with energies typically greater than 15 to 20 GeV. Apart from selections on data 
quality to ensure a well-performing detector and specific trigger selections, a well- 
reconstructed primary vertex in the central detector region is also required. Vari- 
ations in observed jet multiplicities are possible as well due to, for example, jet 
reconstruction thresholds, jet splitting and merging during reconstruction, and ad- 
ditional gluon jets (from initial- and final-state radiation). 

More details on event selection, contributions from background, and sample 
compositions in different analyses (including those concerned with non-standard- 
model signatures), are given in the following chapters. 

4.5. Monte Carlo generation 

A reliable and well-understood Monte Carlo (MC) simulation of signal and back- 
ground processes is a crucial ingredient for any top quark analysis. MC samples are 
needed to understand detector response and acceptance, selection efhciencies and 
kinematic distributions of physical variables and their normalizations. This requires 
both good modeling of the production process from the parton to the hadron level, 
and an accurate simulation of detector response to signatures of interest. 

MC simulations of hadron interactions are based on the factorization theorem 
discussed in Section 13.11 splitting up hadron collisions into universal long distance 
(small Q^) phenomena and perturbatively calculable short distance phenomena. A 
generic example of steps in the simulation of a hadron collision is illustrated in Fig. 
[TBI and described below. 

The non-perturbative Parton Distribution Functions (PDFs) describe the distri- 
bution of the proton (or antiproton) momentum among its partons. The interaction 
of the incoming partons in the hard process of interest is then evaluated based on 
fixed-order (usually LO) matrix elements, yielding the outgoing partons and their 
characteristics, such as their momenta and colors. The ensuing parton showering 
adds higher-order effects through parton splitting into pairs of softer partons (gluon 
radiation, gluon-splitting, photon radiation,. . . ), until non-perturbative hadroniza- 
tion sets in at low , forming color-neutral hadrons from the colored partons, 
which is based on phenomenological models. Unstable particles and resonances are 
then made to decay into the final remnants. 

The colored beam-remnants of the proton and antiproton, other soft multi- 
parton interactions, and color connections to the hard process are added, and all 
form the "underlying event" . Additional soft proton-antiproton collisions from the 
same colliding bunch (minimum-bias events) have to be superimposed in the MC, 
and depend on the instantaneous luminosity. Finally, any overlapping interactions 
from consecutive bunch crossings "leaking" into the current event (pile up) must 
also be considered. 

In principle, different programs and models can be used for every stage in the 
above process, and the best choice may depend on the process to be studied, which 
illustrates the complexity of these MC simulations. Also, the parameters in some 
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Fig. 16: Illustration of a hadron collision [220], indicating some of the steps involved 
in MC simulation. 

models have to be tuned to distributions in data before they can provide adequate 
descriptions of interactions [221]. General purpose, complete event generators such 
as HERWIG [222] and pythia [223] include LO matrix elements for a variety of 
processes, as well as models for showering and hadronization of partons. They are 
widely used, either stand-alone or in combination with other generators that sim- 
ulate the hard process and pass on the information to implement showering and 
hadronization. An introduction and overview of available MC generators can be 
found in Ref. [224], and specific simulation tools for top quark production and de- 
cay are reviewed in Ref. [225]. 

CDF and DO have implemented different simulation chains for their top quark 
analyses, and certain analyses use variants on what is described here. CDF uses 
CTEQ5L [157] PDFs for its generators, while DO employs CTEQ6L1 [94] PDFs. 
The ti signal is generated either with pythia v6.2 [226] (CDF) or ALPGEN v2.1 [227] 
interfaced with pythia v6.3 [228] for parton showering (DO). The latter employs a 
jet-parton matching algorithm to avoid double counting of final states that could 
otherwise be populated from both the hard process and parton showering [229,230]. 
CDF also uses herwig v6.4 [231] for systematic cross checks of modeling of signal. 
For single top MC, DO uses SingleTop [232] based on CompHEP [233], interfaced 
with pythia, while CDF utilizes MadEvent [234] and MadGraph [235], also 
interfaced with pythia. Both experiments approximate t-channel production at 
NLO through a combination of contributing 2 — > 2 and 2^3 processes. Any 
signal cross sections obtained at LO are scaled up to match higher-order theoretical 
prescriptions (see Sections 13.11 and 13. 2p . Most analyses use a top quark mass of 
175 GeV/c2. 

For simulation of VF-|-jets and Z-l-jets background processes, both collabora- 
tions utilize ALPGEN, interfaced with herwig (CDF) or pythia (DO) for parton 
showering, and both apply the above-mentioned jet-parton matching technique. 



Top Quark Measurements at the Tevatron 41 



ALPGEN generates higher final-state parton muhiphcities from 2 ^ n processes, 
and therefore large jet multiplicities, based on exact LO matrix elements, including 
the production of heavy-flavors, which is especially important for analyses using h 
tagging. For the decays of r leptons, both collaborations use TAUOLA [236,237]. 
For h and c hadron decays EvtGen [238] and QQ [239] are used, the latter only by 
CDF. Efii'ects of additional minimum-bias events and pile up are based on pythia at 
CDF, and DO uses zero-bias collider data taken by randomly sampling filled bunch 
crossings, overlaid with the simulated events. 

The generated events are propagated through detector simulation based on 
GEANT [240], which contains a full description of positions, geometrical shapes, 
and types and amounts of material comprising the detectors. The particles in an 
event arc tracked through the detector volume, where they encounter energy loss 
and multiple scattering that depend on particle type, the traversed material and the 
particle momenta, and undergo decay corresponding to their lifetimes. The response 
of the detector's readout electronics to these interactions, including noise and inef- 
ficiencies, is then obtained in a digitization step that yields initial simulated event 
characteristics that are processed with the same reconstruction chain as collider 
data. 

It is not easy to obtain good agreement between any simulation and experimental 
data. For example, object reconstruction, identification and selection efficiencies 
tend to be higher in MC compared to data, and must be corrected via scale factors. 
Such scale factors arc usiially derived from comparison of efficiencies in control 
samples of simulations and data in final states such &s Z ^ (.(. for leptons, and 
7-|-jets for jets. Scale factors can be parametrized in terms of variables sensitive to 
these corrections. Also, energy scales and resolutions for reconstructed objects must 
generally be adjusted. Certain effects are hard to simulate, so that at times only 
data can be used to provide rates for, e.g. jets to mimic lepton signatures. 

Before searching for any signal, the background model must be verified using 
data in control samples that are sufficiently depleted from possible signal, as, for 
example, by requiring a reduced jet multiplicity or no 6 jets to be present. Sometimes 
the shape or normalization of differential distributions has to be corrected, which 
usually reflects not optimally tuned MC or insuflicient precision in the model. 

With increasing Tevatron luminosity, data-based constraints can improve the 
understanding of dominant background processes such as vector boson -|- (heavy- 
flavor) jet production, both in terms of shapes and normalization. This can benefit 
MC simulations and the precision of measurements. For example, the production 
of W bosons has been investigated in terms of associated jet production [241] and 
compared with LO and NLO predictions, or associated heavy-flavor production was 
compared with standard model expectations and found to be in agreement [242, 
243]. Nevertheless, dedicated studies of W boson + c [244,245] or b jet [246,247] 
production have been performed as well, with the most recent results indicating that 
the production rates for these processes are currently underestimated by alpgen. 

A more detailed overview of MC simulations used in top quark analyses at the 
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Tevatron, both for signal and background processes, and the remaining challenges 
can be found in Refs. [248,249]. 

5. Measurements of Top Quark Production 

In this chapter, measurements of top quark production both via the strong and elec- 
troweak interactions are described. Observed rates and mechanisms of production 
are compared with the standard model expectations and used to derive constraints 
on specific extensions of the standard model impacting the properties under con- 
sideration. 

5.1. Top quark pair production cross section 

Measurements of the ti production cross section are important for several reasons. 
They provide a powerful test of the predictions of perturbative QCD calculations at 
high transverse momenta. As shown in Section 13.11 the uncertainties on the tt rate 
predictions have reached the level of «10%, a precision which has already been 
matched by the measurements performed at the Tevatron. 

Deviations from the standard model prediction could arise, for example, from 
novel production mechanisms such as a new resonant production mode in addition 
to the standard model one as discussed in Section 15.71 New physics contributing to 
the electroweak symmetry breaking will probably couple to particles proportional to 
their mass, making the top quark and its strong coupling reflected in its production 
rate a highly interesting probe for such effects. Different top quark decay modes, 
such as the decay via a charged Higgs boson competing with SM decay, as examined 
in Section [6. 51 could cause apparent different production rates amongst the various 
decay channels via modified branching fractions. Contributions of new physics to 
the background samples in the various channels could have similar effects. Analyz- 
ing different decay channels consequently helps not only to improve statistics and 
studies of properties of top events, but is also a sensitive probe for physics beyond 
the standard model. 

To extract a cross section requires good understanding of the reconstruction and 
identification of the involved objects and of the modeling of the contributing back- 
ground processes. By providing selections for samples enriched in top quark signal 
and of well characterized composition, cross section analyses form the foundation 
of all further top quark property analyses. 

Top quark pair production has been studied by now in all possible decay modes 
- the dilepton, lepton -I- jets, all-hadronic and hadronic r channels as defined in 
Section 14. 4( with the tt decay mode being the only exception due to marginal 
branching fraction and challenging separation from background processes. As men- 
tioned before, in the context of this review the term leptons refers to electrons and 
muons alone unless indicated to the contrary. 

The event selections usually require the presence of characteristic objects from 
the top quark decay, namely, leptons, ^t, (heavy-fiavor) jets, or hadronically decay- 
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ing T, with energies typically exceeding 15 to 20 GeV. Further kinematic character- 
istics can be exploited to separate the top quark signal from the various background 
processes. Due to the large mass of the top quark, its decay products tend to be 
very energetic, emitted at central rapidities and non-planar with good angular sep- 
aration. In contrast to this, the jet energy spectrum for background processes with 
jets from gluon radiation is steeply falling. The observed objects are emitted less 
isotropically but more back to back, and mismeasured objects giving rise to tend 
to exhibit characteristic angular correlations with the reconstructed ^t- 

Consequently, additional variables available for top quark signal selection are 
based on the energy present in the event, such as the sum of transverse energies 
(Ht), or the invariant mass of a combination of reconstructed objects. Event shape 
variables such as sphericity and aplanarity, derived from the eigenvalues of the nor- 
malized momentum tensor of the objects considered [250], or centrality, defined as 
the ratio of Ht and the sum of the objects' energies, provide additional discrimi- 
nation. Furthermore, angular relations between reconstructed objects (for example 
A(/)(^T,-^)) and single-object kinematic quantities (such as the jet of highest (lead- 
ing) transverse momentum) are frequently used as well. 

Depending on the ti decay mode considered, the use of b tagging in its different 
forms (see Section H75|) is optional for the event selection. In the dilepton and lepton 
-|- jets channels selections based purely on topological and kinematic characteristics 
suffice for a good signal to background ratio (S/B). Adding b tagging improves sam- 
ple purities but also implies a stronger model dependence by relying on b quarks to 
be present in the final state. The actual extraction of the signal fraction proceeds 
either in a counting experiment or via template fits using the full shape informa- 
tion of the sensitive variable under consideration. While the latter is usually more 
sensitive, it also exhibits a stronger dependence on the MC modeling. Using differ- 
ent methods with different systematic uncertainties to measure the same quantity 
provides a way to assess the model assumptions from different perspectives and to 
check internal consistency, and is beneficial for combinations of increased precision. 
Non-overlapping (orthogonal) sample selections facilitate later combinations of re- 
sults as independent measurements by removing the need to evaluate the correlation 
amongst the measurements from ensemble tests. 

Once the sample composition is measured, the it production cross section is 
calculated as follows: 

-^observed -^background 

^" ^ sB JCdt ' 

where iVobsorved (-^background) is the total (background) number of events, e is the 
it selection efficiency, including detector acceptance, B is the branching fraction for 
the ti decay mode in question and J Cdt is the integrated luminosity of the used 
dataset. 

As illustrated in Section 13.11 the it cross section depends on the top quark 
mass, decreasing by about 0.2 pb for each GeV/c^ increase in the rat mass range 
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from 170 GeV/c^ to 180 GeV/c^. The ti cross section results given in the following 
sections generally refer to a top quark mass of 175 GeV/c^; the few cases where a 
top quark mass of 178 GeV/c^ was assumed will be pointed out explicitly. Especially 
for the recent measurements, a parametrization of the obtained result versus top 
quark mass is provided to allow easy projection to the current world-averaged top 
quark mass. The cross section dependence on the mass can also be turned around 
to provide a measurement of the top quark mass from the cross section, which will 
be discussed further in Section [7.3.41 

In the following subsections the published and latest preliminary Run II results 
will be referenced for the different tt decay modes. Some analyses will be high- 
lighted in more detail. The agreement with theoretical predictions is illustrated in 
the summary section, where results of combinations across channels are also given. 

5.1.1. Dilepton final state 

A typical ti dilepton event selection requires two isolated high-pT leptons of opposite 
charge, and at least two central energetic jets. The dominant physics background 
processes exhibiting both real leptons and arise from diboson {WW, ZZ, WZ) 
production and from Z/7*-|-jets processes, with Z/j* — )■ r+r^,T e, ii. Mis- 
reconstructed "^T caused by experimental resolution in Z/7*+jets events (with 
Z/-f* e+e~/^+^^) contributes to the instrumental background, as does M^-|-jets 
and multijet production (often called QCD), where one or more jets mimic the 
isolated-lepton signature. The physics background processes are usually modeled 
using Monte Carlo simulation, while instrumental background processes (especially 
those involving false isolated-lepton signatures) typically are estimated from data. 
The signal purities of the resulting samples are usually quite good, with a signal to 
background ratio (S/B) typically better than two. 

The sample purity can be further enhanced by means of additional kinematic 
requirements such as Ht to be above a certain threshold, or by selecting events with 
at least one identified b jet. However, this reduces the already limited statistics in 
this channel. To increase the signal acceptance, the reconstruction and isolation 
requirements on the second lepton can be relaxed. If the second lepton is only 
required to be reconstructed as an isolated track (termed lepton -I- track selection), 
especially 1-prong hadronic r decays can then also contribute to the signal. 

In a recent preliminary analysis, CDF determines the it cross section from a 
2.8 fb~^ dataset by requiring two oppositely charged reconstructed isolated leptons 
with Et > 20 GeV, > 25 GeV and at least two jets within |?7| < 2.5 and 
Et > 15 GeV, with the leading jet fulfilling Et > 30 GeV. The tt cross section 
is extracted from the resulting sample once without any additional cuts and once 
after increasing the purity by requiring at least one of the jets to be b tagged. 
The background from Z/7* and diboson WW, ZZ, WZ events is derived from MC, 
while false isolated-lepton signatures are estimated from a dilepton dataset where 
both leptons have the same charge (same sign, "SS"), assuming their contribution 
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Fig. 17: Distributions in jet multiplicity (left) and number of 6 tags (right) in tt dilep- 
ton candidate events, respectively, for > 1 and > 2 jets in 2.8 fb~^ of data analyzed 
by CDF [251]. 



is identical in the opposite sign ("OS") signal selection [251]. 

The untagged sample yields 162 events with a total background contribution 
of 51.9 ± 4.5, where the dominant uncertainties arise from the estimate of false 
leptons and the uncertainty on the jet multiplicity correction factors applied in the 
MC. Requiring at least one h tagged jet, 80 events remain, with an expected total 
background of 4.0 ± 1.7, with the dominant uncertainties arising again from the 
estimates of false leptons and also from uncertainties on the h tag modeling. The 
sample composition is illustrated in Fig.[T71 and detailed in Table [5l The extracted 
cross sections are given in Table [SI together with the other dilepton channel results 
obtained thus far in Run II. 



5.1.2. Lepton + jets final state 

A typical tt lepton + jets event selection requires exactly one isolated high-pT lep- 
ton, and at least three central energetic jets, allowing both lepton + jets and 
lepton + hadronic r signatures to contribute. The dominant physics background in 
this final state arises from W boson + jets production, and the main instrumental 
background comes from QCD multijet production with a jet mimicking the isolated- 
lepton signature. Additional smaller background contributions arise from Z/j* + 
jets, diboson and single top production. While for these smaller background pro- 
cesses shape and normalization are determined commonly from simulation and NLO 
cross sections, W + jets events are usually normalized to data, and their differential 
distributions derived from Monte Carlo. The QCD multijet background's shape and 
normalization are typically derived from data, using, for example, datasets fulfilling 
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Tabic 5: Event yields and sample composition after tt preselections in CDF's 2.8 fb 
dilepton dataset before and after requiring at least one b tagged jet [251]. 



Process 


e+e 


before tag 


-ging 


i+i- 


b tagged 
£+£- 


tt, a = 6.7 pb 


29.2 


21.5 


59.9 


110.6 


65.2 


Z/j* 


9.25 


4.79 


0.52 


14.6 


0.78 


Z/j* ^ T+T- 


2.84 


2.55 


6.62 


12.0 


0.60 


ww e+e~ 


3.05 


2.03 


5.07 


10.2 


0.44 


WZ i+£- 


1.52 


0.72 


0.67 


2.91 


0.09 


zz e+£- 


0.80 


0.40 


0.26 


1.46 


0.10 


Totals (MC) 


46.7 


32.0 


73.0 


151.7 


67.2 


±± (SS) Data 


3.81 


0.00 


6.96 


10.8 


2.00 


Sum Expected 


50.5 ± 1.7 


32.0 ± 1.3 80.0 ± 2.5 


162.5 ± 4.5 


69.2 ±1.7 


+- (OS) Data 


54 


33 


75 


162 


80 



Table 6: ti cross section measurements in dilepton final states performed thus far 
at the Run II Tevatron with their integrated luminosities, data selections {££ = 
dilepton, £+trk = lepton + track) and analysis methods used. The first three results 
have been published; the others are preliminary. The measurement marked with an 
asterisk refers to to* of 178 GeV/c^ rather than 175 GeV/c^, and, unlike the other 
results, incorporates the luminosity uncertainty within the first given uncertainty. 



J£dt 
[fb-i] 


Sel. 


b tag 


(Ttt±(stat.)±(syst.)±(lumi.) 

[pb] 


Ref. 


0.2 


££,£+tik 


no 


7.0^2-4 +1-6 ±0.4 


[252] 


0.2 


££ 


no 


8.6t^;7± 1.1±0.6 


[253] 


0.4 


££,£+tik 


no,yes 


7.4 ± 1.4 ±0.9 ±0.5 


[254] 


0.4 


££ 


no 


8.5j^2.2 io.3 (*) 


[255] 


1.0 


££ 


no 


6.811:? tS]i± 0.4 


[256] 


1.0 


££ 


no 


7.oll:ilS]:«±o.4 


[257] 


1.0 


^+trk 


yes 


5.0t}jl°J±0.3 


[258] 


1.0 


££,£+tTk 


no,yes 


6.2 ± 0.9to[^ ± 0.4 


[259] 


1.0 


£+trk 


yes 


10.1 ±1.8 ±1.1 ±0.6 


[260] 


1.1 


£+trk 


no 


8.3 ± 1.3 ±0.7 ±0.5 


[261] 


2.8 


££ 


yes 


7.8 ±0.9 ±0.7 ±0.4 


[251] 


2.8 


££ 


no 


6.7 ±0.8 ±0.4 ±0.4 


[251] 
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the complete event selection, except the tight lepton isolation, for the background 
shape, and the rate for jets to mimic leptons is derived from data for the normal- 
ization. 

Samples selected with such a basic preselection exhibit S/B ratios below unity, 
around 1/4. Signal purity can be improved significantly via additional topological 
selections or by using b tagging. When no b tagging is used in an analysis, then 
multiple topological and kinematic event properties are usually combined in a mul- 
tivariate discriminant to yield good signal to background separation without relying 
on the presence of b jets in the events, therefore being less model dependent. The 
sample composition can then be determined from a template fit in that sensitive 
variable, providing a higher sensitivity than a plain cut. 

Requiring identified 6 jets to be present in an event is a powerful tool to reject the 
background processes which exhibit little heavy-fiavor content, b tagging algorithms 
based on the long lifetime of B hadrons or reconstruction of soft leptons within jets 
that originate from semileptonic B decays, as discussed in Section |4!3| have been 
deployed for that purpose. Using b tagging, very pure tt samples can be selected 
that have S/B > 10, if at least four jets and at least two identified b jets are required 
in each event. 

The most precise cross section measurement published thus far has been per- 
formed by DO on a 0.9 fb~^ lepton -f jets dataset [262]. Events are selected by 
requiring exactly one isolated electron or muon with Et > 20 GeV, fir > 20 GeV 
for e + jets and > 25 GeV for ^ + jets, at least three jets with |?7| < 2.5 and 
Et > 20 GeV, and leading-jet Et > 40 GeV. Cuts on the azimuthal separation 
between lepton and ^t are applied to suppress background from misreconstructed 
objects. After these selections, the ti signal contributes only about 20% of the total 
sample. The ti cross section is measured using two complementary analyses. 

One approach is based on lifetime b tagging, requiring at least one jet in the event 
to be tagged and determining the ti production rate through a maximum likelihood 
fit to the observed event yields in the different subchannels defined by lepton flavor, 
jet multiplicity and b tag multiplicity. The dominant systematic uncertainties arise 
from uncertainties on tagging efficiencies for b, c, q and gluon jets and the jet energy 
calibration. The second analysis utilizes topological likelihood discriminants for the 
different subchannels based on lepton flavor and jet multiplicity. After applying an 
additional cut on jets of Ht > 120 GeV for three-jet events, five or six different 
variables such as angular object separation, sphericity and aplanarity (which pro- 
vide good discrimination power and are well modeled in MC) are combined into 
discriminants for each subchannel. The sample composition is then determined in 
a maximum likelihood fit of templates of signal and background contributions to 
the observed discriminant distributions. The dominant systematic uncertainties in 
this method arise from uncertainties on the selection efficiencies and the likelihood 
fit uncertainty derived using statistical fiuctuations in the likelihood discriminant 
template shapes. The sample compositions for both analyses are illustrated in Fig. 
[T8l and detailed in Table [71 
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Fig. 18: Topological likelihood distribution for £ + jets tt candidate events when 
three jets are required (left), and jet multiplicity distribution when at least two b 
tagged jets are required (right) in DO's 0.9 fb~^ lepton + jets dataset [262]. 



Table 7: Event yields and sample composition in DO's 0.9 fb^^ lepton + jets dataset 
for both the topological and the b tagging analysis. The tt contribution is based on 
the measured cross section in the respective analysis [262]. 





3 jets, 


3 jets. 


>4 jets, 


>4 jets, 


3 jets. 


>4 jets. 




1 tag 


>2 tags 


1 tag 


>2 tags 


topo 


topo 




147 ±12 


57 ±6 


130 ±10 


66 ±7 


245 ±20 


233 ±19 


-^othcr 


105 ±5 
27 ±2 


10 ±1 
5 ±1 


16 ±2 
8 ±1 


2 ±1 
2 ±1 


1294 ±48 


321 ±30 


-Nmultijct 


27 ±6 


3 ±2 


6 ±3 


±2 


227 ±28 


70 ±12 


total 


306 ±14 


74 ±6 


159 ±11 


69 ±7 


1766 ±59 


624 ±37 


Ndata 


294 


76 


179 


58 


1760 


626 



Both analyses exhibit a statistical correlation of 0.31, as determined from en- 
semble studies and are combined using the best linear unbiased estimate (BLUE) 
approach [276,277]. The resulting cross section is given in Table [51 together with 
the other lepton + jets channel results obtained thus far in Run II. 



5.1.3. All-hadronic final state 

To select all-hadronic ti decays typically requires at least six central energetic jets 
per event, and no isolated energetic leptons or significant to be present. The 
overwhelming background process here is QCD niultijet production, dominating 
over the signal by three orders of magnitude after online selection of events us- 
ing triggers on multiple jets and Ht in the event above a certain threshold. This 
background is usually modeled from the data, as the theoretical description of final 
states with such high jet multiplicities has large uncertainties and datasets even 
more depleted from signal can be easily obtained, for example, by selecting a lower 
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Table 8: ti Icpton + jets cross section measurements performed thus far at the Run II 
Tevatron with their integrated luminosities, data selections and analysis methods 
used. The first eleven results have been published; the others are preliminary. The 
measurements marked with asterisks (*) refer to rrit of 178 GeV/c^ rather than the 
standard 175 GeV/c^. Measurements marked with a double cross (|) include the 
luminosity uncertainty in the systematic uncertainty, while measurements marked 
with a cross (f) have the first uncertainty represent the statistical and systematic 
uncertainties combined in quadrature and the second representing the luminosity. 



fCdt 

L J 


Sol. 


b tag 


(T*f±(stat.)± (syst . ) ± (lumi . ) 

[pi'] 


Ref. 


9 




yes 


r p+1.2+0.9 

''•D-l.l -0.6 \\) 


[ZU/J 


0.2 


^+jets 


yes 




[263] 


0.2 


^+jets 


yes, soft-/x 


5.3 ±3.311:3 {%) 


[264] 


0.2 


^+jets 


no 


6.6 ± 1.1 ± 1.5 (t) 


[265] 


0.2 


^+jets 


no 


6.7t^ill:?±0.4 


[266] 


0.2 


^+jets 


yes 


8.611:^ ±0.6 (t) 


[267] 


0.3 


l+jets 


yes 


8.7±0.9tJ:i 


[268] 


0.3 


e+]ets 


yes 


8.9±i.otl:; 


[208] 


0.4 


i?+jets 


yes 


6.6 ±0.9 ±0.4 (t) 


[219] 


0.4 


£+iets 


no 


6.4tJ:3±0.7±0.4 


[269] 


0.9 


^+jets 


no, yes 


7.4 ±0.5 ±0.5 ±0.5 


[262] 


0.4 


^+jets 


yes, soft-/x 


7.3t?:°±0.4 (t) 


[270] 


0.7 


i+)ets 


yes 


8.5 ± 0.6 ± 1.0 {%) 


[271] 


1.0 


e+]ets 


yes 


8.2±0.5tg:?±0.5 


[257] 


1.7 


l+jets 


yes, soft-e 


7.8 ± 2.4 ± 1.5 ±0.5 


[272] 


2.0 


e+]ets 


yes, soft-/i 


8.7±l.ltg:^±0.6 


[273] 


2.7 


^+jets 


yes 


7.2 ±0.4 ±0.5 ±0.4 


[274] 


2.8 


^+jets 


no 


6.8 ±0.4 ±0.6 ±0.4 


[275] 



jet multiplicity. 

After preselection, signal and background are separated further by applying 
h jet identification and using multivariate discriminants based on topological and 
kinematic event properties. 

CDF has published the most precise cross section analysis in the all-hadronic 
final state to date, based on 1 fb~i of data [278]. Events are required to have at 
least six and at most eight jets with Et > 15 GeV, AR > 0.5 from each other and 
\r]\ < 2, no isolated energetic electrons or muons as used in the leptonic ti analyses, 
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Fig. 19: Left: Number of b tags in 1 fb^^ of all-hadronic candidate events selected 
by CDF versus neural network discriminant (A'out)- Right: Number of tags versus 
jet multiplicity after requiring A'out > 0.94. The ti contributions are normalized to 
the measured cross section of 8.3 pb [278]. 



and divided by y/Hr of the selected jets has to be < 3 VGeV. This yields S/B 
of ~l/370. The signal purity is then increased using a neural network discriminant 
based on variables such as Ht, centrality, aplanarity and minimal/maximal invari- 
ant dijet or trijet mass values of all jet permutations. The signal is modeled using 
MC, and the selected data are used directly for the background, as the expected 
signal contribution in these events is very small. 

At least one of the jets in each event is required to be b tagged, and the sample 
composition is then determined in terms of the number of tags rather than events. 
The average number of tags per signal event for a given neural network cut is 
determined from MC, and is used to derive the ti cross section from the observed 
excess in b tags beyond the background expectation obtained from data. The tagging 
efficiencies of the simulation are corrected for differences relative to data. The cut on 
the neural network discriminant (A'out) is optimized for the highest expected signal 
significance after b tagging, taking both statistical and systematical uncertainties 
of signal and background into account, yielding A'out > 0.94. This cut yields S/B of 
«1/12 before b tagging, and 1/2 after tagging. 

The dominant systematical uncertainty in the measurement arises from the un- 
certainty in jet energy scale, strongly impacting both the preselection of events and 
the input variables for the further neural network selection. The sample composi- 
tion of candidate tags in data is illustrated in Fig. [TH and is detailed in Table [9l 
The resulting cross section is given in Table [TOl together with the other all-hadronic 
channel results obtained thus far in Run II. 
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Tabic 9: Expected and observed yields of tags after requiring TVout > 0.94, with 
uncertainties corresponding to statistical and systematic contributions added in 
quadrature. The corrected background (BG) contribution accounts for the signal 
contamination in the dataset used for its estimate. After tagging, 1020 events remain 
in the signal sample with 1233 tags and an expected background of 846 ± 37 tags. 
The average number of tags in ti events is 0.95 ± 0.07 [278]. 





4 jets 


5 jets 


6 jets 


7 jots 


8 jets 


BG 


16060 ± 575 


2750 ± 92 


536 ± 17 


255 ±8 


146 ±5 


BG (corrected) 


15961 ± 677 


2653 ± 112 


481 ± 20 


223 ± 10 


142 ±7 


it {att = 8.3 pb) 


120 ±20 


266 ± 45 


242 ± 41 


101 ± 17 


38 ±7 


BG + tt 


16081 ± 677 


2919 ± 121 


723 ± 46 


324 ± 20 


180 ± 10 


Data 


16555 


3139 


725 


349 


159 



Table 10: ti all-hadronic cross section measurements performed thus far at the 
Run II Tevatron with their integrated luminosities, data selections and analysis 
methods used. All three results have been published. The measurement marked 
with an asterisk refers to rrit of 178 GeV/c^ rather than the standard 175 GeV/c^. 



J£dt 
[fb-i] 


Sel. 


b tag 


(Tjt±(stat.)± (syst . ) ± (lumi. ) 
[pb] 


Ref. 


0.3 


jets only 


yes 


7.5±2.lt-^:^toi M 


[279] 


0.4 


jets only 


yes 


4.51?:° tl;t± 0.3 


[280] 


1.0 


jets only 


yes 


8.3±1.0li:5±0.5 


[278] 



5.1.4. Hadronic r channels 

By choosing a more inclusive ti event selection, hadronic r decays can be included 
as already mentioned in the discussion of lepton + jets and dilepton channels for 
the T + lepton case. A first measurement without any explicit lepton identification 
has been published by CDF, selecting events with at least four jets, with at least 
one being b tagged, and significant not aligned with any jet. Since events with 
isolated energetic electrons or muons are rejected, the resulting sample is enriched 
in T + jets events [281]. 

Explicit reconstruction of hadronic r decays is far more demanding and usually 
relies on multivariate discriminants. Based on the decay mode (1-prong or 3-prong, 
with or without associated electromagnetic subclusters from neutral pions), different 
discriminants can be deployed, exploiting differences between hadronic t decays and 
jets. For example, isolation in the tracking system and calorimetry, shower shape, 
track multiplicity or correlations between tracks and clusters in the calorimeter can 
be applied. Using such discriminants, the r + jets and r + lepton ti decay modes 
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Fig. 20: Hadronic-T Et spectrum before (left) and after (right) b tagging is applied 
in the preselected 1 ft)~^ £ + t sample analyzed by DO. The tt signal is normalized 
to the SM expectation. The highest Et bin contains all overflows [283]. 



have been studied based on their experimental signature of at least one hadronic 
r candidate, and at least four or two jets, and no or one isolated energetic electron 
or muon for the two channels, respectively, b jet identification is crucial to improve 
sample purity in such analyses. 

A first T + jets cross section analysis has been performed by DO on 0.3 fb^^ of 
data, deploying a preselection and r identification as outlined above, b jet identifi- 
cation and neural networks based on event topology and kinematics provide further 
separation of tt signal and background that is mainly due to QCD multijet produc- 
tion, where jets mimic r decays [282]. While the measurement suffers from large 
statistical uncertainties, it is a proof of principle that will benefit greatly from the 
tenfold increased dataset already in hand. 

DO has also performed a first measurement of the tt production rate in the r 
-I- lepton final state based on 1 fb~^ of data [283]. Events are required to have 
exactly one isolated electron or muon with Et > 15 GeV and < 1.1 or Et > 
20 GeV and \ri\ < 2, respectively, at least one r candidate of opposite charge within 
|?7| < 1, between 15 and 200 GeV, at least two jets within \ri\ < 2.5 with 
Et > 20 GeV and the leading jet above 30 GeV, and at least one of these jets 
identified as b jet. Additional channel-specific background rejection is achieved by 
vetoing events collected in the fiT analysis with invariant mass of the isolated muon 
and a second non-isolated muon between 70 and 100 GeV/c^. In the er channel, 
events are rejected where electron and are aligned in azimuthal angle </) by 
requiring cos(A0(e, -^t)) < 0.9. The resulting sample composition is illustrated in 
Fig. [20l and detailed in Table [TT]- there are significant contributions to signal from 
the lepton -I- jets and dilepton channels. 

Background contributions arise from W boson -I- jets production, Z/j* + jets 
events with Z/j* — > i'^i^ /t^t^ , and diboson production, as described by MC. 
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Table 11: Event yields and sample composition before and after b tagging in the 
preselected 1 fb~^ £ + r sample analyzed by DO. The tt signal is normalized to the 
SM expectation. The uncertainties are purely statistical [283]. 





Ix^forc^ 

fIT 


b tagging 

( T 


after 

fiT 


b tagging 

( T 


W^+jets 


38.0 ± 1.7 


34.1 ±3.5 


2.31 ±0.22 


2.13 ±0.27 


Z/^* — > ee, /i/i 


20.7 ± 1.1 


5.8 ±0.6 


1.09±0.11 


0.38 ±0.05 




19.6 ± 1.2 


7.5 ±0.6 


1.02 ±0.10 


0.54 ±0.06 


Diboson 


2.8 ±0.1 


5.1 ±0.6 


0.21 ±0.01 


0.34 ±0.07 


Multijet 


10.6 ±6.3 


12.7±6.6 


4.52 ±3.01 


-1.27± 1.77 


tt^£ + T 


7.8 ±0.1 


6.67 ±0.1 


5.64 ±0.04 


4.70 ±0.05 




4.3 ±0.1 


0.73 ±0.1 


3.14 ±0.03 


0.47 ±0.07 


tt i+ jets 


12.7 ±0.1 


12.41 ±0.2 


8.40 ±0.11 


7.88 ±0.12 


Total Expected 
Data 


116.6 ±6.8 
104 


85.0 ±7.7 
69 


26.33 ±3.02 
29 


15.17 ± 1.97 
18 



Table 12: tt cross section measurements using hadronic r decays performed thus 
far at the Run II Tevatron, with their integrated luminosities, data selections and 
analysis methods used. The first two results have been published; the others are pre- 
liminary. The measurement marked with a double cross (|) includes the luminosity 
uncertainty in the systematic uncertainty. 



jCdt 
[fb-i] 


Sel. 


b tag 


(T( j ± (st at . ) ± (syst . ) ± (lumi . ) 
[Pb] 


Ref. 


0.2 


e + T 


no 


< 5.2 • C7SM (95% C.L.) 


[216] 


0.3 


+ jets 


yes 


6.0±1.2ir7(« 


[281] 


0.3 


T+jets 


yes 


5.1^3:? ±0.7 ±0.3 


[282] 


0.4 


e + T 


no 


significance at « 1 sd 


[284] 


1.0 


e + T 


yes 


8.3t?;°tl;^±0.5 


[283] 


1.2 


e + T 


yes 




[285] 


2.2 


( +T 


\-(\s 


'•■>-1.2-i.i ±<'-' 


[283] 



The W boson ± jets contribution is normalized to data. Background from QCD 
multijet production is estimated from data where the lepton and r are of same 
charge, corrected for significant contributions from W boson + jets and tt pro- 
duction. The dominant systematic uncertainties in this analysis arise from limited 
background/MC statistics, the rate for jets or leptons to mimic r decays, modeling 
of 6 tagging uncertainties, and jet energy calibration. 
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The corresponding cross section is given in Table \T2\ together with an update 
from an additional 1.2 fb~^ of DO data, as well as other results involving hadronic 
T final states obtained thus far in Run II. 

5.1.5. Summary 

An overview of the current status of tt cross section measurements performed in 
different decay channels is given in Fig.[2T]for CDF and DO, showing good agreement 
among channels, analysis methods and experiments. The theoretical predictions 
discussed in Section [3TT] are also shown for comparison as shaded/hatched bands, 
and also show very good agreement with the measurements. 

CDF combines the results obtained in the Icpton + jets and dilepton channels 
using 1.7 — 2.8 fb"-'^ samples of data, achieving a relative uncertainty on the result 
of 9% [286]. The most precise single measurement is obtained in the lepton + jets 
channel using secondary-vertex b tagging on 2.7 fb~^ of data, yielding a relative 
uncertainty of about 10% [274] . DO combines the results from lepton + jets, dilepton 
and T + lepton channels obtained from approximately 1.0 fb~^ of data, yielding a 
relative uncertainty of about 11% [257]. The most precise single DO measurement, 
with a precision of 11%, has been published in the lepton + jets channel, using 
both secondary-vertex b tagging and kinematic information in 0.9 fb~^ of data in a 
combined result [262] . For comparison, the final Run I combined cross section from 
CDF [287,288] and DO [289] had a precision of ^25%. Unfortunately, no combined 
cross section measurement from both experiments exists to date, unlike for top 
quark mass measurements (see Section l7.3.4p . However, such combination for the 
tt production rate is planned in the future. 

With increasing integrated luminosity, the statistical uncertainties are becoming 
less important and precision is starting to be limited by systematic uncertainties, 
which in turn can also be further constrained with more data. One of the main chal- 
lenges for future measurements will be to study systematic uncertainties in greater 
detail rather than use "conservative" estimates. The most precise single tt produc- 
tion cross section measurement at the Run II Tevatron with the anticipated total of 
8 fb~^ of data can be expected in the lepton + jets channel - probably using both 
kinematical and b tagging information - at a relative precision of ^8%. This will be 
dominated by uncertainties on luminosity and systematic effects. For combination 
of results, the precision may ultimately be driven by the current luminosity uncer- 
tainty of 6% for both experiments. The luminosity uncertainty could be avoided 
by measuring cross section ratios, for example, relative to Z boson productior|H. 
Also, with large datasets, a ratio of ti cross sections measured in different chan- 
nels (such as the lepton -I- jets and the dilepton channels) could be obtained with 
good statistical precision while benefiting from cancellations of common systematic 
uncertainties. 



''After completion of this review, first such measurements became available from CDF [292,293]. 
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Fig. 21: Current status of the tt production cross section measurements by CDF 
(top, [290]) and DO (bottom, [291]) compared with SM predictions (shown as 
shaded/hatched bands). The channels contributing to the combined resuhs are high- 
lighted [257,286]. 
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The precision of the measured ti cross section now matches that of the theoretical 
predictions, which provides stringent tests of perturbative QCD calculations. This 
furthers our understanding of the standard model, which still provides an excellent 
description of all current measurements. Based on the observed production rate 
alone, severe constraints on phenomena beyond the SM become feasible [294]. More 
detailed tests of the standard model predictions for ti production will be described 
in the following sections. Via its mass dependence, the it production rate can also 
be used to test consistency of the SM with the top quark mass measurements 
performed at the Tevatron (with the benefit of easier theoretical interpretation 
of the mass parameter as discussed in Section 17.3. 4p . Studying and comparing all 
available ti final states, including those with hadronic r decays, provides a probe 
for novel contributions that can affect the observed final states in different ways. 
For example, searches for charged Higgs bosons in top quark decays are discussed 
in Section 

5.2. Top quark pair production mechanism 

Top quark pair production at the Tevatron proceeds predominantly via qq annihi- 
lation, as described in Section [3. II The remaining fraction from giuon-gluon fusion 
ifgg) is 15 ± 5%, with the uncertainty mainly reflecting that of the corresponding 
PDFs [116]. 

While the total ti production rate has been studied extensively (see Section [5T|) 
and has been found to be in agreement with the SM expectation, the production 
mechanism has not yet been subject to such scrutiny. A measurement of the frac- 
tion of ti events produced via gluon-gluon fusion fgg — (j{gg — s- ti)/a{pp — > ti) 
provides a test of QCD, and can contribute to a reduction in the uncertainties of 
the corresponding PDFs. In addition, contributions from extensions of the SM to 
ti production could be unveiled [295] , some of which may have remained undetected 
because of the presence of new compensating decay mechanisms [296] . 

CDF performs a first measurement of fgg in a 1 fb~^ b tagged lepton + jets 
dataset [297] . The analysis is based on the fact that soft gluons are emitted with a 
higher probability from gluons than from quarks [90-92,298], and the average num- 
ber of charged particles (tracks) with low transverse momentum should therefore 
be higher in gg ti events than in qq ti events. 

To avoid the large theoretical uncertainties on soft-gluon radiation in the Monte 
Carlo modeling of the multiplicity of soft tracks, W + jets and dijet collider data, 
with well understood production mechanisms, are used to relate the observed soft- 
track multiplicity to the gluon content of a sample (see Fig. \22^ . Templates for 
the soft-track multiplicity distribution in gluon-poor and gluon-enriched events are 
obtained, respectively, from W boson events without jets and dijet events with a 
leading-jet Et of 80—100 GeV. The observed soft-track multiplicity distribution 
in ti candidate events is fitted with these templates. From the fit result, fgg is 
extracted and found to be 7 ± 14%(stat.) ± 7%(syst.), corresponding to a 95% C.L. 
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Fig. 22: Left: Correlation between the soft-track multiplicity observed in collider 
data and the average number of gluons in the corresponding Monte Carlo samples. 
Right: Fit result of the soft-track multiplicity distribution observed in 6-tagged lep- 
ton -|->4 jets tt candidate events with a gluon-poor and a gluon-rich template [297]. 



upper limit of 33%. 

CDF uses a complementary second method to extract fgg from the same 
dataset, based on templates from a neural network using kinematic event prop- 
erties to separate gg tt, qq ti and the dominant W + jets background [299], 
yielding a 95% C.L. upper hmit for fgg of 61%. Combining both results yields 
fgg = 7l^7%(stat. -I- syst.), in good agreement with the SM expectation. 

CDF has also performed a first measurement of fgg in 2 fb~^ of dilepton data, 
based on the variation of the azimuthal correlation of the charged leptons for the 
different tt production modes [300]. This difference arises from the fact that, close 
to threshold, top quark pairs are produced in a '^Si state via qq annihilation and 
in a "'^ 5*0 state via gluon-gluon fusion (see Section [53]) . Consequently, the top quark 
spins tend to be antiparallel for tt production via gluon-giuon fusion and aligned 
for production via qq annihilation, which is reflected in the azimuthal correlation 
of the charged leptons. The relative fraction of tt production via gluon fusion is 
determined in a fit of the observed A(f> distribution in data, with templates for 
gg tt, qq tt and background arising from diboson, Z/j* + jets and W boson 
-I- jets production, yielding fgg = 53j^37%(stat.)^g%(syst.), consistent with the SM. 

5.3. Top quark charge asymmetry 

The strong production of top quark pairs is symmetric under charge conjugation at 
leading order, implying it does not discriminate between top and antitop quarks. 
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Considering that the initial proton-antiproton state at the Tevatron is not an eigen- 
state of charge conjugation, this symmetry is a coincidence. At higher orders, a 
charge asymmetry arises from interference between amphtudes that are symmetric 
and antisymmetric under the exchange of top and antitop quarks [301,302], leading 
to an excess of top over antitop quarks in specific kinematic regions. One resultant 
observable is the integrated forward-backward production asymmetry for inclusive 
ti production at the Tevatron. This is predicted to be 5 — 10% at NLO [301-305], 
implying that top quarks are preferentially emitted in the direction of the incom- 
ing protons. The asymmetry depends strongly on the region of phase space being 
probed, and particularly on the production of additional jets: While the asymmetry 
for exclusive ti production without additional jets is predicted to be 6.4% [303], 
the inclusive tt production with one additional jet exhibits an asymmetry of about 
-7% at LO [303,306], which is reduced drastically to (-1.5 ± 1.5)% at NLO [306]. 

The size of higher order corrections for the tt + jet subprocess illustrates that 
higher-order evaluations of the whole process are still necessary for the total asym- 
metry prediction to converge and correctly describe the partial cancellations of the 
various interference contributions. It should also be noted that in the above the- 
oretical predictions, the top quark decay and possible effects on the asymmetry 
from reconstruction of the final-state objects are not considered. The top quark 
charge asymmetry is also sensitive to extensions of the standard model in tt pro- 
duction involving, e.g., axigluons [304], technicolor [307] or additional neutral Z' 
gauge bosons [308] . Consequently, a measurement of the asymmetry can be used to 
set limits on such processes, particularly for extending the sensitivity of searches 
for ti production via heavy resonances (see Section r5. 7. ip to include not only those 
of narrow width but also wide resonances. 

DO has published the first measurement of the integrated forward-backward 
charge asymmetry in ti production, based on a 0.9 fb~^ b tagged lepton -I- jets 
dataset [309]. The ti system is reconstructed using a constrained kinematic fit, 
where the charged lepton is used to differentiate between the top and the an- 
titop quark. The signed rapidity difference of the top and the antitop quark 
Ay = Ut — yi is used as an observable from which the charge asymmetry is obtained 
a,s Afi) — (Nf — Nb)/{Nf + Nb), with Nf {Nb) being the event yields with positive 
(negative) Ay. The sample composition is determined in a template fit based on a 
multivariate kinematic likelihood discriminant for both signs of the reconstructed 
Ay simultaneously as shown in Fig. [23l 

The resulting measurement is not corrected for reconstruction and acceptance 
due to the limited theoretical knowledge of the shape of the asymmetry. Instead, 
a prescription is provided to model the detector acceptance at the particle level, 
for ease of comparison of any model with the measurement. For comparison with 
the standard model, a slightly more detailed prescription than the one provided 
in Ref. [309] is applied to the prediction from the MC@NLO [310,311] generator, 
and found to be in agreement with the measurement for different jet multiplicities, 
including a change in the sign of the asymmetry, as shown in Table 1131 The domi- 
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Fig. 23: Likelihood-discriminant output distribution for data with >4 jets, overlaid 
with the result of a template fit determining the sample composition for events with 
Ay > (left) and Ay < (right) [309]. 

Table 13: ^^'''^: MC@NLO SM prediction for the observed tt charge asymmetry 
in the DO detector, including uncertainties from acceptance and dilution (misre- 
construction of sign in Ay). A"^^: Uncorrected tt charge asymmetry observed by 
DO [309]. 



Njet 


Ap;- [%] 




AT [%] 


> 4 
= 4 

> 5 


0.8±0.2(stat.)±1.0(acc. 
2.3±0.2(stat.)±1.0(acc. 
-4.9±0.4(stat.)±1.0(acc. 


)±0.0(dil.) 
)±0.1(dil.) 
)±0.2(dil.) 


12±8(stat.)±l(syst.) 
19 ± 9 (stat.) ± 2 (syst.) 
-letlr (stat.) ± 3 (syst.) 



nant systematic uncertainty for the >4 jet sample arises from the relative jet energy 
calibration between simulation and data, and for its subsamples from event migra- 
tion between the subsamples when splitting the sample up into one with exactly 
four and one with at least five jets. However, these systematics are negligible com- 
pared to the statistical uncertainties. The measurement is also used to derive 95% 
C.L. limits on the fraction of tt events that are produced via a specific Z' resonance 
model [307,312] with parity- violating couplings as a function of the resonance mass. 

CDF has obtained two measurements of Af^ based on 1.9 fb^^ 6-tagged lepton 
-I- jets data, using different observables after reconstruction of the ti kinematics in 
a constrained fit [313-315]. CDF chooses a different approach for the measurements 
than DO by providing results both before and after background subtraction and 
correction for acceptance and reconstruction effects. 

The first analysis uses as observable the rapidity difference between the hadron- 
ically and semileptonically decaying (anti-) top quark multiplied by the lepton 
charge [313,315]. This is equivalent to Ay in the measurement reported by DO. After 
background subtraction, asymmetries of 0.119±0.064, 0.132±0.075 and 0.079±0.123 
are observed for jet multiplicities > 4, = 4 and > 5, respectively, which is consistent 
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with the measurement by DO, and which correspond to MC@NLO predictions of 
0.017 ± 0.007, 0.038 ± 0.008 and -0.033 ± 0.012 (errors are statistical). The result 
for the inclusive sample with at least four jets is then corrected for reconstruction 
and acceptance effects, yielding A^^" = 24 ± 13 (stat.) ± 4 (syst.)%. The dominant 
systematic uncertainty comes from the shape uncertainty of the Ay-distribution. 
The result is bigger than expected from NLO predictions, but consistent within 
errors. 

The second analysis measures the charge asymmetry using the product of inverse 
lepton charge and cos 9ti^^^ as observable, where Ot^^^^ is the angle of the top quark 
with the hadronic decay chain relative to the proton beam [314,315]. Since this 
measurement is performed in the laboratory frame rather than the parton rest 
frame, the asymmetry is reduced by about 30% [304]. For > 4 jets, the corrected 
asymmetry is ^f^"^ = 17 ± 7 (stat.) ± 4(syst.)%, with the dominant systematic 
uncertainty arising from background shape and its normalization. This result is 
consistent with the theoretical prediction at the level of two standard deviations 
(sd) for a Gaussian distribution. 

It should be noted that the forward-backward asymmetry in the laboratory 
frame vanishes at the LHC due to the symmetric initial state, in contrast to the 
Tevatron pp collider. The dominance of it production via the charge symmetric 
gluon fusion reaction at = 14 TeV reduces the observable charge asymmetry at 
the LHC. 

5.4. Top quark pair production kinematics 

New contributions to tt production can alter the observed event kinematics, which 
can be exploited in searches for such processes, as described in Section [?771 The basic 
kinematic properties of leptons, jets, and corresponding angular distributions 
are continuously compared to the SM expectation, both in signal-enriched datasets 
and signal-depleted control samples exhibiting features similar to the signal, in all 
studies of top quark properties and especially in the cross section analyses. Thus 
far, no significant deviation from the SM expectation has been found that would be 
indicative of new physics contributions to top quark samples. 

In Run I, CDF and DO observed a slight excess of the tt production rate over 
the SM prediction in the dilepton channel, especially in the efj, final state [316,317]. 
Since some of these events had rather large '^t and lepton-pT, their consistency 
with the SM was questioned and, for example, the kinematic compatibility of these 
events with cascade decays of heavy supersymmetric quarks was pointed out [318]. 
Triggered by this, CDF performed a search for anomalous ti kinematics in Run II, 
based on 0.2 fb^^ dilepton data yielding 13 candidate events [319]. 

A priori four kinematic event variables, including |(t and leading-leptonpT, were 
chosen to quantify any possible deviations of the observed distributions from SM 
predictions. Using a shape comparison based on the Kolmogorov-Smirnov statistic, 
no significant discrepancy was found, and the probability of observing a sample less 
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consistent with the SM was determined to be 1.6%. Including systematic uncertain- 
ties, the p-value was 1.0-4.5%, where the lowest value was obtained by lowering the 
background expectation by one standard deviation. Consequently, presence of pro- 
cesses beyond the SM, with high and high lepton-p^, are not favored by these 
data. 

It is also of great interest to study the kinematic properties of the top quark 
itself and to compare these with the SM expectation. In Run I, DO performed such 
an analysis and found good agreement with the SM [320,321], which was then 
also confirmed through a dedicated study of the px spectrum of the top quark 
by CDF [322]. A corresponding analysis has not yet been published in Run II. 
However, a measurement of the differential tt production cross section da/dMtt has 
been performed by CDF using 1.9 fb~^ of Run II data, as described in Section [5. 7. 3[ 
and also shows good agreement with the SM [323]. 

5.5. Spin correlations in tt production 

Top quark pairs are expected to be produced essentially unpolarized in hadron 
collisions when the incident particles are unpolarized. A small polarization at the 
percent level is induced by QCD processes [324-326], and is perpendicular to the 
production plane, as strong interactions conserve parity. A measurement of this 
effect will be very difficult both at the Tevatron and the LHC, which in return 
suggests to use a corresponding analysis to probe for non-standard contributions 
in it production [327]. An even smaller additional polarization in the production 
plane arises from mixed strong and weak contributions to tt production at order 
a^a [328]. 

While no observable spin polarization in ti production is predicted in the frame- 
work of the standard model, the spins of the top and the antitop quark are expected 
to be correlated [329] . This correlation depends both on the production mode of the 
it pair and the production energy. Close to threshold, the top quark pair is produced 
in a ■^S'l state via qq annihilation and in a ^Sq state via gluon-gluon fusion [330,331]. 
Consequently, the top quark spins are (anti-) parallel and the top quarks have op- 
posite (same) helicities for it production via qq annihilation (gluon-gluon fusion). 
Above threshold, this simple picture becomes more complicated, as effects of orbital 
angular momentum must be taken into account. In the high-energy limit, where the 
top quark mass can be neglected, the conservation of chirality dictates that top 
and antitop quarks be produced with opposite helicities. Since it production at the 
Tevatron proceeds mainly via qq annihilation, as opposed to the LHC where the 
main contribution comes from gluon-gluon fusion, the observable correlation will 
have opposite signs at the two colliders [332]. 

The short lifetime of the top quark (see Section I3.3.2P assures that its spin in- 
formation is passed on to its decay products, and is reflected in their corresponding 
angular distributions. This provides experimental access to spin correlations, and 
a way to check whether the top quark can indeed be considered as a free quark. 
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The resulting indirect limits on the top quark lifetime (see Section [7?^ can provide 
limits on the CKM matrix element \Vtb\, free from the assumption of three quark 
families [333] together with the measurement of the top quark branching fractions 
(Section 16. 2p . In addition, spin correlations probe the dynamics of top quark pro- 
duction and decay for possible contributions from physics beyond the SM. 

The down- type (T3 = — ^) decay products of the W boson from top quark decay 
are most sensitive to the original top quark spin. Their angular distribution in the 
top-quark rest frame is described by 1 -f cos 6*, with 6 being the angle between the 
line of flight of the down-type fermion and the top spin direction. The experimen- 
tal difficulties of distinguishing between jets from up-type and down-type quarks 
(charm tagging would help only in 50% of the cases) can be avoided by focusing on 
the dilepton final state, where the charged (down-type) leptons are clearly identified. 

At the Tevatron, an optimal spin-quantization basis is provided by the "off- 
diagonal" basis [334,335], where the spins of top and antitop quarks produced 
by qq annihilation are fully aligned for all energies, and only the contribution of 
top quark pairs from gluon-gluon fusion leads to a reduction of the correlation. 
The off-diagonal basis is defined via the top quark's velocity /3* and scattering 
angle 6* in the collision rest frame of the incoming partons. The quantization axis 
then forms an angle ip with the proton-antiproton beam axis, defined by tan^ = 
/3*2 sinO* cos6'*/(l - /3*^ sin^ 9*). Consequently, in the limit of /3* ^ (top quark 
production at rest), the spins of top and antitop quarks point along the beam axis 
in the same direction. At very high energies, the spins are aligned with respect to 
the direction of the ti momenta. 

Using as observables the angles 0+ and 0_ of the down-type fermions relative to 
the quantization axis in the rest frame of their respective parent (anti-) top quark, 
the spin correlation is given by [336] : 



where the correlation coefficient k is predicted to be -f 0.88 at Run I of the Tevatron 
when using the off-diagonal basis. Since the distribution is symmetric under ex- 
change of the two angles, an electric charge measurement of the top decay products 
is not necessary. 

DO performed a first search for evidence of spin correlations in tt production in 
Run I, using a 0.1 fb^^ dilepton data sample containing six candidate events [337]. 
From the dependence of a likelihood function on k, at 68% C.L. a lower limit on k 
of —0.25 was extracted. This is in agreement with the standard model expectation 
and disfavors anti-correlation of spins (k = — 1) that would arise from ti production 
via a scalar particle. While this limit is rather weak, it is a proof of principle that 
the analysis can be performed. Unfortunately, there has been no result as yet from 
Run II, although it would greatly benefit from the increase in data already available. 
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Fig. 24: Left: Ht distribution observed in DO's 2.1 fb~^ lepton + jets data with four 
jets and at least three h tags, compared with expected standard model background 
processes and tiH signal scaled up by a factor of 100. Right: Event display {xy- 
view along the proton beam direction) for the triple h tagged event of highest Ht 
(444 GeV). The first three jets have h tags [338]. 



5.6. Search for associated Higgs boson production 

DO performs a first search for associated tt and standard model Higgs boson pro- 
duction in the ttbb final state in 2.1 fb~^ of b tagged lepton + jet data [338]. While 
the observation of a significant signal in this channel is beyond the reach of the 
Tevatron, this analysis can nevertheless contribute to future combinations of the 
Tevatron searches for Higgs bosons at low masses, as favored by the standard model 
(see Section I7.3.4|) . The investigated events exhibit high jet and b tag multiplicities 
that were not studied separately before. It is therefore interesting to search for de- 
viations from the SM predictions that could arise, e.g., from anomalous top-Higgs 
couplings [339] or from a new quark singlet of charge |e [340]. 

The signal signature has the tt lepton + jets event characteristics, but with two 
additional b jets from the Higgs boson decay. The main background arises from 
ti with additional heavy-flavor jet production, but also W boson -|- jets and QCD 
multijet production contribute to the background. For signal discrimination, the 
shape of the Ht distribution of the selected jets is used in events with four or 
at least five jets and one, two or at least three b tags. The observed events in all 
these distributions are consistent with background expectation, which is especially 
interesting for the events with at least three b tags that were studied separately for 
the first time. Figure [24] shows the observed Ht distribution for events with four 
jets and at least three b tags, and an event display for the triple b tagged event of 
highest Ht- 

Since no signal is observed, 95% C.L. limits on ttH production multiplied by 
B{H — > bb) are derived for Higgs boson masses between 105 and 155 GeV/c^. For 
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a Higgs boson mass of 115 GcV/c^, the expected limit is a factor of 45 larger 
than the SM production rate, while the observed limit is a factor of 64 above the 
SM expectation. Optimization of the preselection (currently corresponding to the 
standard it selection) for the flH signal, and of the signal discrimination using 
additional kinematic variables, is underway. 

5.7. Search for top quark pair production beyond the standard 
model 

5.7.1. Search for a narrow-width resonance decaying into ti 

The cixistence of yet undiscovered heavy resonances could be revealed through their 
decays into top quark pairs, which would add a resonant contribution to the stan- 
dard model process. Theories beyond the standard model predict the existence 
of massive Z-like bosons, for example, such as Kaluza-Klein excitations of the 
gluon [341] or of 7 and Z bosons [342], extended gauge theories [343,344], mas- 
sive axigluons with axial vector couplings [345] or topcolor [346,347]. 

The wealth of such models demonstrates the importance of model-independent 
searches. One general way an additional production mode can be observed - pro- 
vided the resonance X decaying to ti is sufficiently heavy and narrow - is to an- 
alyze the ti invariant mass distribution for an excess over expectation. In the cor- 
responding analyses performed at the Tevatron, no significant deviations from the 
SM expectation have thus far been observed, resulting in 95% C.L. upper limits on 
ax ■ B{X ^ tt) as a function of resonance mass Mx- 

These results can be used to set lower mass limits for specific benchmark models 
that provide easy comparison. For example, topcolor [346, 347] provides a dynamic 
electroweak symmetry breaking mechanism via a top quark pair condensate [348] Z' , 
formed by a new strong gauge force, that couples preferentially to the third fermion 
generation. Particularly, a topcolor-assisted technicolor model [307, 312] predicts 
this Z' boson to couple strongly only to the first and third generation of quarks, 
while exhibiting no significant coupling to leptons. This leptophobic and topophyllic 
Z' boson has a significant cross section a{pp Z' ^ tt) that is observable at the 
Tevatron for a variety of masses and widths, and is used as a reference model. 

CDF and DO performed modcl-indcpcndcnt searches for narrow massive vector 
bosons decaying into ti already in Run I in lepton -|- jets datasets of 106 pb~^ and 
130 pb~^, respectively. Using the best kinematic fit to the ti hypothesis in each 
event, the tl invariant mass distribution was reconstructed and no excess observed 
above expectation. The resulting upper limits on ax ■B{X ti) arc turned into 95% 
C.L. mass limits of Mz' > 480 GeV/c^ for CDF [349] and Mz' > 560 GeV/c^ for 
DO [350]. For these results, a width of the Z' (or X) of 1.2% of its mass is assumed, 
which is well below the detector mass resolutions for ti systems. Consequently, the 
results are dominated by detector resolution and independent ofTz' for values below 
the mass resolution of a few percent (« 0.04M^/ for DO in Run I [351]). This kind 
of resonance width is also used in the Run II measurements described below. 
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Fig. 25: Left: Expected and observed tt invariant mass distribution in lepton + jets 
events with four or more jets. Right: Expected and observed 95% C.L. upper limits 
on ax ■ B{X tt) [352]. 



In Run II, both CDF and DO search for a generic heavy resonance X of narrow 
width (Tx = 0.012Mx) compared to the detector mass resolution in h tagged lepton 
+ jets datasets. The ti invariant mass spectrum is reconstructed using either the 
best kinematic fit to the tt production hypothesis (CDF) or directly from the four- 
momenta of the up to four leading jets, the lepton and the neutrino momentum (DO). 
The latter approach was shown to provide greater sensitivity for large resonance 
masses than the previously used constrained kinematic fit, and also allows for the 
inclusion of data with fewer than four jets in case jets merged. As both experiments 
observe no significant deviation from SM expectation, 95% C.L. upper limits on 
<7x ■ B{X -> tt) are given for values of Mx between 450 and 900 GeV/c^ (CDF) and 
Mx between 350 and 1000 GeV/c^ (DO, see Fig. US]) in increments of 50 GeV/c^. 

Both experiments provide 95% C.L. mass limits for the leptophobic topcolor- 
assisted technicolor Z' boson as a benchmark model. With 1 fb^-'^, CDF finds 
Mz' > 720 GeV/c^ (expected limit = 710 GeV/c^) [353], while DO finds Mz' > 760 
GeV/c^ (expected limit = 795 GeV/c^) [352] using 2.1 fb"^ of data, which super- 
sedes a previous result on 0.9 fb~^ of data [354]. CDF also obtains a result on a 
subset of 0.7 fb"^ of the data analyzed above, using an untagged lepton -I- jets sam- 
ple where b tag information only contributes as a way to reduce jet combinatorics 
in a standard model tt matrix element based reconstruction of M^f. This yields a 
slightly better hmit than the analysis on the fuU 1 fb"^ of Mz' > 725 GeV/c^ [355]. 

For future studies, it would be interesting to see how sensitive the observed 
limits are to the assumption of Z boson-like couplings used in the analyses. The 
limits obtained apply to resonances of narrow width only. Wider resonances could 
be detected by studying the tt differential cross section (see Section I5.7.3P or the 
forward-backward charge asymmetry in tt production (see Section [5. 3p . 
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Fig. 26: Left: Expected and observed reconstructed tt invariant mass distribution 
in b tagged lepton + jets events with at least four jets. Right: Corresponding ti dif- 
ferential cross section after background subtraction and unfolding, compared to the 
SM expectation [323]. 

5.7.2. Search for tt production via a massive gluon 

Instead of a new color singlet particle decaying into ti , as described in the previous 
subsection, there could also be a new massive color octet particle G contributing to 
tt production [347]. Such a "massive gluon" production mode would interfere with 
the corresponding standard model production process. Assuming SM top decay, 
CDF has performed a search for a corresponding contribution by comparing the 
ti invariant mass distribution in a 1.9 fb~^ b tagged lepton + jets dataset with the 
standard model expectation. As the largest discrepancy relative to the standard 
model is found to be 1.7 sd for masses and widths of 400 GeV/c^ < Mg < 800 
GeV/c^ and 0.05 < Tq/Mq < 0.5, respectively, 95% C.L. upper and lower limits 
are extracted on the corresponding coupling strengths of such massive gluons [356] . 

5.7.3. Measurement of the tt differential cross section {da / dM^i) 

Since new production mechanisms for top quark pairs could manifest themselves 
in the ti invariant mass distribution as resonances of different widths or, more 
generally, as shape distortions [357], one approach for detecting such contributions 
is to compare the shape of the observed differential ti production cross section 
da/dMft with the SM expectation. 

CDF reconstructs the ti invariant mass spectrum in a 1.9 fb~^ b tagged lep- 
ton -|- jets data sample (see Fig. [26]) by combining the four-vectors of the four 
leading jets, lepton and ^t- After subtracting the background processes, the dis- 
tortions in the reconstructed distribution due to detector effects, object resolutions 
and geometric/kinematic acceptance are corrected for through a regularized unfold- 
ing technique [358] . From the unfolded distribution, the ti differential cross section 
da /dMii is extracted and its shape compared with the SM expectation. The compar- 
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ison shows good agreement with the standard model, yielding an Anderson-Darling 
p-value of 0.45 [323]. 

5.7 A. Search for new heavy top-like quark pair production 

The number of light neutrino species (for m,y < mz/2) has been determined to 
be Ni, — 2.9840 ± 0.0082 based on the invisible Z boson decay width in precision 
electroweak measurements [68] . This rules out a fourth generation of fermions with 
a light neutrino v^. However, the existence of a fourth generation is consistent with 
precision electroweak data for a fermion mass range mz/i < to/4 < 0(1 74 GeV /c ) , 
even without introducing new physics [359,360]. Fourth-generation quark masses up 
to 400 GeV/c^ are compatible with current measurements, and are constrained to 
exhibit small mass splitting, so that decays of an up-type fourth-generation quark 
into Wq {q — d, s, b) are preferred [361]. Such an additional generation would have 
a drastic impact on the phenomenology of the Higgs boson, thereby relaxing the 
mass bounds obtained from the SM up to 750 GeV/c^ at 95% C.L., and altering 
expected kinematics and production rates. 

The existence of a fourth chiral-fermion generation is predicted by various exten- 
sions of the standard model, for example, in an S0(l,13) framework unifying charges 
and spins [362], or in models with flavor democracy [363,364]. Other models that 
add more exotic additional heavy quarks that can decay via Wq have been brought 
up as well [365]. For example, the "beautiful mirrors" model [366] introduces mir- 
ror quark-doublets with the same quantum numbers as their SM counterparts, but 
with vector couplings to the W boson. This addition helps to improve the fit of elec- 
troweak observables by removing the observed discrepancy in the forward-backward 
asymmetry of the b quark. 

CDF performs a search for pair production of heavy top-like quarks {t'P) that 
do not necessarily exhibit SM-like up-type fourth-generation properties in terms of 
charge or spin. The analysis is based on the assumptions that the t' is pair-produced 
via the strong interaction, has a mass larger than that of the top quark, and decays 
promptly into a W boson and a down-type d, s, b quark with 100% branching ratio. 
As a consequence, the t' decay chain is identical to that of the top quark, and t't' 
production can be sought in a lepton -I- jets sample selected solely based on event 
kinematics to not restrict the search to Wb final states by using b tagging. 

The t' signal can be distinguished from SM background using, e.g., the observed 
distributions of total transverse energy Ht in the event based on lepton, jets and 
^T, and the reconstructed t' mass (Afrcco) from the best kinematic fit to the t't' 
hypothesis in each event (see Fig. [27]) ■ Superseding a previously published result 
based on 0.8 fb^^ [368], CDF uses a two-dimensional binned likelihood fit in Ht 
and Afioco to separate SM background and t' signal in 2.8 fb~^ of data [367]. 

Since no evidence for t' production is found, 95% C.L. upper limits on the t't' pro- 
duction cross section (assuming B{t' Wq) — 100%) are derived for 180 GeV/c^ 
< nit' < 500 GeV/c^. Assuming SM couplings, a 95% C.L. lower limit on the fourth- 
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Fig. 27: Left: Ht versus Mjoco distribution observed in data (black points) overlaid 
with the fitted number of standard model background events from ti , T^+jets and 
QCD. Right: Expected and observed 95% C.L. upper limits on the t'i' production 
cross section, assuming 100% B{t' Wq). The dark/light regions represent the 1/2 
sd bands on the expected limit [367]. 

generation t' mass of 311 GeV/c^, based on the calculations in Refs. [116,117], is 
obtained, where the largest systematic uncertainty arises from the jet energy scale. 
The deviation of the observed limit from its expected value above «400 GeV/c^ 
is being investigated. Using a priori defined groups of bins in Ht and Mreco, the 
p- values to observe at least the number of events found in data, given the SM ex- 
pectation, are evaluated. With the smallest p-value being 0.01, the excess in the 
data tails is concluded to be not statistically significant. 

5.7.5. Search for scalar top quark production 

Many processes beyond the standard model exhibit signatures similar to tt events. 
Consequently, tt data samples can, in principle, contain admixtures of such contri- 
butions. 

For example, the Minimal Supersymmetric Standard Model (MSSM) [369] pre- 
dicts that supersymmetric partners of the top quarks, scalar top or "stop" quarks, 
are predominantly produced in pairs via the strong interaction just like SM top 
quarks. The stop-quark pair-production cross section has been calculated at NLO 
supersymmetric-QCD and depends mainly on the stop quark mass and very little 
on other supersymmetric parameters [370]. At a center of mass energy of 1.96 TeV, 
the pair production cross section for the hghtest stop quarks {iiii) of 175 GeV/c^ 
mass is 0.58 pb [371], or roughly 10% of the SM ti production rate (see Section [XT]) . 
The observable final states from stop decays depend strongly on supersymmetric 
parameters, especially the masses of supersymmetric particles in the decay chain. 
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Fig. 28: Left: Multivariate kinematic discriminant distribution for electron +>4 jets 
events in 0.9 fb^^ of DO data and simulated signal and background. The dashed 
black hne shows the expected signal shape for a i\ {Xi) mass of 175 (135) GeV/c^, 
enhanced by a factor of ten. Right: Expected and observed 95% C.L. limits on the 
tit\ production cross section, together with the prediction for different stop and 
chargino mass combinations [376]. 

In Run II, the decay mode ti — > M^v^ has been studied by DO [372] in 0.4 fb^^ of 
data. The decay channel ti cx°, where the lightest neutralino Xi is the lightest 
supersymmetric particle (LSP), was studied by CDF [373] in 0.3 fb~^ and by DO 
in 0.4 and 1 fb^^ of data [374,375]. 95% C.L. mass-exclusion limits on the involved 
supersymmetric particles were provided for both decay channels. 

Another important decay mode for stop is ii Xi^^ where the lightest chargino 
Xi decays to W^x\, resulting in final states identical to those from ti decays, but 
with two additional neutralinos (LSPs) that contribute to ^t- First limits for this 
channel were provided by CDF in Run I on 0.1 fb^^ of lepton + jets data [377]. 

DO performs the first study of this channel in Run II in 0.9 fb~^ b tagged lepton 
+ jets data [3760, assuming a neutralino mass of 50 GeV/c^ (slightly above the 
limit set by LEP [379]) while varying the stop and chargino masses between 145 
and 175 GeV/c^, and 105 and 135 GeV/c^, respectively. 

A possible stop admixture in the sample is searched for by employing a mul- 
tivariate discriminant based on kinematic event properties (see Fig. [28|) . with the 
main challenge being to separate the topologically similar ti background from titi 
signal. Counterintuitively, the additional neutralinos do not provide large differences 
in that can be exploited, but the larger chargino mass compared to that of the 
W boson reduces the phase space for h jets in the event. 

Since no significant signal admixture in the lepton -I- jets dataset is observed, 
95% C.L. upper limits on the tiii production rate are set that are a factor of « 7—12 
above the prediction, as illustrated in Fig. [251 Consequently, the stop quark masses 



■^An updated version of the result has been published after completion of this review, see Ref. [378]. 
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Observed 95% C.L. Observed 95% C.L. 




M(ti) [GeV/cT M(t|) [GeV/c"] 

Fig. 29: Excluded areas observed at 95% C.L. in the neutralino versus stop mass 
plane for various assumed dilepton branching fractions and two different chargino 
masses (left: 105.8 GeV/c^, right: 125.8 GeV/c^), obtained by CDF in 2.7 fb^^ 
of dilepton data. The contributions of e,/^,r to the final state are assumed to be 
equal [380]. 



considered cannot yet be excluded, and this analysis should greatly benefit from 
the increased datasets already in hand. The weaker observed limits relative to their 
expected values are driven by the muon + jets channel. The corresponding excess in 
data was tested with pseudo-datasets to be statistically consistent with the standard 
model expectation. 

CDF searches for a stop admixture in the ti dilepton channel using a 2.7 fb^^ 
dilepton dataset of both b tagged and untagged events [380] . Assuming Xi to be the 
LSP, heavy sfermions, the stop mass below the top mass, and the chargino mass 
smaller than the mass difference between stop and b quark, the decay ii Xi^ ob- 
tains 100% branching fraction. The dilepton final state resulting from xt ~^ Xi^^'^ 
decays is then identical to the ti final state, but with two additional neutralinos 
contributing to -^t- It can be reached through a variety of chargino decay channels, 
resulting in variations of the branching fraction depending on SUSY parameters. 

The stop quark signal is discriminated from standard model background using a 
single quantity, the reconstructed stop mass, in a fit to the observed data distribu- 
tion. The mass is reconstructed from this underconstrained system by treating the 
neutralino and neutrino from each stop decay as one massive pseudo-particle, and 
then applying a standard top mass reconstruction technique in the dilepton channel, 
neutrino weighting (see Section 17.3. 1|) . Since the observed distributions are consis- 
tent with standard model processes, 95% C.L. limits are extracted on the dilepton 
branching ratio in titi production for stop masses 115 — 185 GeV/c^, neutralino 
masses 43.9 — 88.5 GeV/c^, and chargino masses 105.8 — 125.8 GeV/c^, as illus- 
trated in Fig.[29l A 100% branching fraction for ii xt^^ ^^'^ equal contributions 
from e,/i and r to the final state are assumed throughout. 

Model-independent searches for novel admixtures in top quark samples, via the 
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search for anomalous event kinematics, were discussed in Section [STU 
5.8. Single top quark production 

Electroweak production of top quarks without their antiparticles can provide a di- 
rect measure of the |Vtfc| CKM matrix element, test the Wth vertex structure, and 
probe for physics beyond the standard model, such as flavor changing neutral cur- 
rents or new heavy gauge bosons W' (see Section [3T2t . A thorough understanding 
of single top quark production is also important for studies of processes with sim- 
ilar signatures, such as standard model VF-Higgs associated production, for which 
this process constitutes a major background. While the single-top production rate 
is «40% of the strong ti production, the signal extraction from background is very 
challenging because only one top quark signature is present in the final state. Simple 
kinematic selections are insufficient for such an analysis, and sophisticated multi- 
variate techniques have to be deployed. 

For single top quark production at the Tevatron, only s- and t-channel produc- 
tion are relevant, contributing, respectively, O.SS^q'^i pb and l-98_g;22 pb to the 
total rate at NLO [156]. The experimental signature comprises a b jet and the W bo- 
son decay products from the top quark decay. In the s- {tb-) channel, one additional 
b jet arises from the b quark produced together with the top quark. In the t- {tqb-) 
channel, a forward light-quark jet accompanies the production of the top quark, 
sometimes along with another b jet from the gluon splitting into bb (see Fig. [8]). In 
order to suppress multijet background, the W boson is usually required to decay 
leptonically into an electron or muon and corresponding neutrino. Consequently, the 
final state signature of single top quark production contains an energetic isolated 
electron or muon, and two or three jets, with at least one of them being a b jet. 
As usual, additional jets can arise from initial- or final-state radiation. 

5.8.1. Production cross section and Vtb 

Searches for single top quark production were already performed in Run I using 0.1 
fb-i of data by DO [381,382] and by CDF [383,384], yielding upper limits on the 
production rate that were at least a factor of six larger than the SM expectation. 
In Run II, first results were pubhshed using 0.2 fb^^ of data by CDF [385] and 
DO [386,387], where the best observed limit was less than a factor of three greater 
than the SM prediction. Finally, DO published first evidence for single top quark pro- 
duction using 0.9 fb^^ of data [21,22], observing a signal of 3.6 standard deviations 
significance. Preliminary results from CDF based on 2.2 fb~^ confirmed evidence 
for single top quark production [388-391] with an observed signal significance of 3.7 
standard deviations obtained by combining three of these analyses [392jl. 




The analyses yielding first evidence for electroweak top quark production apply 
event selections requiring one energetic isolated electron or muon and -^t- CDF uses 




•^A slightly updated version of the result was published by CDF [23] after completion of this review. 
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Table 14: Expected and observed event yields of the single top selections for e and 
fi, single and double b tagged channels combined (left for DO based on 0.9 fb~^ [21], 
and right for CDF based on 2.2 fb-^ [390]). For the DO result, the overall W + jets 
background includes Z + jets and diboson events. 



Source 


2 jets 3 jets 


4 jets 


Source 


2 jets 


3 jets 


tb 


16±3 8±2 


2±1 


tb 


41.2±5.9 


13.5±1.9 


tqb 


20±4 12±3 


4±1 


tqb 


62.1±9.1 


18.3±2.7 


tt 


59±10 135±26 154±33 


tt 


146.0±20.9 


338.7±48.2 


Wbb 


261±55 120±24 


35±7 


Wbb 


461.6±139.7 


141.1±42.6 


Wcc,Wcj 151±31 85±17 


23±5 


Wcc,Wcj 


395.0±121.8 


108.8±33.5 


Wjj 


119±25 43±9 


12±2 


Wjj 


339.8±56.1 


101.8±16.9 


Multijets 


95±19 77±15 


29±6 


Multijets 


59.5±23.8 


21.3±8.5 


BG Sum 


686±41 460±39 253±38 


Dibosons 


63.2±6.3 


21.5±2.2 


Data 


697 455 


246 


Z+jcts 


26.7±3.9 


11. Oil. 6 








BG Sum 


1491.8±268.6 


754.8±91.3 








Data 


1535 


712 



events with two or three jets and one or two b tags, while DO includes additionally 
events with four jets, where the extra jet arises from initial- or final-state radia- 
tion. The signal acceptances for the s- and t-channel are 2.8% and 1.8% (CDF), 
and 3.2% and 2.1% (DO). The expected and observed event yields are shown in 
Table [T4l The dominant background contributions come from M^-|-jets production, 
it production in the lepton -I- jets or dilepton final states, where one jet or lepton 
is not reconstructed, and from multijet production. The main sources of systematic 
uncertainty are background normalization, jet energy scale, and the modeling of the 
b tagging and triggers. As can be appreciated from the table, the uncertainty on 
the background is larger than the expected signal, which makes advanced analysis 
techniques necessary. 

DO applies three different multivariate analysis techniques to the preselected data 
sample: boosted decision trees (BDT) , Bayesian neural networks (BNN) and matrix 
elements (ME), where the latter two reflect reoptimized studies [22] of previous 
work [21]. Being based on leading-order matrix elements for the description of signal 
and background processes, the ME analysis does not use four-jet events. For each 
analysis, the combined s- and t-channel cross sections are extracted from the peak 
of the Bayesian posterior probability density derived from a binned likelihood of 
the respective discriminants. The results are then combined, yielding: 

(johsf^pp _,tb + X, tqb + X) = A.^tW pb (BDT, 3.4 sd) 

- AAt\l pb (BNN, 3.1 sd) 
= pb (ME, 3.2 sd) ^ ' 

= A.lt\l pb (Combined, 3.6 sd), 

where the uncertainties correspond to the combination of statistical and systematic 
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sources. The observed production rates are in agreement with SM expectation and 
with each other. The significances in the parentheses are obtained from studies of 
large ensembles of pseudo-experiments. The expected sensitivity of the combined 
result is 2.3 standard deviations, indicating that the measurement benefited from 
a statistical upward fluctuation. Separate measurements of the s- and t-channel 
cross sections are also performed with the BDT analysis. The results are Ug = 
1.0 ± 0.9 pb and at = 4.2^14 pb, where the other channel (not measured) is set to 
its SM expectation (cr,, = O.SSt^J;^! pb and at = 1.98l^;22 Pb [156]). The observed 
enhancement in the i-channel with respect to the standard model prediction is not 
statistically significant. 

CDF uses the following multivariate analysis techniques on their preselected 
dataset: neural networks (NN [388]), a likelihood function (LHF [389]), a matrix 
element discriminant (ME [390]) and boosted decision trees (BDT [391]). The results 
are: 

^obs(-p^ ^ t6 + X, tqh + X) = pb (NN, 3.2 sd) 

= 1.8t|^:^ pb (LHF, 2.0 sd) 
= 2.2t°:? pb (ME, 3.4 sd) (28) 

= 2.21:°:^ pb (Combined, 3.7 sd) 
= 1.9t^:? pb (BDT, 2.8 sd), 

where the uncertainties given are both statistical and systematic. The BDT anal- 
ysis became available after the combination of results [392] , and was therefore not 
included in that compilation. The observed results agree with each other and with 
the standard model. The expected sensitivity of the combination is 5.1 standard 
deviations, pointing to a statistical downward fluctuation in the data. 

Figure [3D] shows the discriminant outputs of the two most signiflcant single 
measurements from CDF (ME) and DO (BDT). A graphical summary of the mea- 
surements and a comparison with the standard model expectation is given in Fig. 
[3T] All analyses assume a top quark mass of 175 GeV/c^, B{t Wh) — 100%, and 
the SM ratio for s- to ^-channel cross sections. 

In a recent update, CDF has added 0.5 fb""'^ of data to its single top sample, and 
thereby increased the observed significance for all previous analysis techniques. The 
matrix element analysis yields again the most significant single result, exceeding 
four standard deviations. A combination of the measurements has not yet become 
available. The results for 2.7 fb^^, with combined statistical and systematic uncer- 
tainties, are [393-396]: 

^obs(-pp ^tb + X, tqh + X) = 2.lt°-^ pb (NN, 3.7 sd) 

= 2.0t°;^ pb (LHF, 2.6 sd) 
= 2.7t°:? pb (ME, 4.2 sd) ^ ' 

= 2.4t°;^ pb (BDT, 3.6 sd). 

Since the single top quark production rate is proportional to [Vtb]^, the ob- 
served cross sections can be turned into measurements of \Vtb\ under the following 
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Fig. 30: Multivariate discriminant outputs observed in single-top candidate events 
compared to contributions from signal and background processes. Left: CDF ma- 
trix element discriminant, with yields normalized to their SM predictions [390]. 
Right: DO BDT output in the single-top signal region, with signal normalized to the 
measured cross section [21]. 



assumptions: (i) there are no single top quark production modes beyond the SM, 

(ii) single top quark production and decay are dominated by the Wtb interaction 
(jVtbl ^ \Vtd\, \ Vts\, as indicated by measurements of R described in Section l^?^ . and 

(iii) the Wtb interaction exhibits aV-- A structure and is CP conserving. The latter 
premise allows for anomalous left-handed vector couplings f[ (see Section 13.3. 3p . 
but not for right-handed vector or tensor couplings. Anomalous values {j^ 1) 
do not affect the it production rate or kinematics, nor tb or tqb kinematics, but 
simply rescale the single-top production rate. Consequently, \Vtbfi \ extracted from 
single-top production can be > 1, and constraining the measurement to lie between 
0—1 implies that /f = 1, as predicted by the standard model. The measurements 
of |Vtb/f I and |Vtb| are independent from the number of fermion generations and 
unitarity of the CKM matrix. 

Using the result of the BDT analysis and a positive flat prior for |Vtfc|^, DO ob- 
tains \Vtbfi\ = l-SltoJi- Restricting the prior to [0, 1] yields \Vtb\ = l-OOl^;?^, with 
a corresponding 95% C.L. lower limit of \Vtb\ > 0.68 [22]. CDF uses its combined 
measurement in the same way to obtain a 95% C.L. lower limit of \Vtb\ > 0.66 [392]. 
The matrix element analysis based on 2.7 fb^^ yields \Vtb\ > 0.71 [395]. 

With most of the current measurements giving evidence of > 3 sd for single top 
quark production, the observation at the five standard deviation level seems immi- 
nent at the Run II Tevatron. Extrapolating from the 2.2 fb~^ result, as illustrated 
in Fig. [32l CDF estimates that a > 5 sd significance should be reached by adding 
one more fb^^ of data to the analyses. DO could reach that level in the 2.3 fb~^ 
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Fig. 31: Cross section measurements of first evidence for single top production, and 
the combined results from CDF and DO compared to the SM prediction. 



dataset that is currently being analyzeqj. 

As discussed in Section 13.21 s- and t-channel production exhibit different sen- 
sitivity to physics beyond the standard model. Measuring both rates separately 
provides a valuable tool to check for various exotic model contributions to single 
top quark production. Figure [5^ shows DO's expected sensitivity to physics beyond 
the standard model [164] for the already analyzed data and for the anticipated ac- 
cumulation of 6.8 fb^^ of Run II data. With this increase in integrated luminosity, 
the exclusion of certain models at 95% C.L. should be feasible. With more than 
6 fb^^, a measurement of \Vtb\ with an absolute uncertainty below 0.07 per experi- 
ment should be achievable as well. In addition, further refinements of the analysis 
techniques should facilitate improvements in precision beyond that expected just 
from the accumulation of more data. 



''After submission of this review, both collaborations announced 5 sd observations of single top 
quark production, see Refs. [397,398]. 
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Fig. 32: Left: Extrapolation of the single-top signal significance as a function of 
integrated luminosity from the 2.2 fb^^ CDF analysis result. Right: Expected 95% 
C.L. contours for a simultaneous measurement of the s- and t-channel single-top 
production rate by DO for different integrated luminosities. The standard model 
expectation is shown together with various other models. 



5.8.2. Polarization of the spin of the top quark 

As opposed to top quark pair production via the strong interaction, where the 
top quarks are produced essentially unpolarized (see Section 15. 5p , top quarks pro- 
duced singly via the electroweak interaction are expected to be highly left-handedly 
polarized [160]. The polarization of the top quark is reflected in the kinematic dis- 
tributions of its decay products, providing a test of the V — A structure of the Wtb 
coupling [137,399]. An observation of this polarization would also provide limits on 
the top-quark decay width and |Vfb|, since this would confirm that top quarks decay 
before depolarizing through QCD interactions. 

Both relevant single top quark production mechanisms at the Tevatron (s- and 
i-channel) exhibit up-type-down-type and tb quark lines interconnected by a VF 
boson. Since the W boson couples solely to fermions of left-handed chirality, in their 
rest frame single top quarks are highly polarized along the direction of the down- 
type quark [161,400,401]. (For the contributing 2 — > 2 processes this polarization is 
100%. These diagrams are related to the top decay with hadronic W boson decay 
via crossing symmetry, where the down-type decay products of the W boson exhibit 
optimal analyzing power.) The optimal spin basis for studying the single-top spin 
polarization therefore will use the direction of the down-type quark. 

For s-channel production, predominantly proceeding via ud tb, the antipro- 
ton beam is expected to provide the down-type quark most of the time. Indeed, 
measuring the top quark spin along the direction of the antiproton beam ( "antipro- 
ton basis" ) results in 98% of the top quark spins aligned in that direction. For the 
i-channel, the situation is slightly more complicated, since the down-type quark is 
contained in either the spectator jet or in one of the beams. With the largest con- 
tribution to the total production rate arising from ug — s- tbd, where the down-type 
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quark produces the (light quark) spectator jet, a reasonable choice for the spin basis 
is the spectator jet direction ("spectator basis"). Since the spectator jet is emitted 
in the forward direction, this basis is also still compatible with the cases where the 
down-type quark is in the initial state, resulting in 96% of the top quarks having 
their spins polarized along the light-quark jot direction. 

With the top quark decaying before it hadronizes, its spin information is passed 
on to its decay products. A straightforward observable is the angular distribution 
of the top quark decay product i in the top quark rest-frame: 



where 9\ is the angle between the decay product and spin-quantization axis, is 
the analyzing power describing the correlation between top quark spin and decay 
product and A-\i = (N-^ — Ni)/{N-i + Ni) is the spin asymmetry determining the 
magnitude of the observable angular correlations. The analyzing power is maximal 
(-1-1) for the down- type (T3 = — 5) decay products of the W boson (charged lepton, 
dr or s-quark), making the charged lepton the most sensitive and easily accessible 
spin analyzer. Using the spin quantization axes described above, the expected spin 
asymmetry is 0.96 for the s-channel and 0.93 for the t-channcl [161,400,401]. 

To perform a spin polarization measurement at the Tevatron, single top quark 
production in the f-channel is most promising due to its higher rate relative to 
the s-channel. The required integrated luminosity to observe spin polarization in 
the i-channel at the Tevatron was determined in a study, including effects of jet 
resolution and acceptance [162]. To establish the polarization at a level of five sd, 
~5 fb~^ of data will be needed, which should be available very soon. 

At the LHC, measurements of single-top spin polarization will benefit from the 
expected high-statistics single-top datasets, and optimal spin bases have already 
been explored for the two dominant production modes (t-channel [163], associated 
tW production [402]). Already with the first 2 fb~^, a polarization measurement 
with an uncertainty of 4% should be achievable based on the t-channel production 
alone [403]. 

5.8.3. Search for W bosons 

Electrically charged gauge bosons that are not part of the standard model are 
usually denoted as W . Such bosons are predicted in a variety of extensions of 
the standard model, incorporating larger gauge groups that reduce to the standard 
model at sufficiently low energies [164,404]. 

The most stringent limit to date in a direct search on the mass of such a W 
boson has been set by DO in the leptonic final state {W' £v) using 1 fb~^ of 
Run II data [405]. Assuming the W boson exhibits standard model W boson cou- 
plings to fermions, this search excludes the mass range below 1 TeV/c^ at 95% 
C.L. by studying the tail of the transverse mass [406] spectrum calculated from 




(30) 
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lepton transverse energy and ^t- Indirect W mass constraints are strongly model- 
dependent and range between lower limits of 549 GeV/c^ and 23 TeV/c^, being 
derived from (semi-) leptonic processes as well as from astrophysical and cosmolog- 
ical constraints [47]. 

A direct search for W' bosons in the hadronic final state {W qq') provides 
a less model-dependent measurement since both left-handed and right-handed W 
bosons can be observed in this final state, independent of any assumption about the 
mass of a right-handed neutrino m,y^ in the latter case. In contrast to this, the lep- 
tonic final state is only accessible for a right-handed W boson if the corresponding 
right-handed neutrino is not too massive {m^j^ < mw')- Searches for W' bosons as 
resonant structures in the dijet invariant mass spectrum have been carried out by 
UA2 [407] and at the Run I Tevatron by CDF [408] and DO [409]. 

Focusing on "hadronic" W' searches using third generation quarks in the final 
state, reduces the QCD multijet background compared to the (light) dijet final state 
searches. Such measurements are only sensitive to W bosons with masses above the 
tb threshold of «200 GeV/c^, but the low -mass range is excluded already by the 
current limits on single top quark production [410]. A W' signal would be observed 
as peak in the invariant mass distribution of its tb decay products (as usual, tb 
includes both tb and its charge-conjugate tb). 

Since the W —>■ tb decay mode contributes to s-channel single top production 
(see Section [321), these searches are based on the single-top production cross section 
analyses (see Section [5. 8. 1|) . For left-handed W' bosons, interference occurs with SM 
single top production, which is not the case for right-handed W bosons due to the 
different (right-handed) final state particles. Considering a right-handed W' boson, 
the decay width depends on the mass of the right-handed accompanying neutrino 
in leptonic decays. If m^^ > mw', only qq' final states are accessible, resulting in 
a width reduced by about 25%. This scenario generally results in a more stringent 
mass limit due to the enhanced tb branching fraction. A contribution of the W' 
boson to top quark decay is usually not considered due to its large mass. 

A first search for W' tb was performed by CDF in Run I, based on 0.1 fb^^ 
of lepton -I- jets data. At 95% C.L., lower limits on the mass of a right-handed W' 
boson were obtained, yielding 536 GeV/c^ for mi,j^ <C toh" and 566 GeV/c^ for 
mi,j, > mw' [411]. 

DO published a first search for W' tb in Run II, based on 0.2 fb^^ of lepton 
+ jets data and the corresponding single top cross section result [386]. For a right- 
handed W boson with CKM mixing equal to that of the SM, 95% C.L. lower 
mass limits of 630 GeV/c^ (670 GcV/c^) are obtained for m^j^ < mw' ("t-i/^ > 
mw). In addition, a first corresponding lower mass limit for a left-handed W 
boson is derived, taking the interference with SM production into account, yielding 
610 GeV/c2 [410]. 

Based on the 0.9 fb^^ lepton -I- jets dataset and the analysis from which the first 
evidence for single top production was obtained [21,22], DO obtains further improved 
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Fig. 33: Theoretical prediction at NLO and 95% C.L. limits for (7{pp W') x 
B{W' — > tb) versus mass of the W' boson. Left: Left-handed W boson. Right: 
Right-handed W' boson [412]. Results are from the DO experiment. 



W' mass limits [412]. Using the invariant mass of charged lepton, leading two jets 
and neutrino as a sensitive variable for separating signal and background, the 95% 
C.L. lower mass limit for a left-handed W boson interfering with SM single-top 
production increases to 731 GeV/c^. For a right-handed W' boson, the 95% C.L. 
lower mass limits are 739 GeV/c^ (768 GeV/c^) for m^^ < mw {my^ > mw), as 
illustrated in Fig. [33l The latter two cross section limits correspond to upper limits 
on the W gauge coupling in units of the SM weak coupling of 0.72 (0.68) for a W 
boson mass of 600 GeV/c^. The dominant systematic uncertainties included in these 
limits are the theoretical cross sections (affecting the background normalization) and 
uncertainties on jet energy calibration and b jet simulation (affecting background 
normalization and distribution in the sensitive variable). 

CDF has obtained a preliminary result for their W' tb search, based on 
1.9 fb~^ lepton -I- jets Rim II data, using the invariant mass of the reconstructed 
W boson and the two leading jets as sensitive variable. 95% C.L. lower limits on 
the mass of a right-handed W' boson are found to be 800 GeV/c^ for m^^ < mw' 
and 825 GeV/c^ for m^^ > mw'- Neglecting interference effects the former limit 
is considered to apply for a left-handed W boson as well. The corresponding W' 
gauge coupling in units of the SM weak coupling is found to be <0.68 and <0.63, 
respectively, for a W boson mass of 600 GeV/c^ [413]. 

The more general case of a W boson with an admixture of left- and right-handed 
couplings to SM fermions has not been studied thus far. 



5.8.4. Search for single top production via charged Higgs bosons 

The standard model Higgs sector, with its single Higgs doublet of complex scalar 
fields to break electroweak symmetry and generate masses of weak gauge bosons and 
fermions (see Section [2TT|) can be easily extended to include a second Higgs doublet, 
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Fig. 34: Left: Distribution of the invariant mass of reconstructed W boson and two 
jets in 0.9 fb~^ lepton + jets data, SM background processes, and charged Higgs 
boson signal in a Type-Ill THDM for several to^± , with the expected rate enhanced 
by a factor of 50. Right: 95% C.L. exclusion region in the (m^/i , tan/3) plane for a 
Type-I THDM. If exceeds 50 GeV/c^, the analysis is no longer valid, and no 
limits can be derived as illustrated by the darker area [417]. 



resulting in "Two Higgs Doublet Models" (THDM or 2HDM) [414,415]. In contrast 
to the single neutral scalar CP-even Higgs boson predicted by the SM, THDM 
give rise to five physical scalar Higgs bosons after electroweak symmetry breaking. 
Two of these are charged bosons {H^), providing a unique signature for physics 
beyond the standard model. Three different Higgs-fermion couplings are discerned 
in THDM. Type-I models provide coupling of only one of the Higgs doublets to 
fermions. In Type-II models, each of the doublets couples solely to up-type fermions 
and down-type fermions, respectively, while in Type-Ill models general couplings 
of both Higgs doublets to fermions are allowed. In the latter case, Higgs-mediated 
flavor-changing neutral currents at tree level must be sufficiently suppressed to be 
compatible with experimental limits. This can be achieved through an appropriate 
choice of the Higgs parameters [416]. 

If the charged Higgs boson is heavier than the top quark (m/f± > rrit), its 
production via quark fusion can contribute to single top quark production through 
the decay into third-generation quarks: qq' tb. Due to mass dependent 

couplings of the charged Higgs boson, this decay is dominant in many models. The 
signature of this process is identical to that of s-channel single top-quark production, 
and the search for charged Higgs bosons can be performed similar to that for W 
bosons, with the simplification that interference with the SM production process 
can be neglected. 

DO performs a first direct search for the process qq' th ^ Evbh [417], 

based on the analysis providing first evidence for single top production [21,22] in 
a 0.9 fb~^ lepton + jets dataset. Restricting the jet multiplicity in the events to 
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exactly two jets, corresponding to the s-channel final state, charged Higgs bosons are 
sought in the mass range 180 GeV/c^ < mu± < 300 GeV/c^ for all three types of 
THDM. The sensitive variable used to discriminate the charged Higgs boson signal 
from SM background processes is the invariant mass of the reconstructed W boson 
and the two jets, as illustrated in Fig. [Mj Since no evidence for signal is observed 
in the data, 95% C.L. upper limits on the charged Higgs boson production cross 
section multiplied by branching fraction into third generation quarks are provided 
for all three types of THDM. The dominant systematic uncertainties result from 
the jet energy scale calibration, modeling of the b jet identification and theoretical 
uncertainties in modeling and normalizing the signal. For the Type-I THDM, the 
limits are translated into a 95% C.L. exclusion region in (TO/j±,tan/3) parameter 
space (see Fig. [34|) . where tan/3 is the ratio of the vacuum expectation values for 
the two Higgs doublets. 

More searches for charged Higgs bosons, especially those in top quark decays 
(for m^± < mt), are described in Section [6.51 

5.8.5. Search for single top production through neutral currents 

Single top quark production via flavor changing neutral interactions of light u, c 
quarks and the Z, 7, g gauge bosons is possible in the standard model through 
higher-order radiative corrections, but so strongly suppressed that it cannot be ob- 
served. Consequently, searches for these production mechanisms at tree level probe 
for corresponding anomalous coupling strengths k [164,418] that are predicted by 
various extensions of the standard model [419]. 

The processes involving photon or Z boson exchange have been extensively stud- 
ied at LEP and HERA. At both accelerators top quarks can only be produced singly 
at the available center of mass energies due to the large top quark mass. 

At LEP, single top quark production proceeds via the SM process e"^e^ 
e~Vetb, which can be ignored in the available datasets due to its tiny production 
rate. All four LEP experiments searched for single top production via e'^e~ — > 
tc/tu in both hadronic and semileptonic final states, resulting from the different W 
boson decay modes from the top quark decay. While only the SM decay t — > Wh is 
considered, a possible reduction of its branching ratio due to possible FCNC decays 
is accounted for when deriving the results. Since no evidence for single top quark 
production is observed, 95% C.L. upper limits on the cross section are extracted, and 
corresponding model-dependent upper limits on the anomalous coupling parameters 

and Kz are determined [420-423]. 

Single top quark production at HERA is possible via the charged current SM 
process ep vthX that has a negligible production rate here as well. Both HI and 
ZEUS have searched for the inclusive neutral current production of top quarks in 
ep etX. Because of the large Z boson mass, this reaction is most sensitive to 
couplings involving photons. Due to the large proton momentum fractions needed 
for single top production, the u quark contribution will dominate over that from 
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Fig. 35: 95% C.L. upper limits on the top quark anomalous couplings to photons 
and Z bosons Ktu-y and vtuz [424]. Depicted are the limits from CDF [425,426], 
L3 [422], HI [424] and ZEUS [427]. 



the c quark (see Fig. [5|) , resulting therefore in highest sensitivity to tuj couplings 
at HERA. 

Using 0.1 fb^^ of integrated luminosity, and assuming the SM top quark decay 
t Wb, HI and ZEUS search for single top production both in the leptonic and 
hadronic W boson decay channel. ZEUS observes good agreement with the SM pre- 
diction, and sets 95% C.L. upper limits on single-top production rate and on the 
FCNC magnetic coupling Ktwy and vector coupling vtuZ, neglecting charm contri- 
butions [427]. HI observes five events in the leptonic channel, with an expected SM 
background contribution of 1.31 ± 0.22 events, while the hadronic channel exhibits 
no excess over the standard model prediction. These two channels are compatible 
at the 1.1 sd level, and both a combined single-top cross section with about 2 sd 
significance and 95% C.L. upper limits on the cross section and on Ktu-y (assuming a 
statistical fluctuation in the data) are provided [428] . In a recent preliminary update 
of the measurement in the leptonic channel by HI, using an integrated luminosity 
of 0.5 fb~^, good agreement with the standard model expectation is observed and 
improved limits on the single-top cross section and Ktwy are obtained [424]. 

Further limits on the anomalous couplings and kz (inclusive for u, c contribu- 
tions) have been measured at the Tevatron by CDF via a search for neutral-current 
top quark decays t — > 7(7 and t — > Zq as will be discussed in Section [6.31 The most 
stringent results on anomalous top quark couplings involving photons and Z bosons 
obtained at LEP, HERA and the Tevatron are summarized in Fig. [35l with the ex- 
ception of the latest CDF limit on t ^ Zq decays [429] that constrains kz better 
than the limit from L3. 

Flavor changing neutral-current (FCNC) couplings of top quarks and gluons 
have not been studied as extensively. A constraint on the anomalous gluon coupling 
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Fig. 36: Left: Neural Network (NN) discriminant distribution in 0.2 fb^^ lepton 
+ jets data, with simulated FCNC signal increased by a factor of ten, and SM 
background. The signal distribution represents the sum of tug and teg processes, 
evaluated for Ktqg/h. = 0.03 TeV"-*^. Right: Exclusion contours for anomalous top- 
gluon couplings for different levels of confidence [432]. 

Ktqg/h. of < 0.52 TeV~^, where A gives the scale for new physics, was extracted 
from the observed tt pair production cross section at the Run I Tevatron and a 
possible new physics contribution that could be still accommodated within two sd 
of combined experimental and theoretical uncertainties [430] . Another limit on the 
anomalous gluon coupling was obtained using the single-top production cross section 
limit measured by ZEUS [427]. Neglecting any effects that arise from the different 
final states obtained in the gluon channel compared to the original search (one 
additional light jet is present in the gluon case), at 95% C.L. Ktqg/A < 0.4 TeV~^ 
is obtained [431]. 

DO has performed a first search for single top production via flavor changing 
neutral-current couplings to gluons at a hadron collider, using 0.2 fb~^ lepton -|- 
jets data [432]. The analysis is based on the corresponding search for SM single-top 
production [386,387], but is restricted to events with only one b tagged jet, and 
treats s- and t-channel SM single-top production as background. Since the neutral 
current decays t gu/gc exhibit a negligible branching fraction for Ktqg/h. ^ 
0.2 TeV~^ [433], exclusively the standard model top quark decay can be considered. 

To separate the FCNC signal from the overwhelming SM background, a neural 
network is deployed, with ten input variables based on global event kinematics, 
angular correlations and kinematics of the individual reconstructed objects. The 
resulting data distribution is shown in Fig. 1361 and exhibits good agreement with 
the SM prediction, which provides limits on the FCNC couplings Kjug/A and Kicg/A. 

Systematic uncertainties affecting either the normalization or both normaliza- 
tion and shape of the distributions are taken into account when calculating the 
two-dimensional Bayesian posterior probability density, resulting in the exclusion 
contours for the two couplings for different confidence levels shown in Fig. 1361 The 



84 Marc-Andre Pleier 



c 

> 
m 

CO 
■o 

'■B 
c 
n 
O 



W+Uet ITag CDF II Preliminary 2.2 fb 



300 



200 



100 




-0.5 



0.5 
NN Output 



^ 3 
Q. 



CDF II Preliminary 2.2 fb 



CO 

X 

1° 

D 



1 1 1 1 1 1 1 1 1 1 1 1 

—theory, NLO 

— upper limit @ 95% C 

□ excluded 


L. / 


.^-r<. 1 .... 1 . . 





0.01 0.02 0.03 0.04 



VA [TeV^] 



Fig. 37: Left: Neural network discriminant output in + one b jet events observed 
in 2.2 ib~^ of CDF data (black points) and standard model background processes. 
A FCNC single top signal with the observed 95% C.L. excluded production rate 
is added in red. Right: Extraction of the ntug/-^ limit from the intersection of the 
observed cross section limit and theoretical rate prediction, assuming ntcg = [434] . 



largest normalization uncertainties arise from the background cross section uncer- 
tainties, which includes the uncertainty on the top quark mass for ti and single 
top samples. The largest uncertainties affecting the shape as well arise from jet 
energy scale calibration and b tag modeling. 95% C.L. upper limits on Kt„g/A and 
Ktcg/-^ are obtained by integrating over the other variable, and yield 0.037 TeV~^ 
and 0.15 TeV~^ for tug and teg couplings, respectively. These limits represent a 
significant improvement over previous values by up to an order of magnitude. 

CDF also reports a recent preliminary search for FCNC single top production, 
based on 2.2 fb~^ of data [4340. In distinction to the DO analysis, where 2 — > 2 
teg and tug signal processes are considered, CDF investigates the 2 — > 1 processes 
u{c) + g ^ t. Since also in this analysis only SM top quark decay is considered, 
events with one isolated energetic lepton, and one 6 tagged jet are selected. 

Signal and SM background processes are separated using a Bayesian neural 
network based on 14 input variables containing information from the reconstructed 
objects and event kinematics. In a template fit to the observed distribution in data 
good agreement is found with the SM background, as illustrated in Fig. [371 and a 
95% C.L. upper limit on FCNC single top production via u{c) + .g ^ t of 1.8 pb is 
obtained in accordance with the expected sensitivity. 

Based on LO predictions for the FCNC signal process from TopReX [436] , and 
using NLO fc-factors [437-439], the obtained cross section limit can be converted into 
limits on anomalous gluon couplings. No two-dimensional information is used in this 



An updated version of the result has been pubhshed after completion of this review, see Ref. [435] . 
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analysis for contributions of teg relative to tug signal processes, and one coupling 
is assumed to vanish in deriving the limit for the other coupling. The resulting 95% 
C.L. upper limits are Kt«g/A < 0.025 TeV"! (see ¥ig.^ and Ktcg/^ < 0.105 TeV-^ 
(not shown). 

5.8.6. Anomalous Wtb couplings in single top production 

The couplings between quarks and electroweak gauge bosons were directly scru- 
tinized at LEP [440], with the exception of the top quark. At the Tevatron, the 
couplings of the top quark and the W boson can be studied in measurements of top 
quark decay properties in ti production (see for example Section [6.1|) and via single 
top quark production. Physics beyond the standard model could modify the Lorentz 
structure of the Wtb vertex. Considering a more general extension of the standard 
model Wtb interaction Lagrangian, as discussed in Section [3.3.3l new physics could 
introduce contributions from right-handed vector (//^) and left- and right-handed 
tensor couplings (/^, ff'), in addition to the pure left-handed vector coupling (/f ) 
of the standard model. 

DO has published constraints on such extended Wtb interactions, including the 
first direct limits on the left- and right-handed tensor couplings [441], based on 
the analysis that provided first evidence for single top production in 0.9 fb~^ of 
lepton -I- jets data [21,22]. In the analysis, single top quark production and decay 
are considered to take place only via W bosons, with the dominant contribution 
arising from the Wtb interaction, which is assumed to be CP conserving. Anomalous 
couplings at the Wtb vertex can modify both the total single-top production rate 
and the observed kinematics in the events relative to SM expectation [137,442-444]. 
The latter effect is illustrated in Fig.[35]for the charged-lepton transverse momentum 
distribution. 

Since a simultaneous fit of all four couplings to data is not feasible with the 
available statistics, the SM coupling and one additional anomalous coupling contri- 
bution at a time is considered in varying proportions (with the remaining two other 
anomalous couplings set to zero). The resulting scenarios are denoted as (Li,i?i), 
(ii,i?2) and (Li,i2)- Non-negligible interference effects in the last case are also 
taken into account. For signal discrimination from SM background, boosted deci- 
sion trees are used that are based on the same variables as used in Refs. [21,22], but 
with the lepton pT distribution added. One example distribution for the (Li,i?2) 
case is shown in Fig. 1381 

The decision tree output in data is compared with the various single-top signal 
models in the twelve subchannels defined by lepton flavor (e, ^), jet multiplicity 
(two, three, four) and b tag multiplicity (one, two). This yields a two-dimensional 
Bayesian posterior probability density, depending on \ fi\'^ and the anomalous cou- 
pling |/ano|^ considered in the respective scenarios. Systematic uncertainties are 
taken into account, with dominant contributions arising from background normal- 
ization, modeling of b tagging and jet energy scale calibration. The latter two affect 
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Fig. 38: Left: Charged-lepton pT distribution in 0.9 fb^^ lepton + jets data (two 
jets and one b tag) and the corresponding SM single-top signal and background 
contributions. The effect of each of the four different Wtb couplings on the signal is 
also shown (calculated with the other couplings set to zero), with the normalization 
enhanced by a factor of ten. Right: Boosted decision tree output for the same data, 
signal and background contributions, with the (ii,i?2) scenario couplings overlaid 
(normalization increased by a factor of five) [441]. 



Table 15: Total single-top production rates obtained in three anomalous coupling 
scenarios, together with the corresponding one-dimensional measurements and lim- 
its for the selected couplings [441]. 



Scenario 


Cross Section {tb + tqb) 


Considered Couplings 


(ii,i?i) 


5.2tli pb 


l/f |2 < 2.5 (95% C.L.) 






|/2^|2 < 0.3 (95% C.L.) 


{Li,L2) 


AAtli pb 


|/2^|2 < 0.5 (95% C.L.) 



both normalization and shape of the simulated spectra. The maxima of the likeli- 
hoods in all three considered scenarios yield zero for |/ano|^, and 95% C.L. upper 
limits on these anomalous couplings are provided from the one-dimensional likeli- 
hood projections. These results are summarized together with the measured single- 
top production rates and [/f [^ values obtained from one-dimensional likelihood 
projections in Table [121 The SM Wtb interaction is favored over any anomalous 
alternative studied. This analysis will greatly benefit from the increased statistics 
already in hand and from more expected until the end of Run II. 
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6. Decay Properties of the Top Quark 

The previous chapter demonstrated that no significant deviations from the standard 
model expectations for top quark production via the strong or electroweak interac- 
tion have thus far been observed. In this chapter, decay properties of the top quark 
will be investigated based on tt data, generally assuming that top quark production 
proceeds according to the standard model. 

6.1. Measurement of the W boson helicity in tt decays 

The helicity of the W boson in top quark decays can be used to test the V—A Lorentz 
structure of the Wtb interaction (see Section fS.S.Sp . According to the expectation 
from the standard model, W bosons from top quark decays should be longitudinally 
polarized with a fraction /o « 70% and left-handed with a fraction /_ « 30%. The 
right-handed fraction /+ is strongly suppressed, and below the per mill level [189]. 
For the decay of antitop quarks, the CP conjugate statement is implied, resulting 
in W~ bosons from i decays with either longitudinal or right-handed polarization. 

A pure V + A structure of the Wtb interaction would result in an observa- 
tion of a right-handed fraction /+ — 30%, with negligible left-handed contribution. 
Small V + A admixtures to the SM left-handed charged-current weak interaction 
are predicted, for instance, within SU{2)ii x SU{2)l x U{1)y extensions of the 
SM [445-448]. Such contributions would result in an enhancement of /+ while not 
significantly affecting /q. Since the decay amplitude for longitudinally polarized W 
bosons is proportional to the top quark's Yukawa coupling [174], /o is sensitive to 
the mechanism of EWSB, and would be altered, for example, in topcolor-assisted 
technicolor models [191,449]. 

The radiative decay rate & — > 57 can be used to set indirect limits on the -f A 
admixture in top quark decays to below a few percent [450-453], assuming there 
are no contributions from gluonic penguin diagrams in addition to the electroweak 
ones. This section will discuss the direct measurements of the W boson helicity 
performed at the Tevatron using lepton -I- jets and dilepton datasets. 

Thus far, four analysis techniques have been deployed to extract the W boson 
helicity fractions, based on: 

(i) helicity angle (cos 9*): The helicity of the W boson is reflected in the angular 
distribution cos 9* of its decay products, with 6* being the angle of the down- 
type (T3 = — ^) decay products of the W boson (charged lepton, d or s quark) 
in the W boson rest frame relative to the top quark direction [186-189]: 

= /- • ^(1 - cosr )2 -t- /o • J(l - cos^ 9*) + U ■ ^(1 + cosr )^(31) 

where /_ can be replaced by (1 — /+ — /o). The resulting distributions for each 
helicity fraction and the superposition expected from the standard model are 
shown in Fig. 1391 A measurement of cos 9* provides the most direct measure- 
ment of the W boson helicity, but it requires the reconstruction of the momenta 
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Fig. 39: Helicity angle 008 6** distributions for left-handed, longitudinal and right- 
handed W bosons. The superposition expected from the standard model is shown 
as well. 



of the top quark and W boson, which is challenging and involves using -^t, which 
has a rather poor resolution. 

(ii) charged-lepton px spectrum (p^ ): The helicity of the W boson is correlated with 
the charged-lepton momentum distribution: Since the from W'^ decays are 
left-handed, while the are right-handed, in case of a left-handed W'^ decay 
the are preferentially emitted opposite to the momentum vector of the . 
This leads to a softer p^ spectrum in the laboratory frame. Conversely, the 
are preferentially emitted along the direction of the momentum in case 
of a right-handed decay, leading to a harder p^ spectrum. The £^ from 
longitudinal decay represent an intermediate case (see Fig. |40l) . 

(iii) squared invariant mass of b quark and charged lepton (Mfj^): In the limit of 
mi, — 0, the helicity angle distribution cosO* can be approximated using the 
squared invariant mass of the system composed of the b quark and the charged 
lepton M|j: 

cosr = ^i:f^i-4^.^%-i. (32) 

|p,||pb| m^-M^ 

This way one avoids the challenge of kinematic reconstruction of the top quark 
and the application of by using only momenta measured in the laboratory 
frame. 

(iv) Matrix Element method (ME): The Matrix Element method was originally de- 
veloped by DO [454] , yielding a very precise mass measurement given the limited 
Run I data sample (see also Section I7.3.2p . Using all the available kinematic 
information in each event, a probability for the event to correspond to a tt final 
state as a function of the helicity of the W boson can be calculated, based on 
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Fig. 40: Distributions of reconstructed charged-lepton pr and cosO* (based on 
Eq. [32|) for top quark decays involving left-handed, right-handed and longitudinally 
polarized W bosons [455]. 

the LO matrix element. 

The following subsections will give brief examples for each method, followed by a 
summary of the current status of the measurements. 



6.1.1. and M^j, 

CDF has measured the W boson helicity in a 0.2 fb^^ Run II dataset using the 
charged-lepton {Pt) ^"^^ squared invariant mass of the b quark and charged 
lepton (Mfi^) to approximate cos 9* [455]. The dependence of these observables on 
the W boson helicity for a top quark mass of 175 GeV/c^ after event selection 
and reconstruction is shown in Fig. 1401 Since the world-averaged top and W boson 
masses are used for calculating cos 0* , rather than the corresponding event-by-event 
reconstructed masses that would smear out the distribution due to the larger inher- 
ent uncertainties, values are observed outside of the physical range — 1 < cos 9* < 1. 

For the p^ analysis, a b tagged lepton -I- jets sample is used requiring at least 
three jets, and yielding 57 events of which approximately 2/3 are tt signal. In addi- 
tion a dilepton sample, with a minimum value for the scalar sum of the transverse 
energy of jets, leptons and fir is analyzed, yielding 13 events with a signal fraction 
of «0.79. The analysis based on M|j uses the lepton -I- jets sample alone, requiring 
a fourth jet and a good kinematic fit to the ti hypothesis for a top mass of 175 
GeV/c^. This provides the lepton and appropriate jet to form M^f^, and leaves 31 
events for this analysis with a signal fraction of wO.78. 

For both analyses, the data distributions are fitted separately to p^ and cos 9* 
templates of signal with the different W boson helicities and background. Because 
of limited statistics, the helicity fractions /o and /+ cannot be fitted simultane- 
ously. Consequently, /o or /+ are constrained to their standard model values when 
fitting for /_|_ or /o, respectively. Both analyses are finally combined, taking sta- 
tistical and systematic correlations into account, and yield results consistent with 
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Table 16: Results of individual and combined measurements of /o and /+ using Mf^ 
and py. N indicates the number of events or leptons used. If two uncertainties are 
given, the first is statistical and the second systematic. For the combined results, 
the statistical and systematic combined uncertainty is given. For the p^{££) result, 
an observation of < —0.54 is expected 0.5% of the time for the SM /o value of 
0.7 [455]. 



Analysis 


N 


/o 


/+ 


Ml 


31 


0.99t°j9 ±0.19 


0.23 ±0.16 ±0.08 




26 


-0.54l°f5±0.16 


-0.47 ±0.10 ±0.09 




57 


0.95+!^ ±0.17 


o.iil:gj^±o.io 




83 


0.3lt°J^±0.17 


-o.i8t|^:i2±o.i2 


Combined 




0.74±°;g 


oo+°-2" 


95% C.L. limit 




< 0.95,> 0.18 


< 0.27 



the standard model expectation, as shown in Tabic 1161 The dominant systematic 
uncertainties come from uncertainties on the top quark mass, background shape 
and normalization, effects of initial- and final-state radiation (ISR/FSR) and the 
PDFs. 

CDF has also measured the fraction of right-handed W bosons assuming /q to 
be 0.7, using the Ml method on a 0.7 fb^^ Run II dataset [456]. Using a single and 
double h tagged lepton + jets in addition to a dilepton dataset, /+ is extracted via 
maximum likelihood fits of the Ml distributions in data ioV and V — Ati signal 
Monte Carlo and background contributions. Including uncertainties on the tt sig- 
nal and background cross sections, the lepton + jets sample yields /+ = 0.06 ± 0.08, 
while the dilepton sample gives /+ — —0.19 ± 0.11, corresponding to a compatibility 
of the measurements at the level of 2.3 sd. A combination of these measurements, 
including all systematic uncertainties, yields /+ = — 0.02 ± 0.07(stat ± syst), cor- 
responding to /+ < 0.09 at 95% C.L. The main contributions to the systematic 
uncertainty come from the jet energy scale, background shape and normalization 
and limited Monte Carlo statistics. 

6.1.2. Matrix element method 

DO has used the matrix element method that was originally employed to measure 
the top quark mass [454] to extract the longitudinal W boson helicity fraction from 
0.1 fb^^ of Run I data [457]. The selected lepton + jets event sample corresponds 
to that of the preceding mass analysis [320] and comprises both soft muon b tagged 
events and untagged events, which have additional kinematic requirements, yielding 
a total of 91 events. 

By comparing the measured set of four- vectors in each event with the differential 
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Fig. 41: Two-dimensional probability density observed in 0.1 fb ^ lepton + jets 
Run I data as a function of /o and top quark mass m* [457]. 

cross section for ti signal and the dominant VT+jets background, /o can be extracted 
by fixing /+ to its SM value, allowing the ratio /o//- to vary. The use of both W 
boson decays per signal event increases the statistical sensitivity of the method. 
Since the calculation of signal and background probabilities is based on leading- 
order matrix elements, only events with exactly four jets are accepted, reducing the 
sample to 71 events. In order to increase signal purity, a cut on the background 
probability is applied, leaving 22 events to be analyzed (as in the corresponding 
matrix-element mass analysis [454]), with a signal to background ratio of 12/10. 

To take the dependence of the /o measurement on the top quark mass into ac- 
count, a two-dimensional likelihood, depending on /o and mt, is calculated and cor- 
rected for response deviations from unity for different /o input values (see Fig. I4ip. 
Since statistics are insufficient, a simultaneous optimization for both observables is 
not feasible; instead /o is evaluated by integrating over the top quark mass, for a 
range between 165 and 190 GeV/c^. The maximum in the probability yields the 
central value of the measurement, with the 1 sd uncertainty band corresponding 
to a convolution of statistical and top quark mass uncertainties. Other systematic 
uncertainties from acceptance and linearity of response or jet energy scale are small 
compared to this, yielding the final result: 

/o = 0.56 ± 0.31(stat ® mt) ± 0.07(syst). (33) 



CDF has obtained a preliminary result for /o using the same method on 1.9 fb"^ 
of Run II data, assuming a fixed top quark mass of 175 GeV/c^ [458]. 468 events are 
selected in a 6 tagged lepton -I- jets sample with at least four jets (only the leading 
four are used in the analysis) and a minimum value for Ht, yielding a signal fraction 
of about 0.84. Fixing /+ to its SM value, the longitudinal W boson helicity fraction 
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is found to be /o = 0.64 ± 0.08(stat.) ± 0.07(syst.), with the dominant systematic 
uncertainty coming from the Monte Carlo generator used (pythia versus herwig) 
for the cahbration of the measurement, /o is found to change by =f0.035 for a ±2.5 
GeV/c^ variation in the top quark mass. Thus far, the analysis has not yet been 
extended to measure /o and or /o and mt simultaneously. 

6.1.3. Helicity angle cosO* 

DO has published a first model-independent measurement of the W boson helicity 
fractions by comparing the cos 9* distribution in data to templates of background 
and right-handed, left-handed or longitudinal W bosons in ti signal, using /+ and 
fo as freely floating parameters and /_ = 1 — /+ — /o [459]. In a 1 fb^^ dataset, 
lepton -I- > 4 jets and dilepton events are selected. The signal purity is increased 
in each subsample by a cut on an individually optimized multivariate likelihood 
discriminant based on event kinematics and the output of a neural network 6-tagging 
discriminant. The cut values are chosen in each subsample to yield the best expected 
precision for the helicity measurement. 

The statistical sensitivity of the analysis is further improved by about 20% by 
including the W — *■ qq' decays in the lepton + jets sample in the measurement. 
This is accomplished through picking one of the W boson daughter jets at random 
for the calculation of cos 9* , which introduces a sign ambiguity. Consequently, only 
I cos 9* I is considered which does not permit to discriminate left- from right-handed 
W bosons, but still adds information on /q. 

The four-momenta of the top quarks and W bosons in the lepton + jets sample 
are reconstructed based on the best kinematic fit to a top quark event hypothesis 
for mt = 172.5 GeV/c^, using the leading four jets to obtain cos9* and |cos0*|. 
For the kinematically less constrained dilepton events, a top quark mass of mt = 
172.5 GeV/c^ is assumed and the kinematics solved up to a four-fold ambiguity 
in addition to the two-fold ambiguity from the lepton-jet pairing (only the leading 
two jets are used). Jet and lepton energies in each event are smeared within their 
resolutions to explore the phase space consistent with the observed values. The 
average of the obtained cos 9* values is then used for each charged lepton, providing 
two measurements per event. The resulting distributions are shown in Fig. 1421 Note 
that due to reconstruction effects the shape of the standard model expectation 
differs from the theoretical prediction in Fig. [39l 

A template fit of these distributions yields /o — 0.425±0.166(stat.)±0.102(syst.) 
and /+ — 0.119 ±0.090(stat.)±0.053(syst.). The result is compatible with the stan- 
dard model expectation at 30% C.L. It should be noted that the individual mea- 
surements in the lepton -I- jets and dilepton channels differ by about 2.1 sd [460]. 
The major systematic uncertainties on the measurement are summarized in Ta- 
ble 1171 with the largest uncertainty arising from it signal modeling, evaluated 
through varying the Monte Carlo generators used (pythia versus ALPGEn), from 
changing underlying event models to estimate the effects of gluon radiation, and 
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cos6* 

Fig. 42: Helicity angle distributions in lepton + jets (a,b) and dilepton events (c). 
Points with error bars represent 1 fb~^ of DO data. The solid open histograms show 
the result of the model independent fit described in the text, while the dashed open 
histograms show the standard model expectation. The filled histograms represent 
the background contribution [459]. 

Table f7: Major systematic uncertainties on the simultaneous measurement of /o 
and /+ by DO in 1 fb^^ of data [459]. 



Source 


Uncertainty (/q) 


Uncertainty (/_(-) 


Top mass 


0.009 


0.018 


Jet reconstruction eff. 


0.021 


0.010 


Jet energy calibration 


0.012 


0.019 


6 fragmentation 


0.016 


0.010 


ti model 


0.068 


0.032 


Background model 


0.049 


0.016 


Template statistics 


0.049 


0.025 


Total 


0.102 


0.053 



restricting the samples to contain only one primary vertex to study the sensitiv- 
ity of the measurement to variations in instantaneous luminosity. Constraining 
/o or /+ to their SM values, when fitting for /+ or for /o, respectively, yields 
/o = 0.619±0.090(stat.)±0.052(syst.) and /+ = -0.002±0.047(stat.)±0.047(syst.), 
in agreement with expectations from the SM. 

In a recent prehminary update, DO has added 1.2 fb^^ lepton + jets and 1.7 fb^^ 
dilepton (e/z only) data [461] to the above analysis. The model-independent fit 
for the combined data yields /o = 0.490 ± 0.106(stat.) ± 0.085(syst.) and /+ = 
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CDF II preliminary, 1 .9 fb' | 




Fig. 43: Left: Result of the model-independent W boson helicity fit by DO [461]. 
The ellipses show the 68% and 95% C.L. contours around the measured data point. 
The star shows the SM expectation, and the triangle denotes the physically allowed 
region where fo and /+ sum to < 1. Right: Two-dimensional 95% C.L. exclusion 
area in the (/+, fo) plane measured by CDF [462]. 

0.110 ±0.059(stat.)±0.052(syst.), consistent at 23% C.L. with the standard model 
(see Fig. l43l) . The results from the lepton + jets and dilepton channels remain 
marginally consistent with a p-value of 1.6%. 

CDF has obtained two preliminary results for a model-independent extraction 
of W boson helicity in 1.9 fb~^ of data. These are based on using only the charged 
lepton in b tagged lepton -I- > 4 jets events to obtain cos6'* [462,463]. A combination 
of both results also has become available [464jf|. The measurements are compatible 
with the standard model expectation, with each other, and with the DO measure- 
ments presented above and are summarized together with other results in Table [T51 
The two-dimensional 95% C.L. exclusion area in the (/+,/o) plane measured by 
CDF [462] is shown in Fig. [H 

6.1.4. Summary 

All available measurements of the W boson helicity performed thus far in tt events 
at the Tevatron are compatible with the standard model expectation, and are sum- 
marized in Table [TSl 

The sensitivity of the measurements in Run I and initial studies in Run II 
only allowed model-dependent measurements of any single helicity fraction to be 
performed at a time, while fixing the other fraction to its standard model value. 
However, with the large amount of data available by now, a simultaneous extraction 



^This result has been published after completion of this review, see Ref. [465]. 
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Table 18: W boson helicity measurements in ti events performed thus far at the 
Tevatron, with their integrated luminosities, data selections [Ij = lepton + jets, U 
= dilepton) and analysis methods used. Model independent results are indicated by 
a yes in the "Ind. fit" column. The three analyses using 0.1 fb~^ are from Run I; the 
analyses using more than 1 fb~^ are preliminary. Ref. [464] provides a combination 
of the results from Rcfs. [462,463]. 



J Cdt 

[lU J 


Sel. 


fo 


/+ 


Ind. 

fit 


/+ < 


Method 


Ref. 


0.1 




0.91 ±0.37± 0.13 


0.11 ±0.15 


no 


0.28 


Pt 


[466] 


U.i 


DA DP 




A no _i_ nil 


no 


n 1 Q 
U.io 




[467] 


n 1 
U. 1 


DA 


U.OD ± U.oi 




no 




A/TTT 


[45 (\ 


0.2 




n 74+0-22 
'^-0.34 


o.oot°:?^ 


no 


0.27 


Ml, p'j, 


[455] 


0.2 


ej 




0.00 ±0.13 ±0.07 


no 


0.25 


cose* 


[468] 


0.3 


ej 


OMto.ll ± 0-06 


0.05l[!:o5 ± 0.03 


no 


0.26 


cos 6* 


[469] 


0.4 






0.06 ±0.08 ±0.06 


no 


0.23 


cos 9* 


[470] 


0.7 






-0.02 ±0.07 


no 


0.09 


Ml 


[456] 


1.0 




0.62 ±0.09 ±0.05 
0.43 ±0.17 ±0.10 


0.00 ±0.05 ±0.05 
0.12 ±0.09 ±0.05 


no 
yes 




cose* 


[459] 


1.9 


ij 


0.64 ±0.08 ±0.07 




no 




ME 


[458] 


1.9 




0.59 ±0.11 ±0.04 
0.65 ±0.19 ±0.03 


-0.04 ±0.04 ±0.03 
-0.03 ±0.07 ±0.03 


no 

yes 


0.07 


cos 9* 


[462] 


1.9 




0.66 ±0.10 ±0.06 
0.38 ±0.21 ±0.07 


0.01 ±0.05 ±0.03 
0.15 ±0.10 ±0.05 


no 
yes 


0.12 


COS0* 


[463] 


1.9 




0.62 ±0.11 
0.66 ±0.16 


-0.04 ±0.05 
-0.03 ±0.07 


no 
yes 




cose* 


[464] 


2.7 




0.49 ±0.11 ±0.09 


0.11 ±0.06 ±0.05 


yes 




cos 8* 


[461] 



of /+ and /o for W bosons is possible without constraint, except for unitarity 
(/+ + /o + /- = 1)- Such measurements will also clearly benefit from increased 
luminosity. A combination with the measurement of the single-top production cross 
section in the s- and i-channcl will help to fully specify the tWb coupling [191] (see 
Section [3331§. 

The model-dependent measurements, where one of the helicity fractions is fixed 
to its standard model value, have reached a considerable precision with statistical 



'^A first such measurement was published by DO [471] after completion of this review. 
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approaching systematic uncertainties. 

It is interesting to note that discrepancies of > 2 sd between results from the 
dilepton and lepton + jets samples have been observed both at CDF and DO using 
different analysis techniques. This deserves further scrutiny in future analyses. 



6.2. Measurement of B{t Wb)/B{t Wq) 

As described in Section FS.S.li in the standard model framework the top quark decays 
almost exclusively into a W boson and a b quark due to the dominant corresponding 
CKM matrix element Vtb- The ratio R of the top quark branching fractions can be 
expressed via the CKM matrix elements as: 

^ Bjt^Wb) ^ I Vtb ? 

q—d.s,b 

Measuring R provides therefore the relative size of |Vtfc| compared to \Vtd\ and 
\Vts\, with the current measurements indicating |Vtfc| ^ l^^dl, |Vts|. While a direct 
measurement of the Vtb matrix element is only possible through single top quark 
production, as described in Section !??^ model-dependent constraints on Vtb can also 
be inferred from a measurement of R: Assuming the validity of the standard model, 
specifically the existence of three fermion generations, unitarity of the CKM matrix 
and insignificance of non-W boson decays of the top quark (see Sections 16. 3116. 5p , 
the denominator in Eq. [Ml equals one. R then simplifies to |146p, and is therefore 
strongly constrained to 0.9980 < R< 0.9984 at 90% C.L. by global CKM fits [165]. 

Deviations of R from unity could, for example, be caused by the existence of 
a fourth heavy quark generation, non standard model top quark decays, or non 
standard model background processes. Consequently, precise measurements of R 
probe for physics beyond the standard model, and provide a required ingredient 
for the model- independent direct determination of the \Vtq\ CKM matrix elements 
from electroweak single top production [472]. 

The most precise measurement of R thus far has been performed by DO in 
the lepton -I- jets channel using data corresponding to an integrated luminosity 
of 0.9 fb^^ [473], superseding the previously published measurement based on 
0.2 fb^^ [474]. The tt signal sample composition depends on R in terms of the 
number of b jets present in the sample, as illustrated in Fig. 1441 By comparing the 
event yields with zero, one and two or more b tagged jets, and using a topological 
discriminant to separate ti signal from background in events without b tags, R can 
be extracted simultaneously with the ti production cross section att- This approach 
yields a measurement of att without assuming B{t Wb) = 100%, and exploits 
the different sensitivity of the two measured quantities to systematic uncertainties, 
thereby improving overall precision. 
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Fig. 44: Left: Fractions of events with 0, 1 and >2h tags for ti events with >4 jets 
as a function of R. Right: 68% and 95% C.L. statistical uncertainty contours in the 
R vs. (Jtt plane around the measured point in 0.9 fb^^ of DO data [473]. 



A maximum likelihood fit to the sample composition observed in data gives 



(see Fig. with a correlation of —58% for a top quark mass of 175 GeV/c^, in 
agreement with the standard model prediction. From the measurement 95% C.L. 
limits are extracted, yielding R > 0.79 and |Vffc| > 0.89, the latter being model- 
dependent as mentioned above. R exhibits no significant dependence on the top 
quark mass within ±10 GeV/c^, while cTff varies by =f0.09 pb per ±1 GeV/c^ in the 
same mass range. The total uncertainty on R in this measurement is 9%, dominated 
by the statistical uncertainty of IIq'^qI and the largest systematic uncertainty from 
the b tagging efficiency estimation of to^o47- The cross section measurement yields 
a result similar but not identical to the measurement on the same dataset [262], 
presented in Section 15.1.21 This is due to the assumption of i? = 1 , and slightly 
different event selection in the latter analysis. 

CDF performed the first measurement of R in Run I using both dilepton and 
lepton + jets events on 0.1 fb^^ of data [475], and repeated the analysis in Run II 
on 0.2 fb~^ of data [476], also finding good agreement with the standard model 
expectation. 

All measurements of R performed thus far at the Tevatron are summarized in 




(stat. + syst.) and 

(stat. + syst.) ±0.50 (lumi) pb 



(35) 
(36) 



Fig.m 
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Fig. 45: Summary of the branching ratio R measurements and their total uncer- 
tainties obtained at the Tevatron by CDF [475,476] and DO [473,474]. 



6.3. Search for neutral- current top decays 

Flavor changing neutral interactions of the top quark with a light quark q ~ u,c 
through gauge {Z,j,g) or Higgs (i?") bosons are forbidden at lowest order and 
are suppressed by the GIM mechanism [56] at higher orders in the standard model 
framework. Consequently, the corresponding FCNC top quark decays are expected 
to occur, at most, with branching ratios at 0(10^^^) [419], well out of reach of 
sensitivity of the Tevatron or the LHC. Any observation of such FCNC decays 
would therefore signal physics beyond the standard model. 

Many extensions of the standard model predict the occurrence of FCNC inter- 
actions, affecting both electroweak single top production (see Section [5.8.5p and top 
quark decay. The branching fractions of FCNC top decays can increase by many 
orders of magnitude in such models, as for example in Supersymmetry [477-479], 
additional broken symmetries [480], dynamical EWSB [481,482] including topcolor- 
assisted technicolor [483] or extended Higgs models such as Two Higgs Doublet 
Models [484-487] . Overviews of such models and their impact on top couplings are 
given in Refs. [419,488-490]. 

A search for the top quark FCNC decay t Zq at the Tevatron is considered es- 
pecially interesting due the large top quark mass and very distinct experimental sig- 
nature (see Fig.l46|. It was already suggested in 1989 [491], well before the discovery 
of the top quark. The expected sensitivity for such a branching ratio measurement 
is O(10"2) at the Run H Tevatron and ©(lO""*) at the LHC [492], while the largest 
expected branching fractions from SM extensions reach up to C(10~^) [482] and 
0(10""*) [419], respectively. The best pubhshed hmit before Run II on B{t Zq) 
was obtained at LEP by the L3 Collaboration via a search for single top quark pro- 
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(a) (b) 

Fig. 46: (a) Feynman "penguin" diagram for the FCNC decay i — > Zc with Z 
£~^£~. Including the corresponding diagrams with a d and s quark in the loop, the 
process is nearly cancelled in the standard model, (b) Event signature for top quark 
pairs containing one FCNC t Zq decay and one W — > qq' decay from t Wb, 
resulting in a final state that contains a Z boson and four jets. 



duction, where no significant deviation from the SM background expectation was 
observed, yielding B{t Zq) < 13.7% at 95% C.L. [422]. 

In Run I, the CDF Collaboration performed a search for the FCNC decays 
t Zq and t — > 79 on a dataset with an integrated luminosity of 0.1 fb^"'^ [425]. 
For the t ^ ^q search, a photon is reconstructed as an energy cluster in the electro- 
magnetic calorimeter, either without an associated track or with a single soft track 
(presumably a random overlap) carrying less than 10% of the photon energy point- 
ing to the cluster. Two event signatures are considered, where the W boson from the 
standard-model- like second top decay branches either leptonically into eVe or /ii/^, 
or hadronically into qq' quarks. Consequently, these samples are selected by re- 
quiring either a charged lepton (e or /i), ^t, at least two jets and a photon, or by 
requiring at least four jets and a photon. In both samples, a photon-jet combination 
must yield a mass between 140 and 210 GeV/c^ and the SM-like top decay must 
contain one b tag. 40% of the t ^q acceptance comes from the photon + multijet 
selection, while the lepton -|- photon mode contributes 60%. After all selections, one 
event remains in the leptonic channel and none in the photon -t- multijet channel, 
with an expected background of about 0.5 events mainly from production with 
additional jets in each channel. This translates into a 95% C.L. upper limit on the 
branching fraction of B{t c"f) + B{t uj) < 3.2%. 

In the t Zq search, W qq' decays from the SM-like second top decay are 
considered together with a leptonically decaying Z boson into e+e^ or fj,^ fi~ . Using 
the leptonic W boson decay as well does not substantially increase the acceptance 
and consequently does not improve the limit. The resulting event signature therefore 
contains four jets and two leptons with an invariant mass consistent with that of 
a Z boson, as illustrated in Fig. [46l Since the branching ratio of Z ^ £^£^ is 
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Fig. 47: Expected and observed mass distributions in Z+ > 4 jets events in 
signal samples with and without b tags and in a background-enriched sample used 
to ascertain uncertainties in the background shape and normalization in 1.9 fb""'^ 
of CDF data [429]. The expected FCNC t ~* Zq signal with the observed 95% C.L. 
upper limit on the branching fraction is shown as well. 



small, this search is less sensitive than the t — > 75 one. One Z —t /i^/i^ event 
passes the selection, with an expected background of «0.6 events from Z+ multijet 
and tt production. This corresponds to a 95% C.L. upper limit on the branching 
fraction of B{t cZ) + B{t uZ) < 33%. These measurements can be translated 
into limits on the FCNC couplings and kz at 95% C.L. [426], which are < 
0.176 and k| < 0.533 (see Section [SXS] and Fig. Eg). 

In Run II, the CDF Collaboration has performed a search for the FCNC decay 
t Zq on a dataset with an integrated luminosity of 1.9 fb~^ [429]. Events consis- 
tent with a leptonically decaying Z boson to e+e^ or ii'^ are selected together 
with at least four jets, one of which can be b tagged (see Fig. H5|) . The event selection 
was optimized in a preceding version of this analysis with 1.1 fb~^ of data [493], 
which was a "blind" cut-based counting experiment. By requiring only one well- 
identified lepton for Z reconstruction, while the second lepton can be formed from 
an isolated track, the acceptance is doubled compared to using only fully identified 
leptons. The sensitivity of the search is further increased by dividing the data into 
two subsamples, one b tagged and one not tagged. The best discriminant found 
to separate signal from background is a mass variable that combines the kine- 
matic constraints present in FCNC decays: Two jets in the event have to form a 
W boson, and together with a third jet a top quark, while the Z boson has to 
form a top quark with the fourth jet. Because the event signature does not contain 
neutrinos, the events can be fully reconstructed. The signal fraction in the selected 
dataset is determined via a template fit in signal samples with or without b tags, 
and a background-enriched control sample is used to constrain uncertainties on the 
background shape and normalization (see Fig. I47p . 
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Table 19: Relative signal acceptances and observed 95% C.L. upper limits on the 
branching fractions for several non-SM top quark decay modes as a function of 
assumed top quark mass in a 1.9 fb^^ doubly b tagged lepton + jets dataset [494]. 



Decay 


"^wx/ww [9o] 


Limit [%] 

(175 GeV/c2) 


Limit [%] 
(172.5 GeV/c2) 


Limit [%] 
(170 GeV/c2) 


B{t Zc) 


32 


13 


15 


18 


B{t^gc) 


27 


12 


14 


17 


B{t 7c) 


18 


11 


12 


15 


B{t invisible) 





9 


10 


12 



Since the observed distributions arc consistent with the standard model back- 
ground processes, a 95% C.L. upper limit is extracted on the branching fraction 
B{t cZ) + B{t uZ) of < 3.7%. The expected limit in absence of signal is 5%. 
This is the best limit on B{t — )■ Zq) to date, starting to constrain predictions from 
a dynamic EWSB model [482]. 

6.4. Search for invisible top decays 

Apart from the direct search for t (Z/j) q decays, as described in the previous 
section, one can also perform an indirect search for "invisible" top quark decays by 
comparing the predicted tt production cross section with the observed yield in data. 
In order to be sensitive to novel top decay modes with this method, these decays 
must exhibit a significantly different acceptance from the standard model top quark 
decay. 

Based on a 1.9 fb~^ doubly b tagged lepton + jets dataset, CDF searches for 
deviations of the observed ti production rate from the theoretical prediction [124] 
due to the decays t Zc, t ^ gc, t ^ ■^c and t — > "invisible" states [494]. These 
decays exhibit a relative acceptance T^wx/wwj where one novel and one SM top 
quark decay occur, normalized to the standard model tt decay acceptance, from 
32% down to no acceptance (for decays to invisible states). 

With an observed tt production cross section of 8.8 pb, and a prediction of 
6.7 pb for a top quark mass of 175 GeV/c^, the obtained limits on the novel top 
quark decay modes are all lower than expected, but statistically consistent with 
expectation. The results are summarized in Table [T9l for top quark masses of 170, 
172.5 and 175 GeV/c^. 

6.5. Search for top decays to charged Higgs bosons 

As indicated previously, the standard model incorporates one Higgs doublet of com- 
plex scalar fields to break electroweak symmetry and to generate masses of weak 
gauge bosons and fermions (see Section l2.ip . As a consequence, one obtains a sin- 
gle neutral scalar CP-even particle that still remains to be discovered, the Higgs 
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boson H. An extension of the standard model Higgs sector introduces a second 
Higgs doublet, referred to in Section [5.8.41 as Two Higgs Doublet Models (THDM 
or 2HDM) [414,415]. These models provide five physical scalar Higgs bosons af- 
ter electroweak symmetry breaking, namely two neutral CP-even Higgs bosons H'^ 
and /i", one neutral pseudoscalar CP-odd Higgs particle A", and two charged Higgs 
bosons H^. The observation of charged Higgs bosons would therefore offer clear 
evidence for physics beyond the standard model. 

Three choices of Higgs-fermion couplings are differentiated in THDM. In Type-I 
models only one of the two Higgs doublets couples to fermions, while in Type-II 
models one doublet couples to up-type fermions and the other doublet to down- 
type fermions. Type-III models have general Higgs-fermion Yukawa couplings of 
both Higgs doublets, leading to Higgs-mediated FCNC at tree level, which requires 
tuning of the Higgs parameters to ensure sufficient suppression of FCNC to be 
compatible with current experimental limits. One example of a Type-H THDM is 
the Minimal Supersymmetric Standard Model (MSSM) [369] which is frequently 
used as reference in the analyses described below. The relevant model parameters 
in searches for charged Higgs bosons are the ratio of the vacuum expectation values 
of the two Higgs doublets (tan (3) and the mass of the charged Higgs boson m[j± . 

The inclusive single charged Higgs boson production rate <t{pp — > tH~X) 
reaches a maximum of 0(1 pb) at the Tevatron where the charged Higgs boson 
can be produced via the decay of a top quark. The corresponding inclusive cross 
section for pair production of charged Higgs bosons <j{pp H^H^X) is below 
0(0.1 pb) [495,496]. 

The decay mode t Hh is kinematically accessible if the mass of the charged 
Higgs boson is less than the difference between top and b quark masses mu± < {mt — 
mb), and will then compete with the standard model decay t — > Wb. The distinct 
top quark decay signature provides an additional handle for background suppression 
compared to direct production of charged Higgs bosons. The branching fraction of 
t Hh depends on tan (3 and mu± . As illustrated in Fig. [JH] for the MSSM, the 
branching ratio of i — > Hb increases significantly both for small tan /3 ^ 1 and large 
tan/3 ^ 40, for a given mfj±. The standard model decay is assumed to account for 
the difference of B{t Hh) from unity. For a given tan/3, the branching ratio of 
t Hb decreases with increasing m//± . The decay of H^ is dominated by H^ tv 
for large tan/3, independent of rnui, which would result in an excess of ti events 
in the r decay channel relative to standard model expectation. At small tan/3, the 
decay H^ — > cs is enhanced for small TO/f±, while H^ — )■ t*b dominates for m}j± 
close to the top quark mass. Consequently, searches for charged Higgs bosons focus 
on these three fermionic decay modes. 

In the early 1990s (before the discovery of the top quark) first searches for 
t — > Hb in the H^ — > tv decay mode assuming specific branching fractions were 
performed, and limits derived in the mt versus mfj± parameter space by the UAl 
and UA2 experiments [498, 499] at the CERN SppS colhder, and by CDF at the 
Run I Tevatron [500,501]. All four LEP experiments searched for pair-production of 
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Fig. 48: Branching ratios for top quark and charged Higgs boson decays versus tan/? 
in the MSSM framework as simulated with CPsuperH [497] for m^± = 100 GeV/c^ 
(left) and mjj± = 150 GeV/c^ (right). 



charged Higgs bosons in e+e^ H^H^ , assuming that only the decays rv 
and cs can occur [502-505], as favored by Type-II THDM. The dominant 

background in the resulting three decay modes is pair production of W bosons, yield- 
ing similar final states. 95% C.L. lower mass limits, independent of the decay 
mode, yielded m//± > 78.6 GeV/c^ in a preliminary combination of all four experi- 
ments [506] . This was superseded by a more stringent limit obtained by ALEPH of 
79.3 GeV/c2 [502]. 

Indirect limits on the mass of the charged Higgs boson can be obtained from 
measurements of the 6 — > 57 FCNC process at B factories, since the involved loop 
diagrams are sensitive to contributions from new particles such as . For a Type-H 
THDM scenario, a 95% C.L. lower limit of mjy± > 295 GeV/c^ can be derived [507] 
if the used theoretical description is assumed to be complete. Direct searches are 
less model dependent, and therefore serve as important tools to scan for new physics 
beyond the regions of parameter space excluded by the corresponding direct analyses 
described above. The direct searches for t Hb performed at the Tevatron are 
based on Type-II THDM scenarios. 

After the discovery of the top quark, the first searches for H"^ in top decays tt 
H^W^bb,H^H^bb focused on the decay H"^ ri/, which corresponds to large 
tan/3. CDF published an analysis superseding and extending a previous result [508], 
requiring inclusive final states with ^t, a hadronically decaying tau lepton (t/j), 
and (i) two jets and at least one additional either lepton or jet, or (ii) a second 
energetic [509]. Another search investigated the dilepton channels eTh,fJ.Th, with 
accompanying and at least two jets [510]. DO performed a first analysis in 
Run I based on a "disappearance" search in the lepton -I- jets channel, sensitive to 
fermionic decay modes, by looking for a discrepancy in the event yields relative 
to the standard model predictions [511]. This analysis was complemented through a 
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direct search for tv with a hadronically decaying r reconstructed as narrow 

jet in a dataset with events containing and at least four, but no more than eight, 
jets [512]. All analyses observe good agreement with SM expectation, and provide 
limits in the tan/3, m//± plane. 

It should be noted that the above limits are based on tree level MSSM cal- 
culations of branching fractions depending on tan (3. By now, it has become clear 
that higher-order radiative corrections, which strongly depend on model parame- 
ters, modify these predictions significantly [513,514]. Also, non-fermionic (bosonic) 

decay modes can have non-negligible contributions at small tan /?, as illustrated 
in Fig.[48l affecting the limits derived in that area without taking this into account. 
Independent of these issues, one can still provide upper limits on B{t Hh) based 
on the observed production rate for any specified branching ratio. For example, 
for a purely tauonically decaying charged Higgs boson, 95% C.L. upper limits on 
B{t Hb) are found by CDF to lie between 0.5 and 0.6 for 60 GeV/c^ < m„± < 
160 GeV/c2 [510]. The combined DO result corresponds to B{t Hb) < 0.36 at 
95% C.L. for mH± < 160 GeV/c^ and 0.3 < tan/3 < 150, which is the fuU range 
where the leading-order MSSM calculation is valid [512]. 

CDF published a first search for t Hb in Run II using 0.2 fb^^ of integrated 
luminosity [515]. The search is based on the corresponding ti cross section analy- 
ses [207,216,252] in the topology + jets + £ + X, where £ corresponds to an 
electron or muon and X to either £ (dilepton channel), Th (lepton + r channel) or 
one or more b tagged jets (lepton -I- jets channels). Dropping the assumption of 
B{t Hb) = 0, and avoiding overlaps of the channels, the observed yields can be 
compared with the expected deficits or excesses in the channels relative to the stan- 
dard model prediction, depending on the top quark and H^ branching fractions. 
Apart from the standard model top quark decays, t — > Hb is considered with H^ 
decaying to tv, cs, t*b or Wh°, with h° — > bb. The ti production rate is assumed to 
be not affected by the extension of the Higgs sector. Since no iJ^ signal is observed, 
95% C.L. upper limits on B{t Hb) are obtained, for example, for a tauonic Higgs 
model ^ Tv) = 1) to be 0.4 for 80 GeV/c^ < mjy± < 160 GeV/c^. 95% 

C.L. limits are also obtained in the {m}j± , tan/3) parameter space in the framework 
of the MSSM for certain benchmark settings of parameters [515], taking radiative 
corrections into account. While the excluded area for large tan/3 strongly depends 
on the different benchmarks investigated, this is not the case for small tan/3. An 
example result is shown in Fig. |49l 

DO performs a similar analysis based on the same tt final states in 1 fb^^ of 
integrated luminosity [257]. Two models for the decay mode of the charged Higgs 
boson are studied: (i) a "tauonic" Higgs model with B{H^ tv) — 1 that would 
give an enhancement of the lepton + r channels and a deficit in the lepton -I- jets 
and dilepton channels, and (ii) a "leptophobic" model with B{H^ cs) = 1 that 
would yield an enhancement of the all-hadronic channel and a deficit in all channels 
considered in this analysis. For both model assumptions, good agreement with the 
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Fig. 49: Left: Expected and observed 95% C.L. exclusion limits for charged Higgs 
bosons in the {mu± , tan (3) plane obtained by CDF for the MSSM benchmark sce- 
nario discussed in Ref. [515]. Right: Expected and observed 95% C.L. limits on 
B{t Hb) versus mff± found by DO for a tauonic charged Higgs model using a 
simultaneous fit of B{t Hb) and the tt production cross section [257]. MSSM 
tree-level predictions for several tan /3 values are shown as well. 



standard model prediction is observed, and 95% C.L. upper limits on B{t — > Hb) 
are provided for 80 GeV/c^ < mH± < 155 GeV/c^, yielding 0.16 - 0.20 for the 
tauonic and 0.2 for the leptophobic model. The dominant systematic uncertainties 
arise here from uncertainties on the tt cross section and the luminosity. 

For the tauonic model, an improvement of the obtained limits by about 30% 
in the low mH± range (< 100 GeV/c^) is possible when the ti cross section is 
allowed to float in the fit rather than be fixed to the SM value. The resulting limits 
are displayed in Fig. 2^1 and range from 0.12 to 0.26 in the indicated m^^i range. 
Assuming the standard model scenario of B{t Hb) = 0, a combination of the 
ti cross sections from the analyzed final states is also obtained, as was discussed in 
Section 15.1.51 

As illustrated in Fig. [49l the results obtained in Ref. [515] leave room for im- 
provement, particularly close to tan/3 w 1 and mjy± above the W boson mass. For 
this range of parameters, the MSSM predicts a significant branching fraction for the 
decay H^ cs. CDF has searched for t Hb in this decay channel in 2.2 fb~^ dou- 
ble b tagged lepton + jets data [516]. Both standard model and exotic decay exhibit 
the same final state, but can be distinguished via the dijet invariant mass, where the 
two untagged leading jets are assigned to W^/H^ decay products. A binned likeli- 
hood fit using W^/H^ dijet mass templates yields no significant excess over the SM 
prediction, and 95% C.L. upper limits on B{t Hb) are provided for 90 GeV/c^ 
< < 150 GeV/c^, assuming a leptophobic Higgs model. The limits range from 

8% for mH± = 130 GeV/c^ to 32% for mH± = 90 GeV/c^, complementing the 
analysis by DO described above for the mass range above 100 GeV/c^. 
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A first direct search for charged Higgs boson production in the mass range 
beyond rrit via the process qq' — > tb has been performed by DO, as was 

discussed previously in Section 15.8.41 

7. Top Quark Properties 

In the past two chapters, it was demonstrated that both top quark production 
and decay thus far have been found to be consistent with the standard model 
expectations. No new particles or anomalous couplings have been observed yet. To 
actually confirm the top quark's standard model identity, its fundamental quantum 
numbers need to be measured and their self-consistency in the standard model 
framework needs to be confirmed as well. 

In this chapter, measurements of the top quark's electric charge, lifetime and 
mass performed thus far at the Tevatron are described. Again, top quark pair events 
are used for this, because these events provide higher statistics and favorable sample 
purities compared to single top events. First direct measurements of the Vtb matrix 
element in electroweak single top quark production have already been discussed in 
Section [nXB 

7.1. Top quark electric charge 

The electric charge of quarks can be determined, for example, in electron-positron 
collisions via the ratio of the hadronic cross section to the muon cross section R = 
(j{e^e~ — > hadrons)/(T(e~''e~ — > fi^fj,~), which is proportional to the sum of the 
squared electric charges of the quark flavors accessible at the chosen center of mass 
energy. Due to the large top quark mass, such a direct measurement could not 
yet be performed at past and present electron-positron colliders. Also, a direct 
measurement of photon radiation in ti events at hadron colliders, which would give 
access to the top quark's charge and its electromagnetic coupling, is unrealistic due 
to limited statistics at the Tevatron [517]. Consequently, the top quark is the only 
quark whose fundamental quantum numbers of weak isospin and electric charge 
could thus far be determined only indirectly in the framework of the standard 
model from measurements of its weak isospin partner, the b quark, to be Ta = -|-^, 
Qj = -|-|e (see Section [2?2t . 

Information on the electric charge of the top quark can also be inferred from 
the electric charges of its decay products. However, there is an inherent ambiguity 
in — > — > W^W~bb events when pairing W bosons and b jets, resulting in 
possible charges of |Q|= 2e/3 or 4e/3 for the top quark. An exotic quark with 
charge — 4e/3 being the discovered particle at the Tevatron instead of the standard 
model top quark would be compatible with precision electroweak measurements 
if the right-handed b quark were to mix with the — le/3 charged exotic doublet 
partner of such an exotic top quark. The standard model top quark with charge 
2e/3 would yet have to be discovered in this scenario, due to its large mass of 
271^11 GeV/c2 [366,518,519]. 
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Fig. 50: Left: b/b jet charge distributions obtained in dijet data. Right: Distribu- 
tion of the top quark charge obtained in data, overlaid with the expectations from 
standard model and exotic model [520]. 

DO has published a first measurement discriminating between the 2e/3 and 4e/3 
top quark charge scenarios in a 0.4 fb~^ lepton + jets dataset with > 2 b tagged 
jets [520]. The sample exhibits high signal purity, with the two dominant back- 
ground processes of Wbb and single top production contributing only 5% and 1% 
to the selected events, respectively. Each tt event provides two measurements of 
the absolute value of the top quark charge, one from the leptonic and one from the 
hadronic decay of the W boson in t ^ Wb. The charge of the W boson is determined 
from the (inverse) lepton charge for the leptonic (hadronic) W boson decay. The 
b jet charge, discriminating between b and b jets, is determined using a jet charge 
algorithm, based on the pr weighted average of the charges of the tracks associated 
with the b tagged jet. The corresponding distributions are derived from dijet col- 
lider data [520] (see Fig. [501). The top quark charge observable is then defined as the 
absolute value of the sum of the W boson- and associated b jet charge, where the 
right pairing is determined through a constrained kinematic fit. By comparing the 
distribution in data with the expected shape from the standard model and exotic 
model (see Fig. [50)) . respectively, DO excludes the hypothesis of only exotic quarks 
of charge \Q\= 4e/3 being produced at up to 92% C.L. and limits an exotic quark 
admixture in the sample to at most 80% at the 90% C.L. 

Using a similar analysis approach, CDF obtains a preliminary result on the top 
quark charge using double b tagged lepton -I- jets and b tagged dilepton events in 
a 1.5 fb~^ dataset [521]. The observed 21n(Bayes Factor) is 12, meaning that the 
data favor very strongly the SM top quark hypothesis over the exotic model. 

While the results from CDF and DO are not directly comparable because of 
different statistical approaches in the interpretation of their results, both agree to 
favor the SM top quark charge hypothesis. This is supported by the searches for 
new heavy top- like quark pair production (see Section I5.7.4|) , which are starting to 
exclude additional quark production in the mass range predicted by the alternative 
model. A top quark charge measurement determining the b jet charge in soft-lepton 
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tagged events from the soft-lepton charge rather than the current track-based ap- 
proach has not been performed yet. 

The measurement of the top quark charge in tt-f events will be possible at the 
LHC because of the increased production rate and the reduction of irreducible 
background from photons radiated off the incoming quarks in it production via qq 
annihilation. Consequently, top quark charge measurements via photon radiation in 
ti events are predicted to achieve a precision of 10% [517] . Using the top quark decay 
products to provide an additional charge measurement will help to disentangle the 
measurements of top quark electromagnetic coupling strength from the top quark 
charge in the ttj events. This will help to rule out possible anomalous admixtures 
in the electromagnetic interaction of the top quark. 

7.2. Top quark lifetime 

The lifetime of the top quark Tt — hT^^ « h{1.3GeV)~^ is approximately 
5 • 10^^^ s in the framework of the standard model, as discussed in Section 13.3.21 
Consequently, the production and decay vertices of the top quark are separated 
by C'(10~^^)m, which is orders of magnitude below the spatial resolution of any 
detector. Also, the top quark width is narrower than the experimental resolution 
at both the Tcvatron and the LHC. Consequently, a direct measurement of the top 
quark lifetime or its width will be limited by detector resolution. A measurement 
of the top quark lifetime is still useful to confirm the standard model nature of the 
top quark, and exclude new top quark production channels through long-lived par- 
ticles. A measurable lifetime of the top quark itself would imply a correspondingly 
small Vtb matrix element and render single top quark production at the Tevatron 
undetectable, in contradiction with the observed evidence described in Section [5.81 

CDF has set limits on the top quark lifetime and width using two approaches. 
One analysis uses a b tagged lepton -I- jets dataset of 0.3 fb^^ [522] to measure the 
impact parameter (smallest distance) between the top quark production vertex and 
the lepton track from the leptonic W boson decay in the plane orthogonal to the 
beam direction. Fitting the obtained distribution with signal Monte Carlo templates 
for cTf between and 500 /im and a background template, the template for /xm 
describes the data best, which translates into a 95% C.L. limit on cTt of < 52.5 /im. 

The second analysis uses a b tagged lepton -I- jets dataset of 1 fb^^ [523] to 
reconstruct the top quark mass in each event using a kinematic fit. The observed 
distribution is compared in a fit to ti signal Monte Carlo templates of different 
widths for a top quark mass of 175 GeV/c^ and background templates. From the 
fit result, at 95% C.L. the top quark width is found to be smaller than 13.1 GeV, 
corresponding to a lower limit on the top quark lifetime of 5 • 10~^^ s. 

7.3. Top quark mass 

The top quark is set apart from all other known fundamental particles by its large 
mass. Being the only particle with its Yukawa coupling close to unity also raises the 
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question whether it plays a special role in the process of mass generation. Since the 
lifetime of the top quark is so short (see Sections 13.3.21 and 17. 2p . unlike the other 
quarks, it does not hadronize and properties like its mass can be determined directly 
without the complication of having a quark embedded in a hadron. Being a sensitive 
probe for physics beyond the standard model, it is also important to measure its 
mass in different decay modes that could be affected differently by novel physics 
contributions. 

In the framework of the standard model, the top quark mass is a free parameter. 
As discussed in Section I2.3i its precise determination together with a precise W 
boson mass measurement provides a test of the self-consistency of the framework and 
constrains the mass of the yet undiscovered Higgs boson (or other new particles [524, 
525]) via electroweak radiative corrections. 

Measurements of the top quark mass have been performed thus far only in 
ti events, and mainly in the dilepton, lepton + jets and all-hadronic final states. 
A complete kinematic reconstruction of the tt pair from the observed objects in 
the event can be performed in the all-hadronic final state where no neutrinos are 
present. Assuming arises solely from the escaping neutrino in the lepton -I- jets 
channel, a kinematic fit can be performed here constraining the invariant mass of the 
charged lepton and neutrino to that of the W boson, yielding a twofold ambiguity 
for the neutrino's longitudinal momentum solution. Because of the two neutrinos 
contributing to in the dilepton final state, a direct kinematic reconstruction 
of the ti event is not possible without adding more information or making other 
assumptions about the kinematics of the objects in the event. 

Since the assignment of partons to reconstructed objects in an event is not 
definite, combinatorial ambiguities arise in all channels. Depending on the analysis 
technique, either all combinations are used to extract the top quark mass, or the 
best combinations are selected based on, for example, the lowest in ^ kinematic fit 
to the tt event hypothesis. Identification of b jets can be used to reduce the number 
of combinations to consider in the lepton + jets and the all-hadronic channels. 
Even if both b jets from the tt decay are identified, four combinations remain in 
the lepton -I- jets channel (including the neutrino pz ambiguity) and six in the all- 
hadronic channel. In these channels, usually at least one b tagged jet is required to 
increase sample purity, or the data can be split on the basis of b tag multiplicity, 
and therefore purity, to optimize overall sensitivity. 

The techniques used in top quark mass analyses can be divided into three cate- 
gories: 

(i) Template Method (TM): Observables sensitive to the mass of the top quark 
such as the reconstructed top quark mass mj-cco or Ht are evaluated in the 
dataset under consideration. The resulting distribution is then compared in a 
fit with expected contributions from ti signal (with varying top quark masses) 
and background processes. 

(ii) Matrix Element Method (ME): Based on the leading-order matrix ele- 
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ments of contributing signal and background processes, the four-vectors of the 
reconstructed objects in each event define a probabiHty density as a function 
of the top quark mass. The total likelihood for the event sample is given as the 
product of the individual event likelihoods. This method is also referred to as 
the Dynamical Likelihood Method (DLM). 
(iii) Ideogram Method (ID): An event-by-event likelihood depending on the as- 
sumed top quark mass is formed based on a constrained fit of the event kine- 
matics, taking all object permutations and possible background contributions 
into account. As mentioned above, this kind of constrained fit is only possible 
in the all-hadronic and lepton -I- jets channels. 

Naturally, the analyses most sensitive to the top quark mass are also very sensitive 
to the jet energy scale (JES) calibration. The systematic uncertainty due to the 
external jet energy calibration (see Section [4.3.3|) then is usually the dominant sys- 
tematic uncertainty in such analyses. This can be reduced in decay channels where 
at least one W boson decays "hadronically" by using the well-measured W boson 
mass to constrain in-situ the jet energy calibration in top quark events [526,527]. 
Determining such an overall scale factor for jet energies, absorbs a large part of 
this uncertainty into an uncertainty scaling with ti sample statistics, while residual 
uncertainties, for example, due to rj and pr dependence of JES corrections or differ- 
ences between light-quark and 6-quark JES still remain. By performing an analysis 
simultaneously in the dilepton and all-hadronic and/or lepton -I- jets channels, the 
in-situ JES calibration can also be transferred to the dilepton channel [528,529]. 

Another approach to reduce the dependence of top quark mass measurements on 
JES is to utilize observables with minimal JES dependence that are still correlated 
with the top quark mass, such as the mean px of the charged lepton from the W 
boson decay (pfp) or the mean transverse decay length of 6 jets L^y in ti events [530]. 
While such measurements are statistically limited at the Tevatron, their uncertain- 
ties are basically uncorrelated with those of other statistically more sensitive anal- 
yses. This helps to reduce the overall uncertainty on mt when all measurements are 
combined. Also, the observed signal event yield can provide an additional constraint 
on mt via the mass dependence of the tt production cross section [531]. 

Performance, calibration and statistical uncertainty derivation of each mass anal- 
ysis are checked using sets of simulated pseudo-experiments (ensemble tests) , based 
on mean and rms of the extracted mass and pull distributions. 

Measurements of mt were pioneered in Run I, based on 0.1 fb^^ of data 
in the dilepton [316,532-535], lepton + jets [320,321,454,536,537] and all- 
hadronic [538, 539] channels. A combination of the Run I results yields mt — 
178.0 ± 2.7(stat.) ± 3.3(syst.) GeV/c^ [158]. Also in Rim II, results come mainly 
from the dilepton, lepton -I- jets and all-hadronic channels, with the most precise 
measurements being from lepton -I- jets samples. One analysis uses an inclusive + 
jets signature, vetoing energetic isolated leptons and thereby enhancing the r -I- jets 
signal contribution to 44% [540] . This result is listed together with measurements in 
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the all-hadronic channel in Section [7. 3. 31 and is consistent with the world- averaged 
top quark mass. A top quark mass measurement with explicit hadronic r recon- 
struction has not been performed thus far, but given the recent progress in the 
corresponding cross section analyses discussed in Section I5.1.4[ this could still be 
feasible at the Tevatron. 

In the following sections, the latest (final) Run I results are given, along with the 
current preliminary and published Run II top quark mass measurements using vari- 
ous analysis techniques for each of the three main decay channels. Some of the most 
precise analyses entering the world-averaged top quark mass will be highlighted. A 
more detailed review of top quark mass analysis techniques pursued at the Tevatron 
can be found in Ref. [541]. The final section presents the current world average and 
some of its implications. 

7.3.1. Dilepton final state 

Analyses in the dilepton final state are performed either based on the matrix ele- 
ment or the template method. For the template approach, additional assumptions 
on the kinematics of the involved objects are made in order to solve the otherwise 
underconstrained system kinematics. Assuming several top quark masses, the con- 
sistency of the observed event kinematics can be used to obtain weights for each 
event as a function of mi, based on input parton distribution functions and the ob- 
served charged lepton energies ( "matrix weighting" M), or using simulated neutrino 
pseudo-rapidity or azimuthal angle ("neutrino weighting": v^i, v^,), or ti longitudi- 
nal momentum (p**) distributions [534,542]. Top quark mass estimators are derived 
from the obtained weight distributions, such as the peak mass position, or mean 
and rms of the distributions. These values are then used in a template fit to data 
to obtain the most likely mt from the sample. 

The most precise top quark mass result in the dilepton channel entering the 
world average is from DO, and has a precision of 2.2% [543]. It combines results 
from neutrino weighting (i/^) obtained on 1 fb~^ in the dielectron, dimuon and 
lepton -I- track channels [544] with a measurement in the e/x channel using 2.8 fb^^ 
of data and the matrix element method [543]. 

The matrix element method evaluates the probability density for each event 
(Pevt) with measured object four-vectors x to originate from ti production, depend- 
ing on the top quark mass, or from the dominant background arising from Z ^ tt 
+ jets production in the following linear combination, based on the known expected 
signal fraction in the sample /: 

Pcvt(2:;"it) = / • Psig{x;mt) + (!-/)• -Pbkg(a;). (37) 

Psig and Pbkg are the signal and background probability densities for ti and Z ^ tt 
+ jets production, based on the leading order matrix element for qq — > ti and the 
VECBOS [555] parametrization of the Z-production matrix element, respectively. 
The probability densities are calculated by integrating over all unknown quantities. 
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Tabic 20: Top quark mass measurements performed thus far at the Tcvatron in the 
dilepton channel with their integrated luminosities, data selections {£i = dilepton, 
^+trk = lepton + track, NN = neural network) and analysis methods used. The two 
analyses using 0.1 fb~^ are from Run I, and the references marked with an asterisk 

correspond to preliminary results. 



JCdt 
[fb-i] 


Selection 


Method 


rrit ± (stat.) ± (syst.) 
[GeV/c2] 


Ref. 


0.1 


a 


TM:i/,, 


167.4 ± 10.3 ±4.8 


[532,533] 


0.1 


u 


TM:7W,i/^ 


168.4 ±12.3 ±3.7 


[535] 


0.3 


u 


ME 


165.2 ±6.1 ±3.4 


[545,546] 


0.4 




TM:i.^,z/^,p** 


170.1 ±6.0 ±4.1 


[542] 


0.4 




TM:A^,f^ 


178.1 ±6.7 ±4.8 


[547] 


1.0 


u 


ME 


164.5 ±3.9 ±3.9 


[548] 


1.0 


a 


TM:M 


175.2 ±6.1 ±3.4 


[549]* 


1.0 






176.0 ±5.3 ±2.0 


[544]* 


1.2 


u 


TM:p** 


169.7tt:9±3.1 


[531] 


1.2 


u 


TM:p**®atf 


170.7t;^:^±2.6±2.4(th.) 


[531] 


1.8 


11 


TM:p^ 


156119 ± 4.6 


[550]* 


1.8 


a 


TU.:Vn 


n2.Qtll ± 3.6 


[551]* 


1.8 


a 


ME 


170.4 ±3.1 ±3.0 


[552]* 


2.0 


«e NN 


ME 


171.2 ±2.7 ±2.9 


[553] 


2.8 


£+trk 




165.lt|J±3.1 


[554]* 


2.8 


ejj, 


ME 


172.9 ±3.6 ±2.3 


[543]* 


2.8 


e/i/££,£+trk 


ME/TM:i/^ 


174.4 ±3.2 ±2.1 


[543]* 



such as the unmeasured neutrino energies and all parton states that can lead to the 
X objects observed in the detector: 

Psig(a;;mt) = l/CTobs(mt) / V dqidq2f{qi)f{q2) ^'^^^ ^"^^ d^6W{x,y), 

Jqi,g2,ya 91Q2S 
y± i-j^ 5y flavors 

(38) 

where qi and 52 are the momentum fractions of the colliding partons from the 
proton and antiproton, f{qi) the corresponding PDFs, M is the matrix element 
for the signal process yielding the partonic final state y, s is the squared center- 
of-mass energy and rf$f, a six-body phase space clement. The transfer function 
W{x,y), which incorporates the detector resolution, describes the probability for 
a final state x in the detector to be reconstructed from the partonic state y. The 
two possible permutations from the unknown jet-parton assignment are summed 
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over, and the probability is normalized to the expected observable production rate 
o'obsiiTi't)- The calculation of Pbkg(a^) proceeds in an analogous way, except that 
there is no dependence on rrit. 

The top quark mass of an event sample can be obtained by maximizing the 
total likelihood function, which is the product of the individual event probabilities, 
with respect to mt . In this way, each event contributes according to its quality and 
inherent resolution. While the ME technique exploits the full kinematic information 
available, and usually yields the statistically most sensitive measurements, it is also 
computationally intensive because of the involved multidimensional integrations. 
The result of the 2.8 fb~^ analysis is given together with other mass measure- 
ments performed in the dilepton channel at the Tevatron in Table [20l The dominant 
systematic uncertainty in this analysis arises, as expected, from systematic uncer- 
tainties on the JES calibration, and response differences between light quarks and 
b quarks. 

7.3.2. Lepton + jets final state 

The precision of the world-averaged top quark mass is driven by measurements in 
the lepton -f- jets channel that provides the best compromise between sample purity 
and signal statistics. In this channel, all three analysis methods (TM, ME, ID) 
have been deployed, with the most precise results consistently coming from matrix 
element analyses, starting with its first application at DO in Run I [454]. 

For the current world-averaged mass, CDF and DO contribute measurements 
in the lepton -t- jets channel based on the matrix element method, as discussed 
in the previous section, simultaneously fitting the top quark mass and an overall 
in-situ JES scale factor to the data. Events with one energetic isolated lepton, large 

and exactly four jets are selected so as to best comply with the leading order 
matrix element used in the calculations. Also, at least one of the jets is required 
to be b tagged. Using datasets of 2.7 and 2.2 fb-\ for CDF [556] and DO [557], 
respectively, both measure the top quark mass with a precision of 1.0%, incidentally 
even yielding the same mass value of 172.2 GeV/c^. A comparison of the invariant 
dijet and three-jet mass distributions based on the permutation of largest weight in 
data and simulation is shown in Fig. [51] for a 1 fb^^ subset [558] of DO's 2.2 fb^^ 
ME analysis. 

The largest systematic uncertainties on the measurement by CDF arise from 
the tt MC generator used to calibrate the result (pythia versus herwig) and the 
residual JES uncertainty. For DO, the dominant uncertainty comes from the 6-jet 
over light -jet calorimeter-response ratio and the signal modeling uncertainty, taking 
the impact of extra jets into account based on the observed four to at least five 
jet event ratio in data. Both experiments are in the process of streamlining their 
methods used to assess systematic uncertainties as well as exploring new sources of 
uncertainties that start to become important at the current level of precision [567]. 
Examples for the latter are differences arising from using NLO rather than LO MC 
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m2j(GeV) m3j(GeV) 

Fig. 51: Dijet (left) and three-jet (right) invariant mass distributions in 1 fb^^ of 
DO data compared with simulated ti signal (for mt — 170 GeV/c^) and background 
contributions [558]. The results are for jet permutations of largest weight using the 
ME method. 



generators, or non-perturbative QCD effects such as color reconnection [184, 185]. 

Table [21] summarizes the latest (final) Run I results and the current preliminary 
or published Run II top quark mass measurements in the lepton + jets channel. 

7.3.3. All-hadronic final state 

Analyses in the all-hadronic channel have been performed so far using template and 
ideogram methods, and exhibit comparable sensitivity. The precision in this channel 
by now is similar to that in the dilepton final state, mainly as a result of using 
in-situ JES calibration, thereby reducing the otherwise overwhelming systematic 
uncertainty from the external JES. 

CDF's best measurement in the all-hadronic channel entering the world average 
has been performed using a template method on 2.1 fb^^ of data, yielding a precision 
of 2.4% [568]. Similar to the analysis [278] described in Section [5.1.3[ events are 
required to have between six and eight energetic central jets, no isolated energetic 
lepton or significant ^t, and have to pass a selection based on the output of a 
neural network discriminant. The dataset is split into subsamples with exactly one 
and exactly two b tagged jets. 

The leading six jets in each event define the tt signal and multijet background 
templates for the reconstructed mass of the top quark and of the W boson, with 
assignments based on a kinematic fit, where the permutation with lowest is 
selected for further study. While the signal templates depend on both rrit and the 
JES scale factor, the background templates are assumed to not depend on the top 
quark mass, and no JES dependence is considered either. The measurement is then 
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Table 21: Top quark mass measurements performed thus far at the Tevatron in the 
lepton + jets channel with their integrated luminosities, data selections (£+jets = 
lepton + jets, ££ — dilepton) and analysis methods used. The two analyses using 
0.1 fb~^ are from Run I, and the references marked with an asterisk correspond 
to preliminary results. Measurements marked with a cross contain the uncertainty 
from the in-situ JES calibration within the quoted statistical uncertainty. 



J Cdt 

[fb-i] 


Selection 


Method 


rrit ± (stat.) ± (syst.) 
[GeV/c2] 


Ref. 


0.1 


£+jets 


TM:mroco 


176.1 ±5.1 ±5.3 


[537] 


0.1 


£+jets 


ME 


180.1 ±3.6 ±3.9 


[454] 


0.3 


£+jets 


DLM 




[526,559] 


0.3 


£+jets 


TM:mrcco ® JES 


173.5li^-^± 1.3t 


[526,527] 


0.4 


£+jets 


ME © JES 


170.3ttUli^ 


[560] 


0.4 


£+jets 


ME ® JES 


m.2t',itlV 


[560] 


0.4 


£+jets 


ID ® JES 


173.7 ±4.4t|Jt 


[561] 


0.7 


£+jets 




180.7t^^;^±8.6 


[562] 


1.0 


£+jets 


ME e JES 


171.5± 1.8± 1.1^ 


[558] 


1.0 


£+jets 


ME ® JES 


170.8 ± 2.2 ±1.4t 


[563] 


1.0 


£+jets 


TMimioco 


168.9 ±2.2 ±4.2 


[564]* 


1.2 


£+jets 


ME e JES 


173.0 ± 1.9 ± l.Ot 


[557]* 


1.7 


£+jets 


DLM © JES 


171.6 ± 2.0 ±1.3^ 


[565]* 


1.9 


£+jets 


TM:L,„,p^ 


175.3 ±6.2 ±3.0 


[530]* 


1.9 


£+jets 

U 


TM:mrcco © JES 
TM-.i^,,, Ht 


171.9± 1.7± l.lt 


[528,529] 


1.9 


£+jets 


ME ® JES 


172.7± 1.8± 1.2t 


[566] 


2.2 


£+jets 


ME ® JES 


172.2 ± 1.0 ± 1.4 


[557]* 


2.7 


£+jets 


ME © JES 


172.2 ± 1.3 ± l.Ot 


[556]* 



performed in a two-dimensional fit to mt and in-situ JES scale factor, using these 
templates and the observed distribution in data in both subsamples, and applying 
a Gaussian constraint from the external JES calibration. The reconstructed mass 
distributions obtained in data with two b tags, overlaid with templates for expected 
background and tt signal (for nit — 177 GeV/c^ and unchanged JES with respect 
to the external calibration, as obtained from the fit) are shown in Fig. [52l 

The result of this analysis is given in Table [22l together with other top quark 
mass measurements performed in the all-hadronic channel at the Tevatron. The 
dominant systematic uncertainties in this analysis arise from uncertainties on shape 
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Fig. 52: Reconstructed mass of the W boson (left) and top quark (right) obtained 
in 2.1 fb^^ of double b tagged CDF data compared to the expected contributions 
from multijet background and 177 GeV/c^ ti signal, using the JES corresponding 
to the best fit [568]. 



Table 22: Top quark mass measurements performed thus far at the Tevatron in the 
all-hadronic channel with their integrated luminosities, data selections and analysis 
methods used. The two analyses using 0.1 fb~^ are from Run I, and the references 
marked with an asterisk correspond to preliminary measurements. Results marked 
with a cross contain the uncertainty from the in-situ JES calibration within the 
quoted statistical uncertainty. 



JCdt 

[fb-i] 


Selection 


Method 


nit ± (stat.) ± (syst.) 
[GeV/c2] 


Ref. 


0.1 
0.1 
0.3 
0.3 
0.9 
1.0 
1.9 
2.1 


jets only 
jets only 
^T+jets 
jets only 
jets only 
jets only 
jets only 
jets only 


TM:mroco 
TM:mroco 
TM:Ht 
ID 

TM:ME e JES 

TM:mi.oco 
ID e JES 
TM:mroco © JES 


178.5 ± 13.7 ±7.7 
186 ± 10 ±5.7 
172.3j:j^''g^ ± 10.8 
177.1 ±4.9 ±4.7 

171.1 ± 3. 7± 2. it 
174.0 ±2.2 ±4.8 

165.2 ± 4.4 ± 1.9t 
176.9 ± 3.8 ± 1.7t 


[539] 

[537,538] 

[540] 

[570] 

[571,572] 
[278] 
[569]* 
[568]* 



and normalization of the background templates and residual JES uncertainties. The 
compatibility of this measurement with that using the ideogram method [569] on 
an overlapping dataset of almost the same size, which yields a central value ~12 
GeV/c^ lower, is currently being investigated. 
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7.3.4. World-averaged top quark mass 

The top quark mass has been measured at the Tevatron in the three main decay 
channels using various methods, as described in the past sections. In the lepton 
+ jets channel, precisions of 1% are achieved in single measurements, while in the 
dilepton and all-hadronic channels the precision is ~2%. The different methods 
assume standard model tt production and decay, but still exhibit differences in the 
strength of their model dependence. While ME methods provide the best sensitivity, 
they also are strongly model dependent through their implemented matrix elements. 
Template methods purely relying on the measured event kinematics are more robust 
with respect to possible deviations from the standard model, but in general exhibit 
lower sensitivity. 

Pursuing mass measurements in all ti decay channels with different methods is a 
valuable test of the self-consistency of the standard model assumptions, and can also 
be used to probe for new phenomena [296] . While no top quark mass measurement 
has been performed thus far in tt decay modes involving hadronic r decays, the 
progress in the corresponding cross section analyses discussed in Section [5.1. 41 indi- 
cates this could still be possible at the Tevatron, completing the tt decay channels 
for measuring TOj. 

CDF and DO have combined their recent preliminary Run II results with their 
measurements obtained in Run I, respectively. Based on the results of highest sen- 
sitivity in the dilepton, lepton -I- jets and all-hadronic channels (CDF), and in the 
dilepton and lepton -|- jets channels (DO), both experiments attain a total precision 
of 0.9% on their combined measurements, respectively. Based on analyses using up 
to 2.7 ib-\ CDF gets 172.4 ± 1.0(stat.)± 1.3(syst.) GeV/c^ [573], while DO obtains 
172.8 ±0.9(stat.)±1.3(syst.) GeV/c^ [574] using analyses on up to 2.8 fb^i of data. 

The results of both experiments are in very good agreement, and their overall 
combination yields 172.4 ± 0.7(stat.) ± 1.0(syst.) GeV/c^ [12], corresponding to an 
overall precision of 0.7%, as illustrated in Fig. [53] Combining the separate results 
from the all-hadronic, the lepton -I- jets and the dilepton channels, yields 177.5 ± 
4.0 GeV/c2, 172.2 ± 1.2 GeV/c^ and 171.5 ± 2.6 GeV/c^, respectively. These results 
are consistent with each other, and have probabilities of at least 17% between 
any two of the channels. 

All these combinations are calculated using the BLUE method [276,277] and as- 
sume Gaussian systematic uncertainties, with their correlations properly accounted 
for. The different sources of systematic uncertainties are broken down into twelve 
orthogonal categories. Six of them deal with uncertainties related to the JES, while 
others address signal and background modeling, fitting procedures, specifics of MC 
generation and lepton energy scale. The main contributions to the 1.0 GeY/c^ sys- 
tematic error on the world average of are (in units of GeV/c^): total JES (±0.8), 
signal, background and MC model (±0.3 each) and lepton scale and fitting proce- 
dure (±0.1 each). 

Having reached a precision of 0.7%, the world- averaged top quark mass is now 
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Fig. 53: Measurements of the mass of the top quark used as mput to the current 
prehminary world average [12]. 

hmited by the systematic uncertainties that in turn are dominated by JES-related 
uncertainties. Further improvements on the JES can be expected since the increasing 
integrated luminosity will help constrain the corresponding uncertainties better, 
especially the significant contribution from in-situ JES calibration. While a final 
absolute top quark mass uncertainty of Anit < 1 GeV/c^ should be achievable by 
the end of Run II, it will still require a significant effort to determine the contributing 
systematic uncertainties consistently among the experiments, and evaluate any new 
contributions that should be considered at this level of precision. 

This measurement marks the most precise determination of a quark mass and 
will certainly provide a legacy well into the LHC era, where it will serve as an 
important calibration signal until large datasets can produce more refined mea- 
surements. However, improving the precision by another order of magnitude can 
only be expected from a threshold scan of ti production at a future linear e^e~ 
collider [575-577]. 

Before the impact of the current top quark mass measurement is discussed, 
it should be noted that the value of this SM parameter depends on the defining 
convention. For instance, the MS calculation gives a value lower than the pole 
mass of the top quark propagator by wlO GeV/c^ at 0{al). The pole mass itself 
exhibits an intrinsic ambiguity of ©(Aqcd) ~ 0.2 GeV (see for example Ref. [578] 
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Fig. 54: Mass dependences of the tt production cross sections measured by DO in the 
lepton + jets channel (left, [262]), and in a combination of lepton + jets, dilepton 
and r + lepton channels (right, [257]), compared with the theoretical prediction 
by Moch and Uwer [127], based on CTEQ6.6M [126] PDFs. The previous world- 
averaged mt — 172.6±1.4 GeV/c^ [581], as well as the 68% C.L. contours of the joint 
likelihoods resulting from the convolutions of measurement and prediction [579] , are 
also shown. 

and references therein). 

The top quark mass measurements described in this review are usually inter- 
preted as representing the pole mass. However, they are calibrated using LO MC 
simulations with higher orders approximated by parton showers, where the top mass 
parameter does not follow a theoretically well-defined convention. Hence, calcula- 
tions and predictions using the measured mass as the pole mass should be taken 
with a grain of salt. 

DO has conducted consistency checks of the compatibility of the direct top quark 
mass measurements at the Tevatron with the pole mass extracted from the it pro- 
duction rate [262, 579] . Comparing the measured tt production cross section with SM 
predictions derived at NLO, including soft-gluon resummations that are performed 
in a well-defined renormalization scheme using the top quark pole mass, constrains 
the mass of the top quark. The cross section measurements depend less on MC mod- 
eling of signal kinematics than direct mass measurements do. For cross sections, the 
MC is mainly needed for determining the signal acceptance, which is expected to 
be rather insensitive to higher order corrections: A comparison of NLO and LO 
predictions shows that higher order corrections affect more the normalization than 
the shape of the relevant kinematic distributions [580]. 

In a recent analysis, DO uses two tt cross section measurements to extract con- 
straints on the top quark mass by comparing measurement with theoretical pre- 
dictions [579]. One result is obtained in the lepton + jets channel based on the 
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combination of a counting experiment using b tagging and an analysis utilizing a 
topological multivariate discriminant for 0.9 fb~^ of data [262], as discussed in Sec- 
tion [51121 The second result combines measurements in the lepton + jets, dilepton 
and T + lepton channels obtained from approximately 1.0 fb"^ of data [257] (see 
Section [5. 1.5p . Being based on different analysis techniques and different final states, 
both results exhibit a different dependence on the top quark mass, as illustrated in 
Fig. [Ml 

One theoretical prediction used for the comparison is that of Moch and 
Uwer [127], and is based on CTEQ6.6M [126] PDFs (see Section [SI]). A joint 
likelihood depending on mt and is obtained as the product of the likelihood 
functions of the measurement, including its total experimental uncertainty, and the 
theoretical prediction, including scale and PDF uncertainties. The contour of the 
joint likelihood's smallest region containing 68% of its integral is also shown in Fig. 
[54lfor both measurements. By integrating over the tt production rate, the top quark 
mass can be extracted. For the lepton + jets channel measurement a top quark mass 
of 171.2^0 2 GeV/c^ is obtained. The combined lepton + jets, dilepton and r + lep- 
ton measurement yields 169.6^5 5 GeV/c^, which includes an additional systematic 
uncertainty of 1 GeV/c^ due to a smaller mass dependence range available from 
this measurement. Both results are in good agreement with the world-averaged mt 
obtained from the complementary direct measurements. 

The current world-averaged nit of 172.4 ± 1.2 GeV/c^ [12] is also in good 
agreement with its predicted value of 179lg^ GeV/c^ [14] in the framework of 
the standard model, based on precision electroweak data, as discussed in Sec- 
tion 12.31 The top quark mass measurement together with that of the W boson 
mass (80.398 ± 0.025 GeV/c^ [29]) can be used to obtain limits on the Higgs boson 
mass via the radiative corrections on the W boson mass in a global electroweak fit. 
This yields uih = 84^26 GeV/c^ [14], as illustrated in Fig. [55l where the uncertain- 
ties are only from experiment. To demonstrate the impact of the improvements in 
precision for the measurements of both mt and mw since the beginning of Run II, 
the corresponding fit results in spring 2004 [582] are also shown in Fig. \S5[ The 
current resulting 95% C.L. upper limit on the Higgs boson mass is 154 GeV/c^, 
which includes both experimental and theoretical uncertainties. 

The direct searches for the standard model Higgs boson at LEP provide a 95% 
C.L. lower bound of 114.4 GeV/c^ [583], as also illustrated in Fig.[55l CDF and DO 
have recently excluded a SM Higgs boson of 170 GeV/c^ mass at 95% C.L., and a 
mass range of about 165 to 175 GeV/c^ at 90% C.L. [584] based on 3 fb^^ of data. 
This is not reflected in Fig. [55l 

8. Summary 

More than thirteen years after its discovery, the properties of the top quark are 
being studied at the Tevatron with unprecedented precision by the CDF and DO 
collaborations. The Tevatron is operating very smoothly, and has already delivered 
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Fig. 55: Mass constraints on the Higgs boson in 2004 (top, [582]), using the Tevatron 
Run I mass combination of mt — 178.0 ±4.3 GeV/c^ [158], and now (bottom, [14]), 
using the current prehminary world-averaged result rrit =172.4 ± 1.2 GeV/c^ [12]. 

more than 4 fb~^ of integrated luminosity to each experiment. Before the end of 
Run II, it is expected to deliver 2 fb~^ per year. CDF and DO have exploited 
these increasing luminosities together with novel advanced analysis techniques to 
improve upon previous measurements, but also to explore top quark properties 
that were not accessible before such as its electroweak production, which provides 
first direct measurements of the |Vtfc| CKM matrix element. Thus far all results are 
consistent with standard model expectations and between the experiments, which 
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constrains specific extensions of the standard model impacting the properties under 
consideration. 

Top quark pair production is well established in the lepton + jets, dilepton and 
all-hadronic final states, and these channels are also used to study further properties 
of the top quark, such as its mass. The top quark signal is also being established 
in final states involving hadronically decaying r leptons. The observed production 
rates in all final states are consistent with each other and with the standard model 
expectation. No novel contributions to ti production and signal samples have yet 
been observed, and corresponding constraints are derived on tt production via reso- 
nances or massive gluons, and possible contributions of a fourth fermion generation 
or scalar top quarks to the selected signal samples. The kinematics of the observed 
events also agree with the SM prediction. Since contributions from physics beyond 
the standard model could affect the observed ti final states differently also via the 
top quark decay, for example, as a result of top quark decays into charged Higgs 
bosons, corresponding limits arc derived as well. 

The ti production cross section has been measured to be «7.3 pb for rrit = 
175 GeV/c^ with a precision of 10%, matching the uncertainties of the theoretical 
predictions. This provides stringent tests of the corresponding perturbative QCD 
calculations. Ultimately, the precision of the cross section at the Tcvatron might 
reach the 6% level, dominated by the uncertainty on the integrated luminosity. First 
measurements have been performed to determine contributions of qq annihilation 
and ghion-giuon fusion to ti production, and arc fomid to be consistent with QCD 
predictions. Also, higher order effects such as the top quark charge asymmetry, 
measured for the first time, agree with SM expectation within statistical precision. 

First evidence for electroweak single top quark production has been found by 
both experiments, and observation at the five standard deviation level appears im- 
minent. The observed production rates are consistent with the standard model 
expectation of f«3 pb and provide first direct measurements of the CKM matrix 
clement \Vtb\- The most stringent 95% C.L. lower limit is found to be \Vtb\ > 0.71. 
Searches for contributions to single top production via mechanisms beyond the SM, 
e.g., mediated by W or charged Higgs bosons, or flavor changing neutral interac- 
tions with gluons, are also performed. The lack of any significant deviations from 
the standard model has been used to set corresponding stringent limits on such 
processes. Sensitivity of single top production to the Wtb vertex structure has pro- 
vided constraints on more general Wtb interactions, including first limits on left- 
and right-handed tensor couplings. 

The decay properties of top quarks have been studied using fl samples provid- 
ing both sufficient statistics and sample purity. The W boson helicity in tt decays 
can now be measured in a model independent way by extracting the fractions of 
left-handed and longitudinally polarized W bosons simultaneously. This provides 
additional information about the Wtb vertex structure that can be used together 
with constraints from single top quark production to fully specify Wtb couplings. 
The W boson helicity measurements are found to be consistent with SM expecta- 
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tion, and further studies in the lepton + jets and dilepton final states wih scrutinize 
the observed > 2 sd discrepancies between the two channels. The ratio of decays of 
t to 6 quarks versus any quarks, i.e., R = B{t Wb)/B{t Wq) = 0.97lo;o8, has 
reached a precision of 9%, confirming the expectation of dominant t Wb decay 
assumed in most analyses. Top quark decays beyond the standard model, mediated 
through neutral currents, invisible decay modes, or decays to charged Higgs bosons 
have been sought, but have yielded only upper limits on such processes. 

Measurements of fundamental top quark properties such as its charge, lifetime 
and mass based on ti final states thus far confirm the standard model nature of 
the top quark. First measurements of charge and lifetime are consistent with expec- 
tation, and help to constrain ideas beyond the SM. The top quark mass has been 
measured in lepton + jets, dilepton and all-hadronic final states, yielding consistent 
results among the channels and between the CDF and DO experiments. Combining 
the results yields rrit = 172.4 ± 1.2 GeV/c^, which marks the most precise measure- 
ment of the mass of any quark with a precision of 0.7%. By the end of Rim II, a 
measurement with an absolute precision of < 1 GeV/c^ should be achievable. This 
result will provide an important calibration at the LHC, until there is sufficient 
luminosity for further refined measurements. 

While these mass measurements are usually interpreted as representing the polo 
mass of the top quark, it should be recognized that their calibration through current 
Monte Carlo simulations raises certain ambiguities of interpretation. Nevertheless, 
indirect mass measurements utilizing the mass dependence of tt production based 
on the pole mass are consistent with the direct measurements. 

Utilizing radiative corrections to the W boson mass in a global electroweak fit 
to data that includes the world-averaged top quark and W boson masses as inputs, 
the mass of the yet to be observed standard model Higgs boson can be constrained. 
This provides a 95% C.L. upper limit on the Higgs boson mass of 154 GeV/c^. 

The large mass of the top quark does not only render it an ideal window to new 
physics, but this most massive known fundamental object also has a lifetime that is 
short compared to hadronization times. Consequently, observables sensitive to the 
top quark spin can be accessed undisturbed by hadronization processes. Spin-related 
measurements have not as yet been performed in Run II, but it was shown in Run I 
that a measurement of tt spin correlations is feasible, and this measurement will 
greatly benefit from the increased statistics. The observation of single top quark 
production in addition should enable first studies of the polarization of top quarks 
when produced via the electroweak interaction. 

For precision measurements such as those of top quark mass and pair prodiiction 
cross section, the study of systematic uncertainties (consistently across experiments) 
and evaluation of any possible new or smaller contributions not considered in the 
past become a high priority. Other measurements, particularly involving single top 
quark production, will remain statistically limited throughout Run II. The LHC 
will be a "top factory" , producing millions of top quarks per year thanks to the two 
orders of magnitude increased production cross sections and enhanced luminosity 
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relative to the Tevatron. A broad top quark physics program is in preparation at 
the LHC [403, 578] that will complement and further expand that of the Tevatron. 

A great many analyses are being pursued at the Tevatron, characterizing top 
quark data samples both as signal and as background contributions for other pos- 
sible processes of similar signature that still remain to be studied. The top quark 
serves as a probe into new physics, both in production and decay, that could appear 
in the form of new particles or as modified couplings relative to the standard model. 
While thus far all measurements are in agreement with the standard model, there is 
still much room for new physics to be explored both at the Tevatron and the LHC. 
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